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ABSTRACT 

An  interactive  computerized  experimental  procedure  is 
presented  for  obtaining  the  elastic-plastic  fracture 
toughneaa  <Jic>  for  a  material.  The  proceaa  employs  a  bend 
type,  single  specimen  approach  by  using  the  measured 
unloading  crack  opening  displacement  (COD)  compliance  to 
obtain  the  crack  length.  Data  are  taken,  manipulated  and 
plotted  during  the  test  to  provide  the  operator  with  the 
necessary  information  to  control  the  testing  machine.  The 
program  automatically  plots  the  J  versus  crack  extension 
curve  and  determines  Jxc  upon  completion  of  the  test.  A 
second  program  ia  presented  to  assist  the  operator  in 
fatigue  precracking  the  apecimen.  A  third  program  is  used 
to  recall,  liat  and  plot  the  data  that  were  stored  on 
diskette  during  the  test. 
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I .   INTRODUCTION 

A.   GENERAL 

The  toughness  of  a  material  ia  a  Measure  of  the  energy 
per  unit  volume  required  to  cause  fracture.  It  can  be 
quantified  by  simply  taking  the  area  under  the  true  stress 
strain  curve  from  a  simple  tension  test  CRef .  11 .  The 
practical  usefulness  of  this  measurement  is  limited  however, 
because  it  provides  no  information  on  the  material's  ability 
to  resist  crack  propagation.  Since  all  real  structures 
contain  notches  or  cracks  of  some  sort,  by  design  or 
otherwise ,  which  act  aa  stress  concentrators,  another 
measurement  of  toughness  was  developed  called  fracture 
toughness. 

The  fracture  toughness  <KC)  of  a  material  ia  a  measure 
of  the  resistance  of  the  material  to  fracture  in  the 
presence  of  a  defect.  K  ia  defined  aa  the  plsne  strain 
fracture  toughness  and  is  used  to  quantify  a  material's 
resistance  to  crack  propagation  CRef.  2].  The  American 
Society  for  Teating  and  Mater iala  (ASTM)  developed 
Standard  E399-83  CRef.  33  which  eatablishes  precise 
detaila  for  accepted  testing  procedures.  The  standard 
also  provides  rigid  criteria  for  the  apecimen  aize  and 
geometry  to  enaure  plane  atrain  conditions  exist  at  the 
tip  of  the  crack. 


The  J- integral  method,  developed  by  Rice  In  1968, 
provldea  a  meana  of  including  plaatic  deformation  at  the 
crack  tip  CRef.  43.  J  ia  a  measure  of  the  intensity  of  the 
entire  elastic-plastic  atreaa-atrain  field  that  aurrounda 
the  crack  tip.  Jjc  ia  defined  aa  the  value  of  J  at  the 
onaet  of  cracking.  J- integral  teating  haa  gained  acceptance 
aa  a  valuable  means  of  evaluating  a  metallic  material 'a 
resistance  to  crack  initiation  aa  demonstrated  by 
establishment  of  ASTH  Standard  E813-81  CRef.  53.  The 
measured  value  of  Jjc  haa  practical  applicationa  in 
measuring  the  effects  of  metallurgical  variables,  heat 
treatments,  and  weldmenta  and  ia  being  uaed  by  researchers 
today  to  make  auch  quantitative  distinctiona  CRef.  5-7] . 

B.   Jic  VERSUS  Kic  TESTING 

The  procedurea  eatabliahed  by  ASTM  Standard  E399-83  for 
obtaining  Kic  impoae  minimum  specimen  aize  requirements 
that  cauae  difficulty  when  teating  materials  that  poaaeaa 
high  toughneaa.  The  apecimen  thickneaa  muat  be  sufficient 
to  cauae  plane  atrain  conditiona  to  exist  CRef.  33. 
Hertzberg  CRef.  23  preaenta  the  empirical  equation 
developed  by  Brown  and  Srawley  CRef.  83 : 

t  and  a  >  2.5  «  <Kj;c/Ya> 
where 

t  *  apecimen  thickneaa 

a  3  crack  length 
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Kic  3  plane  strain  fractura  toughness 

Ya  3  yield  strength 
A  similar  analyaia  for  J-integral  teat  specimens  reveals 
that: 

b  or  B  >  25  •  (Jic/Fa> 
where 

b  =  initial  uncracked  ligament 

B  =  specimen  thickness 

JXc  3  anticipated  value  if  Jxc 

recheck  when  teat  is  completed 

Fa  ■  flow  atreaa  C(Ya+Uta) /23 

Hertzberg  takea  thia  reault  one  atep  further  by  showing  that 
JXc  equala  (Ktc)2/e  (where  E  equala  the  material's  elaatic 
modulua)  for  Jxc-  BY  dividing  the  two  equationa  one  aeea 
that: 

CKxc  apecimen  aize)  >  20  •  (JXc  specimen  alze) 
thus,  becauae  the  J-integral  allows  plaatic  deformation,  one 
of  ita  moat  aignificant  advantagea  ia  that  a  aatlafactory 
apecimen  can  be  obtained  at  one- twentieth  (1/20)  of  the 
thickneaa  of  a  Kxc  apecimen.  Thia  not  only  meana  smaller 
specimens  may  be  uaed  but,  also,  smaller  teating  machinee 
are  required  and  therefore  coata  are  reduced.  CRef .  23 

Although   Jxc   and  Kic  provide  a   different   meaaure   of 
toughness,  Jjc  can  be  uaed  to  eatimate  Kjc.    ASTM   Standard 
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£813-81    CRef.    53   provides   limited   guidelines   on   the 
relationship  between  the  two. 

C.   PURPOSE  AND  SCOPE  OF  RESEARCH 

The  purpose  of  this  thesis  is  to  develop  the  procedures 
required  to  obtain  Jjc  from  J-integral  test  results  using 
equipment  currently  available  within  the  Materials  Science 
group  of  the  Mechanical  Engineering  Department  at  the  Naval 
Postgraduate  School  (NPS) .  A  list  of  the  »a]or  equipment 
follows: 

a.  Series  810  MTS  machine 

b.  HP-9826  computer  with  internal  disk  drive 
and  320K  random  access  memory 

c.  HP-3497A  data  acquisition  system 

d.  HP-3437A  voltmeter 

e.  HP-2671G  thermal  printer 

f .  HP-7225B  graphics  plotter 
Figure  1  shows  the  system  layout .  > 

The  test  method  developed  is  to  comply  with  ASTM 
Standard  E813-81  CRef.  5],  be  easy  for  others  to  use,  and 
require  a  minimum  of  operator  interaction.  It  is 
specifically  designed  for  use  with  three  point  bend 
specimens  using  a  single  vice  multiple  specimen  approach. 
The  method  developed  and  presented  here  covers  the  entire 
test,  including  specimen  preparation,  precracking,  testing, 
and  validation  of  results. 
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II.   NATURE  OF  THE  PROBLEM 

A.   J- INTEGRAL 

The  J- integral  is  a  mathematical  expreaaion  used  to 
characterize  the  near  tip  stress-strain  environment  of 
cracked  elastic-plastic  bodiea  CRef .  9] .  It  la  a  path 
Independent  line  integral  that  aurrounda  the  crack  tip  for 
non- linear  plastic  metallic  mater iala  CRef. 53 .  The  equation 
delveloped  by  Rice  ia  CRef.  43 : 

J  *  [   <W  dy  -  T'Ou/ax)  da> 
where 

W  *  loading  work  per  unit  volume,  or  for  elaatic 

bodiea,  atrain  energy 
I""  ■  the  path  of  the  integral,  which  encloaea  the 

crack  tip 
da  *  increment  of  the  contour  path 
T  *  outward  traction  vector  on  da 
u  =  diaplacement  vector  at  da 
x,y,z  =  rectangular  coordinatea 

T'(du/dx)da  ■  the  rate  of  work  input  from  the  atreaa 
field  into  the  area  encloaed  by  r  . 
Figure    2   providea   a   pictorial   repreaentation   of   thia 
equation. 

The  mathematical  expreaaion  developed  by  Rice  CRef.  43 
haa   been   manipulated   and   reduced  to   aiapler   forma   for 
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specific  specimen  geometries.  Rice,  Paris,  and  Merkei  CRef . 
10]  reduced  the  J- integral  expression  for  a  three  point  bend 
specimen  to: 

2  •  A 


j  = 

b  »  B 

where 

J  =  resistance  to  crack  initiation 
A  =  area  under  load  versus  load  line 

displacement  curve 
b  =  uncracked  ligament 
B  =  specimen  thickness 
This  thesis  is  not  intended  to  redevelop  the  J- integral   but 
rather   apply   it,   therefore  only  the  result   is   presented 
here. 

More  recent  work  by  Ernst,   Paris,   and  Landes  CRef.  11] 
corrects   the  calculation  of  J  due  to  crack  growth  by   means 
of  the  following  equation: 

J(i*l)  *  < Ji*<n/b) i*Ai, i+i>  ■  <1- <Y/b> « <ai+i-ai> ) 
where 

J<i*l)  =  current  value  of  J 
Ji  =  last  J  computed 

n  *  1  for  three  point  bend  type  specimens 
Ai,i*i  =  area  enclosed  by  the  actual  test 
record  and  lines  of  constant 
displacement  di  and  d^+i 
Y    ■  1  for  three  point  bend  type  specimens 
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b  *  current  uncracked  ligament 
W  =  specimen  width 
ai  =  laat  crack  length 
ai*l  a  current  crack  length 
Correcting   for   crack  growth  during  the  teat   la   necessary 
when   uaing   a   aingle  specimen  approach   aa   uaed   by   this 
method . 

The  blunting  line  defined  aa: 

J  a  2  »  Fa  »  delta_a 
where 

Fa  =  flow  atreaa 

delta,  a  =  change  in  crack  length 
approximates  apparent  crack  advance  due  to  crack  tip 
blunting  in  the  abaence  of  slow- stable  crack  tearing.  It  is 
plotted  on  the  J  veraua  crack  extenaion  plot.  A  linear 
regression  fit  through  the  applicable  data  points  is  then 
plotted.  The  intersection  of  the  two  lines  is  by  definition 
Jjc  (see  Figure  3>.CRef.  SI 

B.   PREVIOUS  WORK 

The  early  experimental  methoda  uaed  to  obtain  Jjc 
required  the  preparation  of  aeveral  apecimena.  Each 
specimen  received  a  different  load.  The  specimens  were  then 
heat  tinted  (to  identify  crack  length)  and  broken  apart. 
Each  specimen  provided  one  data  point  on  the  J  versus  crack 
extension  curve,  commonly  referred  to  as  the  J-R  curve  CRef . 
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9  and  113 .  Checks  were  then  performed  on  the  validity  of 
the  data.  A  minimum  of  four  valid  data  points  were  required 
to  obtain  a  linear  least  squares  fit  and  hence  Jic  [Ref. 
5,10  and  12] . 

Researchers  over  the  last  five  years  have  developed  a 
much  improved  test  method  whereby  a  single  specimen  can  be 
used.  Today's  computerized  data  acquisition  equipment  and 
data  manipulation  are  ideally  suited  for  the  task.  By 
relating  the  compliance  (how  much  a  specimen  bends  under  a 
given  load)  to  the  crack  length,  all  the  data  required  to 
obtain  the  J-R  curve  can  be  collected  by  using  a  single 
specimen  CRef. 11-16] .  Clarke  [Ref. 16]  states  that  prior  to 
the  use  of  computers  the  time  required  for  data  reduction  to 
obtain  the  J-R  curve  was  four  to  five  hours.  Using 
computerized,  single  specimen,  unloading  compliance  methods 
it  is  now  possible  to  to  obtain  the  J-R  curve  in  a  matter  of 
twenty  to  thirty  minutes  after  the  test  starts.  The  obvious 
savings  in  specimen  preparation  and  data  reduction  have  made 
the  J-integral  a  valuable  means  of  analysis  where  no 
practical  means  previously  existed. 

C.   TEST  CRITERIA 

The  test  method  developed  in  this  thesis  is  specifically 
designed  for  use  on  a  three  point  bend  specimen.  It 
utilizes  the  single  specimen  approach  by  relating  the  crack 
opening   displacement  (COD)  compliance  to  the  crack   length. 
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The  following  equations  are  currently  uaed  by  DTNSRDC.  For 
bend  specimens  with  a  apan  to  width  (S/W)  ratio  of  4.0,  Dr. 
J.  Joyce  CRef .  173  of  the  United  States  Naval  Academy 
determined: 

a/W  *  0.998265  -  3.81662  »  U  -  1.80596  »  u2  ♦  32.31038  »  U3 
-  44.15665  »  U4  -  52.67876  «  U5 


where 


a  *  crack  length 
W  =  specimen  width 


U  * 


/4  »  Be  »  E  •    W  «  compl  U/2  ♦  1 


where 

Be  *  effective  thickness  =  Bmax  -  (Bmax  -  Bmin)2/Bmax 

E  =  material's  elastic  modulus 

S  =  span  between  supports 

W  =  specimen  width 

compl  *  compliance  =  1/alope  of  the  load  versus  COD   line 
for  the  specific  unloading. 
The  reason   for  using  the  single  specimen   approach   has 
already   been  presented.    The  reasons  a  bend  type   specimen 
and  COD  type  unloading  compliance  are  uaed  are: 

a.  Three   point  bend  specimen  geometry  is  better   suited 
than  compact  tension  specimens  for  testing  weldments. 

b.  Three   point  bend  specimen  geometry  is  also  easier  to 
machine  than  compact-tension  specimens. 
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c.  COD  voltage  readings  are  leas  susceptible  to  noise 
than  linear  variable  differential  transformer  (LVDT) 
voltages  used  to  measure  load  line  displacement. 

d.  The  COD  compliance  equation  affords  0.1'i  fitting 
accuracy  and  allows  comparison  of  results  with 
DTNSRDC  test  data  CRef.  171. 
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III.   TEST  METHOD 

A.   OVERVIEW 

A  J-integral  teat,  can  be  divided  into  three  diatinct 
steps,  1)  pre-test,  2)  test,  and  3)  poat-teat.  Each  of 
theae  are  deacribed  in  greater  detail  in  thia  chapter  and 
preaented  aa  a  atep  by  atep  procedure  in  Appendix  A. 

After  the  specimen  la  Machined  it  must  be  precracked 
uaing  cyclic  fatigue  loading  to  minimize  reaidual  atresses 
and  create  a  aharp  crack  tip  CRef .  53 .  A  computer  program 
was  developed  to  asaiat  the  operator  in  accomplishing  this. 
The  program  computes  the  loading  parameters  for  precracking 
in  accordance  with  ASTM  Standarda  E399-S3  and  ASTM  E813-S1 
CRef.  3  and  53 .  The  operator  uaea  thia  information  to  run 
the  MTS  machine. 

During  the  J-integral  teat  a  aecond  computer  program 
collects  and  recorda  the  load,  COD,  and  load  point 
displacement  (deflection) .  Theae  data  are  used  by  the 
program  to  determine  J  and  the  amount  of  crack  extenaion  for 
each  operator  controlled  unloading  cycle.  An  unload  cycle 
conaiata  of  decreasing  the  load  by  10*  of  the  expected 
maximum  load  and  reloading  to  the  original  load  point.  The 
program  providea  the  operator  with  a  real-time  graphical 
repreaentation  of  the  any  two  parametera  (load  versus  COD  or 
deflection)   during   the  teat.    It  will  alao  plot   the   J-R 
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curve  and  compute  Jjc  upon  completion  of  taking  all  the  teat 
data.  Data  are  atored  on  diskette  to  serve  na  a  record  and 
allow  future  recall  listing  or  plotting. 

The  post-teat  process  includes  heat  tinting  the  specimen 
to  aid  in  visual  verification  of  crack  extension  upon 
completion  of  the  test  as  required  by  ASTM  Standard  ES13-81 
CRef .  53 . 

B.   SPECIMEN  PREPARATION 

Figure  4  shows  specifically  the  dimensions  used  in 
preparing  eight  apecimens  from  two  different  one  inch  thick 
welded  plates  of  HY-80.  One  plate  was  preheated  while  the 
other  was  not.  It  was  necessary  to  reduce  the  width  to 
approximately  0.95  inches  to  remove  the  distortion  (bow) 
created  by  the  weld.  Four  apecimens  were  prepared  from 
each  plate,  the  machined  notch  was  placed  in  the  weld  metal 
for  two  and  in  the  base  metal  for  the  other  two.  The  clip 
gage  (used  to  measure  COD)  is  attached  to  the  top  of  the 
apecimen  by  means  of  two  cut  razor  blades.  The  blades  are 
held  down  by  small  straps  which  allow  adjustment  of  the 
initial  gap.  The  initial  gap  is  set  to  0.475  inches,  the 
calibrated  zero  aetting  for  the  clip  gage,  by  meana  of  a 
machined  gage  tool  that  reats  in  the  apecimen 'a  machined 
notch.     Thia  method  ia  also  uaed  by  DTNSRDC. 

Two  additional  holes  are  drilled  and  tapped  for 
retaining  a  small  bracket  on  the  face  of  the  apecimen.    The 
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bracket  is  used  to  provide  a  sensing  point  for  a  linear 
variable  differential  transformer  (LVDT) .  The  LVDT  signal 
is  used  to  calculate  the  area  under  the  load  versus 
displacement  curve.  At  present  the  external  LVDT  is  unable 
to  control  the  MTS  machine  so  the  current  method  uses  the 
MTS'  internal  LVDT  to  measure  deflection  via  actuator 
movement. 

Clarke  CRef .  16]  states  that  side  grooves  improve  the 
accuracy  of  calculating  the  crack  length  when  using 
unloading  compliance  procedures.  The  side  grooves  reduce 
plastic  deformation  on  the  surface  and  aid  in  propagating  a 
straighter  crack  front.  DTNSRDC  has  found  20X  side 
grooving  produces  the  best  results,  hence  the  same  was  used 
for  this  work. 

After  the  specimen  is  prepared  a  test  record  sheet  (see 
Appendix  B)  is  completed.  Notch  and  side  groove  depths  can 
be  measured  using  a  dial  indicator  with  a  needle  tip.  The 
record  provides  the  operator  with  ready  access  to  the 
information  necessary  to  run  the  precracking  and  test 
programs. 

C.   PRECRACKING 

Before  the  J-integral  test  can  be  performed  the 
specimen  must  be  precracked  to  establish  a  sharp  crack  tip. 
ASTM  Standards  ES13-81  and  E399-83  establish  parameters  that 
must  be  met  during  precracking  to  prevent  over-stressing  the 
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specimen  or  propagating  the  crack  through  the  entire  width. 
Precracking  ia  accomplished  with  the  aid  of  a  program  called 
"PRECRACK"  which  will  compute  loading  requirements  based 
upon  these   paramenters . 

The  program,  written  in  HP-Basic,  affords  flexibilty  by 
allowing  the  use  of  different  bend  specimen  sizes  and 
materials.  It  is  easy  to  use  due  to  its  interactive  design 
and  easy  to  modify  should  such  a  need  arise.  A  listing  of 
"PRECRACK"  is  provided  in  Appendix  D. 

The  operator  is  requested  by  the  program  to  enter 
specific  specimen  data.  The  program  will  compute  the  Initial 
maximum  and  minimum  loads  and  change  in  the  associated 
stress  intensity  factor  (delta  K) .  If  this  delta  K  exceeds 
40  kai'lnl^  then  new  maximum  and  minimum  loads  are 
computed.  The  operator  ia  responsible  for  ensuring  the  NTS 
machine  ia  operating  within  the  prescribed  limits.  The 
operator  can  obtain  a  hard  copy  of  the  entered  and  computed 
data  for  record  purpoaea  (see  Appendix  G) . 

The  NTS  machine  produces  cyclic  loading  by  means  of  an 
inverted  (compressive  load)  haver  sine  wave. 

The  operator  must  periodically  interrupt  the  cyclic 
loading  and  measure  the  current  crack  length  by  manually 
loading  and  unloading  the  specimen  (determine  the  specimen's 
compliance) .  The  program  computes  the  new  crack  length  and 
new  loading  limits.  During  the  last  SOX  of  precracking: 

Pmax  <  .4Pl 
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where 

Piax  =  maximum  load 

PL  -  C(4/3)  »  (B  »  b2  «  Fa  /  (2  »  Ha))] 
where 

B  ■  net  thickness  (Bmin  if  apecimen  haa  aide  groovea) 

b  =  current  uncracked  ligament  (W-a) 

Fa  =  flow  atreaa  £ CYa+Uta) /2) 

Ha  =  half  apan 
and  Delta_K/E  <  0.005  mml/2  (0.001  inl'2>  whichever, 
producea  the  leaat  load  CRef.  5] .  The  program  will 
automatically  begin  uaing  theae  parametera  when  the  moat 
recently  calculated  crack  length  exceeda  thia  50*  of  the 
required  extenaion  to  reach  a/W  =  0.6.  Thia  value  ia  uaed 
becauae  it  ia  the  optimum  crack  length  with  which  to  begin 
the  J- integral  teat  CRef.  5] . 

D.   J- INTEGRAL  TEST 

The  J- integral  teat  ia  performed  by  slowly  applying  a 
load  to  the  apecimen  and  periodically  unloading  it  by  10*  of 
the  expected  maximum  load.  The  NTS  machine  ia  set  up  in 
atroke  control  with  a  very  amall  inverted  ramp  function 
reaulting  in  a  loading  rate  of  0.005  inch/minute.  The 
operator  muat  again  interrupt  the  MTS  program  and  manually 
unload  the  apecimen.  Each  unloading  cycle  (unload  and 
reload)  providea  the  neceaaary  information  to  compute  the 
current    crack   length.   The  crack  length  la  then   uaed   to 
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determine  the  value  of  J  up  to  that  point.  The  computer 
program  "J_INTGL"  provides  the  operator  with  pertinent 
information  during  the  teat.  Once  the  teat  ia  completed , 
the  operator  can  obtain  the  J-R  curve  and  value  of  Jjc  (see 
Figure  4  and  Appendix  H) .  The  program  will  uae  only  data 
pointa  between  the  offsets,  (blunting  line  ♦  0.15mm)  to 
(blunting  line  ♦  1.5mm)  aa  required  ASTM  Standard  E813-81 
when  computing  the  linear  regreaaion  line.  A  minimum  of 
four  data  pointa  muat  be  within  the  offaets  and  one  must  lie 
near  (within  1/3  the  overall  range  of  delta_a)  the  blunting 
line.  Specific  acceptability  criteria  for  data  is  provided 
in  ASTM  Standard  E813-81  CRef.  51. 

The  program  "J_INTGL"  ia  the  foundation  of  the  test 
method ,  the  single  specimen  approach  would  not  be  possible 
without  it.  The  program  collects,  manipulates,  displays  and 
recorda  all  pertinent  data  during  the  test.  A  flow  chart  of 
"J_INTGL"  ia  provided  in  appendix  C.  A  Hating  of  MJ_INTGL" 
ia  provided  in  Appendix  E. 

Deflection  measurements  ar&  preaently  taken  uaing  the 
LVDT  installed  in  the  NTS  machine.  A  correction  is 
therefore  required  to  account  for  the  system's  compliance 
while  loading  a  apecimen.  The  correction  waa  obtained  by 
loading  the  system  with  a  two  inch  thick  piece  of  HY130  and 
recording  the  stroke  reading  for  various  loada.  This 
process   waa  checked  three  times  with  consistent  (within  5X) 
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reaulta.    A   plot,  of  the  data  and  the  4th  order   ploynomial 
curve  fit  ia  included  aa  Figure  5. 

E.   POST  TEST  VERIFICATION 

The  reaulta  tust  be  verified  upon  completion  of  the 
teat.  ASTM  Standard  E813-S1  preaenta  the  following  criteria 
that  auat  be  aet  before  the  teat  can  be  conaidered  valid: 

1.  Total  crack  extenaion  (visually  Measured)  must  be 
within  15*  of  that  recorded  during  the  teat. 

2.  bi  and  Bain  >  25  »  Jjc  /  Fa  (all  terns  are  aa 
defined  earlier  except  bi  ■  post -fatigue  uncracked 
ligament) . 

3.  The  alope  of  the  linear  regreaaion  fit  to  the  data 
auat  be  leaa  than  flow  atreaa  (dJ/da  <  Fa). 
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IV.   RESULTS 

A.   TEST  RESULTS 

Specimen  HY80-5B,  prepared  from  the  base  metal  of  the 
preheated  welded  plate,  was  precracked, and  tested  resulting 
in  a  Jic  value  of  142  kJ/m2.  Figure  8  is  a  tracing  of  a 
photograph  of  the  final  fracture  surface  after  heat  tinting. 
All  post  test  checks  were  completed  and  yielded  satisfactory 
results.  The  Jic  value  is  in  good  agreement  with  previous 
tests  of  HY80  plate  as  conducted  by  DTNSRDC  resulting  in  an 
average  Jic  of  144  kJ/m2  CRef .  6] . 

Plots  of  J-R  curve,  load  versus  COD,  load  versus 
deflection  and  the  are  included  as  Figures  5,  6  and  7, 
respectively.  The  load  versus  deflection  curve  was  obtained 
using  the  program  "J_DATA" ,  thus  proving  that  the  data  is 
stored  on  disk,  and  is  able  to  be  recalled. 

Data  was  taken  at  a  nearly  constant  rate  of  2.5  data 
points  per  second  (150  per  minute)  while  the  program  was 
executing. 

Test  reults  were  verified  following  the  guidelines 
provided  by  ASTM  Standard  E813-81 .  The  fracture  surface  was 
enlarged  by  taking  a  photograph,  a  tracing  of  which  is 
included  as  Figure  8.  The  observed  crack  length  was  14. 4* 
longer  than  the  last  recorded  crack  length.  This  is  within 
the   acceptable   limit   (recall  ASTH   requires   <   15*)   but 
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because  the  clip  gage  reached  ita  limit  detector  and  caused 
the  test  MTS  to  stop  between  unloadinga,  it  is  possible  some 
crack  extension  occurred  after  the  last  unloading  cycle. 
This  would  help  to  explain  the  difference. 

Additional  information  on  this  specific  test  is  provided 
by  being  used  in  the  example  calculations  shown  in  Appendix 
A.  The  exemplary  computer  output  for  each  of  the  programs 
is  also  from  specimen  HYS0-5B.   (See  Appendixes  G,  H  and  I). 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

The  method  developed  and  presented  by  this  work  produces 
accurate  and  reliable  results  that  meet  the  specifications 
established  by  ASTM  Standard  ES13-81.  Precracking  specimen 
HY30-5B  took  approximately  2  hours  and  the  J-test  required 
approximately  2.5  hours  including  several  pauses  to  manually 
verify  several  calculations.  The  time  required  for  the 
processes  will  most  likely  decrease  aa  operator  experience 
and  confidence  increases. 

All  three  computer  programs  have  run  successfully  and 
perform  as  intended. 

B.  RECOMMENDATIONS 

As  with  all  first  versions  of  computer  programs  aa 
extensive  as  "J_INTGL",  changes  and  improvements  are  to  be 
expected.  The  following  programming  changes  may  be  of  value 
to  future  operators: 

1.  Display  pertinent  information  (i.e.  current  crack 
length  in  both  English  and  metric  units  (currently 
only  metric  is  provided) . 

2.  Provide  a  means  of  storing  the  data  arrays  during  the 
test   without   being   required  to   rerun   the   entire 
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program  (re-entering  specimen  data,   conversion  data, 
etc. ) 
3.   Provide   a  means  of  eliminating  a  specific   J-R   data 
point    from   being   used   in   obtaining   the   linear 
regression   fit.     This   may   be   programmed   to   be 
accomplished   automatically   by  applying   Chauvenet's 
Criterion  to  the  data. 
Improvements  in  the  current  means  of  conducting  the  test 
can   also   be  made,    such  as  making  the   MTS   control   both 
loading   and   unloading  using  the  same  rate,   requiring   the 
operator   to   do   no  more  than  press  a   button   or   throw   a 
switch.    Lastly,   pursue  the  use  of  the  external  LVDT,  this 
will   eliminate   the   need   for   correcting   the   deflection 
reading  and  speed  up  the  rate  of  data  acquisition. 
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Figure  4.   Specimen  Drawing 
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Figure  8.   Tracing  of  HY80-5B  Fracture  Surface 


37 


APPENDIX  A 
TESTING  METHODS 

A.  INTRODUCTION 

The  steps  required  to  perform  a  J-intagral  teat  uaing  a 
aingle  specimen  approach  can  be  divided  into  three  major 
categories,  pre- test,  test,  and  post  teat.  Each  of  these 
categories  may  be  further  divided  as  shown  below. 

1.  Pre-teat 

a.  Machine  apecimen 

b.  Precrack  specimen  by  fatigue  loading 

2 .  Teat 

a.   Run  J-integral  teat 

3.  Poat-teat 

a.  Heat  tint  apecimen  and  break  apart 

b.  Measure  crack  extension  and  verify  results 
Each   proceaa   ia  described  in  greater   detail   in   this 

appendix.  A  working  knowledge  of  the  equipment  is 
preauppoaed  and  awareneaa  of  all  aafety  precautiona  ia 
required  before  operating  the  MTS . 

B.  PRE-TEST 

a.   Machine  Specimen 

ASTM  Standard  E813-81  providea  very  specific  guidance  on 
apecimen   dimenaiona.    All   dimenaiona  are  related   to   the 
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width  (W) .  Fatigue  precracking  ia  enhanced  if  the  radius  at 
the  tip  of  the  notch  ia  kept  to  approximately  0.003  inchea. 
Thia  waa  accomplished  by  machining  with  a  apecial  cutter 
which  waa  uaed  to  make  only  the  final  paaaea  in  the  notch. 
Figure  2  can  be  conaulted  for  the  actual  dimenaiona  uaed  to 
machine  the  apecimena  uaed  in  thia  work. 

b.   Precracking 

The   following  atep  by  atep  guide  ia  provided  to   assist 
operators . 

1.  Before  turning  on  the  MTS  machine  ensure  the  proper 
range  card  is  installed  in  the  strain  aignal  conditioner  for 
the  clip  gage  to  be  used. 

2.  Turn  on  power  to  all  equipment  with  the  computer  being 
awitched  on  laat . 

3.  Load  and  run  "PRECRACK"  program  on  the  HP -9826  computer. 

4.  Select  range  acalea  on  MTS; 

example:    10  v  ■  4000  lbs  20*  of  full  load  (range  #3) 

10  v  =  0.02  in  10%  of  full  clip  gage  (range  #4) 
10  v  ■  2.5  in  SOX  of  full  stroke  (range  #2) 

5.  Preaa  kO-key  on  computer  and  enter  specimen  data  as 
requeated. 

6.  Preaa  k2-key  and  enter  conversion  data  aa  requeated. 

7.  Preaa  k5-key  to  obtain  a  table  of  COD  readings  for 
various  loada  and  crack  lengtha. 

8.  Set  READOUT  knoba  on  the  MTS  Controller  to  be: 
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top  =  load 
middle  =  strain 
bottom  =  stroke 
Set   all   to   (♦)  polarity.    See  Figure  Al .    This   ensures 
channels   of   the  NTS  machine  are  properly  aligned  with   the 
data  acquisition  system. 

9.  Set  MTS  limit  detector  switches  based  on  the  information 
just  generated  by  the  program.  Limit  switch  settings  are 
determined  by  the  simple  formula: 

Limit  switch  ■  L  »  10  /  FSR 
where 

L  =  desired  limit  in  English  units 
FSR  =  full  scale  range  of  current  operating 
being  used  in  English  units 
example  "PRECRACK"  determined: 

Pmax  =  1668  lbs  set  load  limit  to  1600  lbs 
for  specimen  HY80-5B 
load  limits: 

upper  limit  small  positive  load  for  example  40  lbs 

(should  never  be  positive) 
♦40  lbs  »  10  v  /  4000  lbs  =  ••■0.1  v 
therefore  set  upper  limit  knob  to  +0.1  v 

lower  limit  1600  lbs  maximum  load 

-1600  lbs  »  10  v  /  4000  lbs  =  -4.0  v 
therefore  set  lower  limit  knob  to  -4.0  v 
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strain  limits: 

upper  limit  0.01  in  beyond  maximum  expected  during 

precracking 
■►0.01  in  »  10  v  /  0.02  in  =  *5.0  v 

lower  limit  some  small  value  for  example  -0.1  v 

(should  never  be  negative) 
stroke  limits: 

upper   limit  1.5  in  to  allow  clearance  for   installing 

specimen   and   yet  stop  the  MTS   if   thia 
limit  is  reached 
♦1.5  in  «  10  v  /  2.5  in  =  +6.0  v 

lower   limit  0.25  in  of  compressive  actuator   movement 

is  more  than  sufficient  for  precracking 
-0.25  in  *  10  v  /  2.5  in  =■    -1.0  v 

10.  Adjust  the  upper  roller  supports  as  required  to  obtain  a 
span-to-width  ratio  <S/W)  of  4.0. 

11.  Install  specimen  (with  the  clip  gage  already  installed) 
using  alignment  collars  to  properly  center. 

12.  Zero  the  load  and  strain  voltages. 

13.  Reset  and  clear  all  limit  switches  (all  green  lights  on 
controller  panel  should  be  showing) . 

14.  Using  load  control  energize  hydraulics  and  apply  a  low 
(<100  lbs)  compressive  load  to  the  specimen  by  SLOWLY 
turning  the  SETPOINT  knob  counter-clockwise. 
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15.  Slowly  decraae  the  compressive  load  on  the  specimen  to 
between  5  and  10  lba  by  turning  the  SETPOINT  knob  clockwise 
and  then  zero  the  stroke  voltage  reading. 

YOU  ARE  NOW  READY  TO  OBTAIN  AN  UNLOADING  COD  COMPLIANCE 
MEASUREMENT  OF  THE  CRACK  LENGTH 

16.  Press  k4-key  on  the  computer  and  proceed  with  taking 
data. 

17.  Slowly,  and  as  smoothly  as  possible  increase  the  load  to 
approximately  500  lbs  and  return  to  approximately  20  lba 
ending  the  loading  cycle  with  a  alight  (10  lbs)  increase  to 
trip  the  computer  out  of  the  data-taking  loop. 

The  computer  will  display  all  pertinent  crack  length 
information  and  ask  if  you  want  another  measurement.  You 
may  repeat  this  process  as  often  as  you  like,  just  be  sure 
the  maximun  load  is  lower  than  that  required  to  extend  the 
crack . 

Id.  Make  necessary  adjustments  to  the  MTS  such  that  the 
strain  limit  switch  will  trip  the  MTS  when  the  crack  extends 
to  50*  of  the  desired  extension.  This  is  determined  by 
interpolating  results  from  the  table  obtained  earlier.  (See 
Step  7) . 

NOTE:  As  the  crack  propagates  the  acceptable  maximum 
load  decreases.  For  the  specimen,  HY80-5B,  used  to  test 
this  program  a  Pmax  ■  1000  lba  was  auitable  to  initiate  the 
crack  and  propagate  it  to  50%  of  the  desired  extension 
(as  0.506  in) . 
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Example: 

50*  crack  extension  *  crack  length  (a)  of  0.506  in 
from  table: 

COD  ■  1.7  volts  at  1000  lbs  and  a  a  0.506  in 
Therefore  set  UPPER  STRAIN  LIMIT  SWITCH  to  +  1.7  v. 
YOU  ARE  NOW  READY  TO  BEGIN  CYCLIC  LOADING 

19.  Setup  NTS: 
DISPLAY: 

a.   Setup  to  read  maximum  and  minimum  load  in  lbs. 
FUNCTION  GENERATOR: 

a.  CONTROL  MODE:  Remote 

b.  RATE  1:  10  Hz 

c.  RATE  2:  N/A 

d.  HOLD  BUTTON:  In 

e.  BREAKPOINT:  N/A 

f .  HAVER  SINE:  In 

g.  INVERT:  In 

CONTROLLER  : 

a.  SPAN  l:  0.0 

b.  SPAN  2:  O.O 

COUNTER: 

a.  COUNTER  INPUT:  Oscillator 

b.  COUNTER  MULTIPLIER:  X10O 

c.  Press  RESET 

20.  Slowly  increase  the  SETPOINT  to  Pmin. 

21.  Press  the   RUN  button  on  the  CONTROLLER  and  then  release 
the  HOLD  button  of  the  FUNCTION  GENERATOR. 

22.  Slowly   increase  the   span   of  the  cyclic  loading   using 
SPAN  1  knob  to  the  desired  Pmax. 
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The  MTS  will  continue  cyclic  loading  until  it  is 
manually  stopped  or  a  limit  detector  ia  tripped.  Most 
likely  the  COD  limit  detector  will  trip  at  which  time  the 
machine  will  will  atop  and  unload  the  specimen. 

23.  A  new  crack  length  measurement  can  now  be  obtained 
using  the  same  manual  process  as  described  in  Step  16  and 
Step  17. 

24.  Readjust  limit  switches  and  resume  fatigue  loading  as 
done  in  Steps  18  through  22. 

25.  Continue  this  process  until  the  crack  reaches  the 
desired  length. 

C.   J- INTEGRAL  TEST 

The  following  steps  are  required  to  run  a  J- integral 
test.  They  are  very  similar  to  those  used  in  precracking. 
The  major  change  is  that  a  ramp  function  under  stroke 
control  is  now  used  to  control  the  MTS  machine. 

1.  Before  turning  on  the  MTS  machine  ensure  the  proper 
range  card  ia  installed  in  the  strain  signal  conditioner  for 
the  clip  gage  to  be  used. 

2.  Turn  on  power  to  all  equipment  with  the  computer  being 
switched  on  last . 

3.  Load  and  run  "J_INTGL"  program  on  the  HP-9826  computer. 

4.  Select  range  scales  on  MTS; 

example:     10  v  =  4000  lbs  20*  of  full  load  (range  #3) 
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10  v  =  0.02  in  10%  of  full  clipgage  (range  #4) 
10  v  ■  0.5  in  10*  of  full  stroke  (range  #2) 

5.  Press  kO-key  on  computer  and  enter  specimen  data  as 
requested . 

6.  Press  kl-key  and  enter  conversion  data  as  requested. 

7.  Press  k3-*key  and  follow  instructions  for  creating  a  data 
file  on  an  initialized  data  diskette.  See  HP-Mannual  on  how 
to  initialize  a  diskette  if  necessary. 

8.  Press  k3-key  to  setup  the  plotter. 

9.  Press  k5-key  and  obtain  the  maximum  load  you  may  apply 
while  measuring  aO. 

YOU  ARE  NOW  READY  TO  SETUP  THE  MTS  MACHINE 

10.  Set  READOUT  knobs  on  the  NTS  Controller  to  be: 

top  =  load 
middle  =  strain 
bottom  ■  stroke 
Set  all  to  (O  polarity.   See  Figure  Al . 

11.  Set  NTS  limit  detectors  switches, 
example: 

Pmax  =  S95  lbs  for  HYS0-5B  (see  Appendix  H) 
load  limits: 

upper  limit  small  positive  load  (40  lbs) 

(should  never  be  positive) 
>40  lbs  »  10  v  /  4000  lbs  =  +0.1  v 

lower  limit  895  lbs  compression  maximum 

-895  lbs  *  10  v  /  4000  lbs  =  -2.125  v 
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strain  limits: 

upper  limit. want  to  be  able  to  read  full  scale 

on  range  #2   (0.1  in) 

*0.1in  »  10  v  /  0.1  in  a  +10  v 

lower  limit  a  very  snail  value  set  to  -0.1  v 

C  should  never  be  negative) 

stroke  limits: 

upper   limit  full   scale  (+10.0  v)   on   this   most 

sensitive  scale 

lower   limit  a  maximum   compression  of  0.25   in   is 

sufficient  for  the  test  and  will   stop 
the   actuator  before  it  can  cause   any 
serious  damage  to  the  specimen  holder 
-0.25  in  »  10  v  /  0.5  in  »  -5.0  v 

12.  Adjust  the  upper  roller  supports  ab  required  to  obtain  a 
span-to-width  ratio  of  4.0. 

13.  Install  the  specimen  (with  the  clip  gage  already 
installed)  using  the  alignment  collars  to  properly  center. 

14.  Zero  the  load  and  strain  voltages. 

15.  Reset  and  clear  all  limit  switched  (all  green  lights  on 
controller  panel  should  be  showing) . 

16.  Using  load  control  energize  hydraulics  and  apply  a  low 
(<100  lbs)  compressive  load  to  the  specimen  by  SLOWLY 
turning  the  SETP0INT  knob  counter-clockwise. 

17.  Decrease  the  load  to  between  5  and  10  lbs  and  zero  the 
stroke  voltage. 
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18.  Adjust  the  MTS  to  the  following  settings: 
DISPLAY: 

8.   Setup  to  resd  loed  in  lbs  snd  stroke  in  inches. 
FUNCTION  GENERATOR: 

s.   CONTROL  MODE:  Remote 

b.  RATE  l:  6000  (sec) 

ss   determined  by:   desire  0.005  in/min  loading  rate 
therefore; 

0.5  in  (full  scale)  /  0.005  in/mis  •  60  sec/min  =  6000  sec 

c.  RATE  2:  N/A 

d.  HOLD  BUTTON:  In 

e.  BREAKPOINT:  N/A 

f.  RAMP  BUTTON:  In 

g.  INVERT:  IN 

CONTROLLER : 

a.  SPAN  l:  10.0  (full  scale) 

b.  SPAN  2:  O.O 

COUNTER:  N/A 

YOU  ARE  NOW  READY  TO  MEASURE  THE  INITIAL  CRACK  LENGTH 

19.  The  computer  will  begin  taking  data  after  the  ^CONTINUE' 
key  is  pressed.  The  MTS  will  begin  loading  the  specimen 
after  the  MTS'  RUN  switch  is  pressed  and  the  HOLD  button  is 
released. 

20.  When  the  desired  maximum  load  (approximately  500  lbs)  is 
reached;  depress  the  HOLD  button,  switch  to  load  control, 
and  manually  decrease  the  load  (to  approximately  20  lbs)  by 
SLOWLY  turning  the  SETPOINT  knob  clockwise.  As  with 
precrscking  you  must  end  the  unloading  with  a  small  increase 
in  load  (10  lbs  from  where  you  stop)  to  signal  the  computer 
you  have  completed  the  load/unload  cycle.   The  computer  will 
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plot  the  data  aa  you  go,  enaure  you  remain  on  the  linear 
portion  of  the  loading  curve. 

21.  The  computer-  will  display  all  the  pertinent  crack  length 
information.  You  may  repeat  this  operation  by  switching 
back  to  atroke  control  and  repeating  Steps  19  and  20  as 
often  aa  you  like  aa  long  aa  you  remain  on  the  linear- 
portion  of  the  loading  curve. 

22.  After  the  computer  haa  stored  the  desired  initial  crack 
length  (a0>  measurement  on  data  diskette,  you  are  ready  to 
begin  propagating  the  crack. 

23.  Enaure  you  are  in  atroke  control  with  near  zero  (<1C 
lbs)  load  and  the  RUN  light  is  on  with  the  HOLD  button 
depressed. 

24.  Press  the  k7-key  on  the  computer  follow  instructions  and 
proceed  to  collecting  data  (data  acquisition  system  will  be 
continuoaly  switching  channela  and  the  plotter  will  be 
plotting),  now  releaae  the  HOLD  button  on  the  MTS. 

25.  When  you  reach  a  point  where  you  wiah  to  unload  the 
specimen;  press  the  HOLD  button,  switch  to  load  control,  and 
slowly  decreaae  the  load.  The  computer  will  beep  when  you 
have  unloaded  the  apecimen  by  10%  of  the  expected  maximum 
load.  Switch  back  to  stroke  control  and  release  the  HOLD 
button  to  reload  the  apecimen. 

26.  When  the  load  returna  to  the  point  from  which  the 
unloading   began  the  computer  will  diaplay  and  printout   all 
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the  information  it  has  determined  with  regards  to  the 
current  unloading  cycle. 

27.  This  process  of  switching  control  modes  (Steps  25  and 
Step  26)  are  repeated  for  every  unloading  cycle  until  the 
you  decide  to  atop  by  pressing  the  k9-key  on  the  computer  or 
one  of  the  limit  detectors  shuts  off  the  MTS  machine. 
Remember,  you  must  have  atleast  four  data  points  between  the 
offsets  one  of  which  lies  near  the  blunting  line  on  the  J-R 
curve . 

28.  The  kS-key  interrupts  the  computer  program  and  tells  it 
to  stop  taking  data  and  write  the  data  it  has  collected  on 
the  data  diskette. 

29.  The  computer  will  then  display  a  table  of  the  collected 
J  versus  delta_a  and  ask  if  you  desire  a  plot.  Simply  follow 
the  directions  from  the  computer  and  you  will  obtain  the  J-R 
curve  and  the  value  of  Jic- 

D.   POST -TEST 

The  post- test  includes  verifying  the  final  crack  length, 
and  validating  the  other  parameters  as  discussed  in  ASTM 
Standard  E813-81  such  aai 

bi  and  Bmin  >  25  »  Jic  /  Fa 
where 

bi  =  uncracked  ligament  after  precracking 
Bmin  *  minimum  specimen  thickness 
Fs  =  flow  atreaa 
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and  that  the  alopa  of  the  linear  regresalon  fit  to  the  data 
ia  leaa  than  the  flow  atreaa  (dJ/da  <  Fa) .  Thi3  data  ia 
obtained  directly  from  the  " J_INTGL"  program  output . 

Heat  tinting  the  specimen,  10  minutea  at  300C  for  moat 
steels,  aide  in  the  viaual  identification  of  the  amount  of 
crack  extenaion.  After  the  apecimen  ia  removed  from  the 
oven  and  air  cooled,  place  it  in  liquid  nitrogen  to  assist 
in  obtaining  a  brittle  fracture.  After  breaking  the 
apecimen  apart,  place  the  piecea  in  alcohol  while  returning 
to  room  temperature.  Thia  will  prevent  water  vapor  from 
condensing  and  freezing  on  the  aurface  and  minimize 
oxidation. 

A  photographic  enlargement  of  the  aurface  can  now  be 
obtained  and  the  crack  length  can  eaaily  be  verified.  The 
final  crack  extenaion  muat  be  within  15%  of  that  recorded 
during  the  teat. 
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Figure  Al .   MTS  Control  Panel 
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APPENDIX  B 
TEST  RECORD  SHEET 

THREE  POINT  BEND  SPECIMEN 
TEST  REPORT 


SPECIMEN  I.O.. 
Data 


HATERIAL 


Operator. 


ORIENTATION. 
Taat  T«»d 


Baaaa 


Enia 


(Load  pt.) 
I* 


L  - 


Ha  ■ 


Baax  • 


Bain  ■ 


Utj 


Fa  ■ 


Polaaon'a  ratio  • 


Elaatic  Hodulua  • 


CoiMnu: 
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APPENDIX  C 
J  INTGL  FLOW  CHART 
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"7F" 


Initial^ 
a 


WHILE 

Flagl=0 

Flag2=0 


Data_acq 


TEST: 
If  Lodiff  -25.0 

Flagl=l 
If  Lodiff  25. 0 

Flag2-1 


Correct 
defl 


Slope 
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Main_prog 
Line  *  9380 


Line  #  9400 


Line  #  9480 


End  test 


I     3S         13      t  . 


Stor  arrays 


R  curve 


Line  #10430 
(calc  area) 

Line  #10440 
(reassign 
variables) 


Take  data 


-^j  Correct 
defl  " 


\/_ 


Slope 




Setup _plotter 
viz     g 


Ext ernal_pl otter 


V       YES 
(X 


L=L+1 


11  =  0 


No   crack  ext 
11=0 


5t* 


YES 


BEEP 
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APPENDIX  E 


PRECRACK  LISTING 
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280 

290 
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310 
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330 

340 
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360 

370 

380 
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400 

410 
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510 
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550 
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580 

590 


PRECRACKING  PROGRAM  FOR  THREE  °OINT 
BEND  SPECIMEN.  THIS  PROGRAM  INCORPORATES 
THE  USE  OF  THE  FOLLOWING  EQUIPMENT: 

1 .  HP-9826  COMPUTER 

2.  HP-3497A  DATA  ACQUISITION  SYSTEM 

3.  HP-2671G  THERMAL  PRINTER 

4.  NTS  MACHINE 

THE  PROGRAM  USES  BOTH  ENGLISH  AND  METRIC 
UNITS.  ALL  METRIC  VARIABLES  ARE  DENOTED 
UITH  AN  in  IN  THEIR  NAME  ie  Wm-WIDTH  IN 
METRIC  UNITS  (mm  IN  THIS  CASE). 

THIS  PROGRAM  HAS  WRITTEN  BY: 

LCDR  W.  K.  TRITCHLER,  USN 
VERSION  1.0     DATE  11/33 


COM  /J_intgl/  F? 

DIM  Lodm(2000) .Def lra<2000).Codm<2000)  these  arrays  will  be  dimensioned 

required  in  the  New  a  routine. 

Lodm<I>  -  LOAD  IN  NEW TONS 

Deflm(I)  -  DEFLECTION  IN  mm 

Codm<I)  -  CRACK  OPENING  DISPLACEMENT  IN  mm 

DIM  Coef<6>.G(6).A0<!00> 

DUMP  DEVICE  IS  706 

Dvm-709 

CLEAR  Dvm 

PRINT  CHRSU2)       !  CLEARS  CRT 

PRINT  — 

PRECRACK ING  PROGRAM" 
VERSION  1 .0" 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 


WRITTEN  BY  LCDR  W.K.  TRITCHLER* 
PLEASE  SELECT  A  k-KEY" 


Begin:  ! 
ON  KEY 
ON  KEY 
ON  KEY 
ON  KEY 


0  LABEL  "SPEC  DATA"  GOSUB  Spec_data 
2  LABEL  "CONV_DATA"  GOSUB  Conv_data 

4  LABEL  ~NEW_A"  GOSUB  New_a 

5  LABEL  Ma_TABLE"  GOSUB  A_table 


Spin: 


!  WAITS  FOR  KEY  INPUT 


GOTO  Spin 


Spec_data:    !  INPUTS  REQUIRED  SPECIMEN  DATA 
PRINT  CHR$<12)       !  CLEARS  CRT 
PRINT  "INPUT  SPECIMEN  IDENTIFICATION  <ID):' 
INPUT  DS 
PRINT  "" 
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PRINT   "J-INTEGRAL    TEST    Tor    ";DS 

PRINT   "" 

PRINT  "INPUT  SPECIMEN  DATA  USING  ENGLISH  UNITS  (LBS. 

PRINT  ••" 

PRINT  "WHAT  IS  THE  SPECIMEN'S  WIDTH.  NOTCH  LENGTH.  Bmax. 


600 
610 
620 
630 
640 
?•• 

650 

660 

670 

680 

630 

700 

IMATE  STRENGTH?" 

INPUT  E.Ys.Ut 
PRINT  CHRSU2) 
Fs-<Ys+Uts>/2 

PRINT  "USING  ENGLISH  UNITS  I in.lbs.psi>" 
PRINT  "THE  FOLLOWING  HAS  BEEN  ENTERED  for: 
PRINT  "WIDTH  -";W;"    NOTCH  LENGTH  -";A 
PRINT  "Bmax  -":Bmax:"        Bmm  «";Bmm 
PRINT  "HALF  SPAN  »":Hs;"    POISSON'S  RATIO 
PRINT  "YIELD  STRENGTH  -";Ys;" 
PRINT  "ELASTIC  MODULUS  -";E;" 


INPUT  W, A. Bmax. Bmm 
PRINT  "" 
PRINT  "WHAT  IS  THE  SPECIMEN'S  HALF  SPAN  AND  POISSON'S 

INPUT  Hs.Nu. 
PRINT  "" 
PRINT  "WHAT  IS  THE  SPECIMEN'S  ELASTIC  MODULUS.  YIELD 


INCHES,  etc.)" 
Bmm  (THICKNESS) 

RATIO?" 
STRENGTH  AND  ULT 


;Nu 

ULTIMATE  TS  -";Uts 
FLOW  STRESS  -";Fs 
PRINT  — 
PRINT  "THIS  PROGRAM  ALLOWS  NO  EASY  RECOVERY,  IF  ANY  DATA  IS  INCORRECT  YCU 


710 
720 
730 
740 
750 
760 
770 
780 
790 
800 
810 

MUST  RE 

830 

840 

850 

860 

870 

880 

890 

900 

910 

920 

930 

940 

PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 


!  CLEARS  CRT 


";DS 


!  REPEATS  Spec_data 


950 

960 

970 

980 

990 

1000 

1010 

1020 

1030 

T040 

1050 

1060 

1070 

1080 

1090 

1100 

1110 

1120 

1130 

1140 

1150 

1160 


PRINT  "EVERYTHING." 

PRINT  "" 

PRINT  "IS  ALL  DATA  CORRECT?" 

INPUT  HS 

IF  HS-"NO"  OR  HS-"N"  THEN  550 

PRINT  "" 

PRINT  "DO  YOU  WANT  A  HARD  COPY?" 

INPUT  HS 

IF  HS-"N"  OR  HS-"N"  THEN  1060 

j 

PRINTER  IS  706  !  TOGGLES  PRINTER  TO  THERMAL  PRINTER 

PRINT  "" 

USING  ENGLISH  UNITS  (in.lbs.psi)  THE  FOLLOWING  INFORMATIONS- 
HAS  BEEN  ENTERED  FOR:  ";DS 


WIDTH  -";W:" 
Bmax  -";Bmax:" 
HALF  SPAN  »":Hs: 
YIELD  STRENGTH  ■ 
ELASTIC  MODULUS 


PRINT 

PRINTER  IS  1 
t 

PRINT  CHRS<12> 
» 

!  CONVERSION  TO  METRIC  UNITS 

Conlen-25.4 

Wm-W»Conien 

Am-A»Conien 

Bmaxra"Bmax«Conlen 

Bmmm-Bmin»Conien 

Hsm-Hs»Con 1  en 

Conpress-6. 89475293 

Em-£«Conpress 


NOTCH  LENGTH  -";A 
Bmm  -";Bmm 
"  POISSON'S  RATIO  >";Nu 
";Ys;"  ULTIMATE  TS  -";Uts 
-";E;"  FLOW  STRESS  -";Fs 


♦  TOGGLES  PRINTER  BACK  TO  CRT 
!  CLEARS  CRT 


!  1  in  -  25.40mm 


!  tpst  -  6.89475293kPa  <HP  PROGRAMING  BOOK) 
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It  70 
1180 
H90 
1200 
1210 
DS:" 
1220 
1230 
1240 
1250 
12S0 
1270 
1230 
1290 
1300 
1310 
1320 
1330 
1240 
1350 
1360 
1370 
t380 
1390 
1400 
1410 
1420 
1430 
1440 
1450 
1460 
1470 
1480 
1490 
1500 

1510  Neu_a: 
i 


Ysm- Ys*Conpress 

Utsm-Uts**>onpress 

Fsm»Fs»Conpress 
i 

PRINT  "CONVERTING  THIS  DATA  TO  METRIC  UNITS  (mm.N.KPa) 


for  ";DS;"  YIEL 


PRINT  "WIDTH 
PRINT  "Bmax  ■ 


":14m:"    NOTCH  LENGTH  -";Am 

:6maxm:"       Btnin  »":Bnuntn 
PRINT  "HALF  SPAN  >";Hsm:"    POISSCN'S  RATIO  -"jNu 
PRINT  "YIELD  STRENGTH  -";Ysm:"   UTILMATE  TS  -":Utsn 
PRINT  "ELASTIC  MODULUS  -";Em;"   FLOW  STRESS  -";Fsm 
PRINT  "" 

PRINT  "DO  YOU  WISH  TO  HAVE  A  PRINT  OUT  OF  THIS  DATA?" 
INPUT  G3 

IF  GS»MN"  CR  GS-"NO"  THEN  1410 
t 

PRINTER  IS  706  !  TOGGLES  TO  THERMAL  PRINTER 

PRINT  "CONVERTING  THIS  DATA  TO  METRIC  UNITS  (mm.N.KPa)  FOR: 

PRINT  "WIDTH  -":Wm:"      NOTCH  LENGTH  -";Am 

PRINT  "Bmax  "";Bmaxm:"    Bmin  "";Bmmm 

PRINT  "YIELD  STRENGTH  -";Ysm;"   ULTIMATE  TS  «":Utsm 

PRINT  "ELASTIC  MODULUS  »";Em;"   FLOW  STRESS  -";Fsm 

PRINT  "" 

PRINTER  IS  1 
r 

PRINT  CHRSU2) 
t 

Flag3-0 

Al  =A 

GOSUB  Params 

RETURN 


":D£ 


!  TOGGLES  PRINTER  TO  CRT 
!  CLEAR  CRT 


!  THIS  ROUTINE  COMPUTES  NEW  CRACK  LENGTH  AND  NEW  PARAMETERS 


1520 
1530 
154& 
1550 
1560 
1570 
1580 
1590 
1600 
16t0 
1620 
1630 
1640 
1650 
1660 
1670 
1630 


Beffra-Bma*m-(Bmaxni-Bminm>~2/Bma*m     !  B  EFFECTIVE  FROM  ASTM  E813-81 

D2-1  !  INITIALIZES  COUNTER 

PRINT  CHRSU2)    !  CLEARS  CRT 

PRINT  "THIS  ROUTINE  WILL  OBTAIN  rHJRRENT  CRACK  LENGTH" 

PRINT  "  AND  NEW  PARAMETERS  FOR  PRECRACKING." 

PRINT  "BEFORE  YOU  CAN  RUN  THIS  ROUTINE  YOU  MUST:" 


ENTER  SPECIMEN  DATA" 
ENTER  CONVERSION  DATA" 
SETUP  MTS  FOR  MANUAL  LOADING  OF  THE"; 
SPECIMEN." 


ARE  YOU  READY?* 


PRINT  "    1. 

PRINT  "    2. 

PRINT  "    3. 

PRINT  " 

PRINT 

PRINT 

INPUT  HS 

IF  H$-"N"  OR  HS-"NOM  THEN  2450   !  RETURNS  TO  Ic-KEY  SELECTION 

PRINT  CHRSU2)       !  CLEARS  CRT 

PRINT  "YOU  MAY  LOAD  AND  UNLOAD  SPECIMEN  AS  OFTEN  AS   YOU  LIKE  —  BUT  YOU  M 
UST  REMAIN  "; 

1690   PRINT  "BELOW  CURRENT  PMAX  OF" :Pmax ;"Lbs" 
1700   DIM  Lod<2000).Def M2000) ,Cod<2000> 
1710   DIM  Lodm(2000).Deflm<2000),Codm(2000) 

1720  PRINT  "WHEN  YOU  ARE  READY  TO  BEGIN  PRESS  THE  'ENTER'  OR  'CONTINUE'  KEY." 
1730   INPUT  HS 
1740   PRINT  CHRS<12)  !  CLEARS  CRT 
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USED  AS  COUNTER 
USED  AS  INDEX 


1750  K-0 

T760   I-l 

1770  Flaql-0 

1780  Flag2-0 

1790     WHILE  FUgl-0 

1800        GOSUB  Data_acq 

1810     END  WHILE 

1820     Flag2-0         ! 

1830     WHILE  Flag2=0 

1840        GOSUB  Data_acq 

1850     END  WHILE 

1360   ! 

1370   !  THIS  SECTION  CAN  BE  USED  TO  VIEW  THE  COLLECTED  DATA  3Y  REMOVING  THE  REMA 


!  ESTABLISHES  ONE  LOAD-UNLOAD  DATA  TAKING  LOOP 


RESETS  FLAG2 


RK  CHARACTER  < ! > 


1880 
1890 
1900 
1910 
1920 
1930 
1940 
1950 
1960 
1970 
1980 
1990 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2080 
2090 
2100 
2110 
2120 
2130 
2140 
2150 
2160 
2170 
2180 
2T90 
2200 
2210 
2220 
2230 
2240 
2250 
2260 
2270 
2280 
2290 
2300 
2310 
2320 


FOR  1-1 
!  PRINT 
!  IMAGE 
NEXT  I 
WAIT  5 
■1 


TO  K 

USING  1960;I.LodnHl)  ,Deflm<I>  ,Codm<I) 

X.3D,4X,6D.D,4X.3D.7D.7X,3D.7D 


Slope 
Crack_length 


!  SLOPE 
!  USES 


K1 

K3»< 

GOSUB 

GOSUB 

» 

A0(D2)-CraclcIm/25.4 

A1-A0<D2) 
i 

!  FIND  95%  CONFIDENCE 
Stiffm-Stiffnn-T9-U4 
GOSUB  Crack_Iength 
Amax-Cracklm/25.4 
Stiffra-Stiffa»-2»T9«U4 
GOSUB  Crack .length 
Anun-Crack  lm/25.4 


•  I /COD  COMPLIANCE 
COD  COMPLIANCE 


!  CONVERTS  TO  INCHES 


INTERVAL  FOR  CRACK  LENGTH 


!  CONVERTS  TO  INCHES 


CONVERTS  TO  INCHES 


j 

PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
INPUT 
IF  G$- 

PRINTER 

PRINT  " 

PRINT  " 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT  " 

PRINTER 
» 


CHRS<12>     ! 
"THE  FOLLOWING 


CLEARS  CRT 
INFORMATION  HAS 


BEEN  OBTAINED:" 


CRACK  LENGTH  -":A0<D2> :"in     FROM  MEASUREMENT  *";D2 
WITH  952  CONFIDENCE  INTERVAL  FROM" ;Anun;"  TO";Amax 


USING  COD  COMPLIANCE  0F:";Compl 
WITH  A  CORRELATION  OF:";Correl ;' 


"in/Lb    '";Complm',"mm/H' 
FROM" ;K; "DATA  POINTS" 


"DO 
G$ 
N" 


YOU  WANT  A  HARD  COPY?" 
OR  GS-"NQ"  THEN  2340 


IS  706 

M 

THE  FOLLOWING 
CRACK  LENGTH  - 


•  TOGGLES  TO  THERMAL  PRINTER 

INFORMATION  HAS  BEEN  OBTAINED: 
";A0<D2>;"in  FROM  MEASUREMENT 


*";D2 
"WITH  952C0NFIDENCE  INTERVAL  FROH";Amin: "  TO" ; Amax : " i n" 
"USING  COD  COMPLIANCE  OF:" jCompl ;" in/Lb  -" ;Complm:"mm/N" 
"WITH  A  CORRELATION  OF:" :Correl ;"FROM":K;"  DATA  POINTS." 


IS  1 


!  TOGGLES  PRINTER  TO  CRT 
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DO  YOU  WANT  ANOTHER  MEASUREMENT?' 


2330 

2340 

2350 

23G0 

2370 

2380 

2390 

2400 

2410   IF  A0(D2-!)>.4-0eIa+A  THEN  Fiag3-1 

F  PRECRACKING 


PRINT  ' 

PRINT  ' 

INPUT  FS 

D2-D2+1       !  INCREMENTS  COUNTER 

IF  F$-"Y"  OR  F$-"YES"  THEN  1670 

» 

IF  A0<D2-1)>0pta  THEN  GOSUB  End_precrack   !  PRECRACKING  IS  COMPLETE 
t 


!  SETS  FLAG  INDICATING  IN  LAST  60/.  0 


GOSUB  Params 
i 

RETURN 


2420 
2430 
2440 
2450 
2460 
2470 
2430 
2490  Data_acg:    !  DATA  ACQUISITION  ROUTINE 

OUTPUT  Dvm;"VR5  AF1  AL3" 

; 

OUTPUT  Dv»;°flll  VT1" 
ENTER  Dvm;Lod<I> 
Lodm<I)-Lod(I)«Slm(t) 
f 

OUTPUT  Dvm:"AI2  VTV 

ENTER    Dvm:Cod(I> 

Codm<I)-Cod<I)«Slm<2> 
; 

OUTPUT  Dvn:"AI3  VT1" 
ENTER  Dvm:DefKI> 
DefIm<I)-Defl(I)»Slm<3) 


2500 
2510 
2520 
2530 
2540 
2550 
2560 
2570 
2580 
2590 
2600 
2610 
2620 
2630 
2640 
2650 
2660 
2670 
2680 
2690 
2700 
2710 
2720 
2730 
2740 
2750 
2760 
2770 
2780 
2790 


!  READS  LOAD 


!  READS  COD 


!  READS  DEFLECTION 


TEST  DATA  TO  BREAK  OUT  OF  DATA  ACQ  LOOP 


IF  1-1  THEN  2720  ■    !  SKIPS  IF  FIRST  DATA  POINT 
Lodxffm-Lodm(I)-Lodm<I-1) 

IF  Lodiffm<-25.0  THEN  Flagl-1 

IF  Lodiffn>25.0  THEN  Flag2-1 
Lodiffm-0 

!  INCREMENTS  INDEX 


!  MUST  BE  UNLOADING 
!  MUST  BE  RELOADING, 


1-1*1 
i 

k»K+i 

» 

RETURN 
| 


INCREMENT  COUNTER 


2800  Craclc_length:    !  CALCULATES  CRACK  LENGTH  USING  HUDAK  (DTNSRDC  PROG) 

2810  !   NOTE:  THIS  ROUTINE  USES  COD  COMPLIANCE  METHOD  VICE  DEFLECTION  COMPLIA 

NCE. 

2820  FOR  I-t  TO  6 

2830  READ  Coef(I) 

2840  NEXT  I 

2850  Complm-1/Stiffm 

2860  C-2»Um-Beffm-£ni«Complm/<Hsin-1000> 

2870  C-1/<SQR<C>+1) 

2880  G<T)-C 

2890  FOR  1-2  TO  5 

2900  G<I)-G<I-1)-G<1) 
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2910 
2920 
2930 
2940 
2950 
2960 
2970 
2980 
2990 
3000 
3010 
3020 
3030 
3040 
3050 
30S0 
3070 
3080 
3090 
3100 
3110 
3120 
3130 
3140 
3150 
3160 
3170 
3180 
3190 
3200 
3210 
3220 
3230 
3240 
3250 
3260 
3270 
3230 
3290 
3300 
3310 
3320 
3330 
3340 
3350 
3360 
3370 

N 

3330 
3390 
3400 
3410 
3420 
3430 
3440 
3450 
3*60 
3470 
3480 
3490 


next  r 

F'Coef ( t  ) 
FOR  1-1  TO  5 

F-F+Coef(IH)-G(I) 
NEXT  I 

Crack  Ini-F-Wm 
RESTORE 

DATA  0.998265,-3.81662.-1 .80596,32.31038, -44 . 15665 . -52 .67376 
i 

RETURN 


Conv_data: 

PRINT  CHR5(12) 

PRINT  "INPUT:  LOAD  RANGE 

PRINT  •"• 

PRINT  "   (BY  PERCENT 

INPUT  Lr.Cr.Sr 
PRINT  •- 

PRINT  "LOAD  RANGE  -  ":lr 
PRINT  "STRAIN  <COD>  RANGE 
PRINT  "STROKE 
Lmax-20000 
Cmajc«.20O 
Smax-5 .0 

V— 10 
i 


INPUTS  CONVERSION  DATA  FOR  MEASURING  DEVICES 
!  CLEARS  CRT 
STRAIN  RANGE. 


STROKE  RANGE" 


le.  10.20.50  OR  100)" 


,#:Cr 
(DEFLECTION)  RANGE  -  ":Sr 
!  MAX  LOAD  RANGE  IN  (lbs). 
!  MAX  COD  RANGE  IN  (in). 
!  MAX  STROKE  RANGE  IN  (in) 
I    MAX  VOLTAGE  IN  (volts). 


!  MAXIMUM  OPERATING  POINTS  ARE: 

Olmax-Lmajt-Lr/100 

Ocmax»Cmajt-Cr/100»<-1.0)   !  -1  CHANGES  SIGN 

Osmajt-Soax-Sr/100 

j 

!  CONVERSION  CONSTANTS  ARE:  (ENGLISH  UNITS) 

Sl(1)-Olmax/V 

Sl(2)-0cmax/V 

Sl(3)-0smax/V 

♦  CONVERTING  TO  METRIC: 

Conlen-25.4 

Conib-4. 448221615 

Slmd  )-Sl(1Xonlb 

Slm(2)-Si(2)-ConIen 

Sltn(3)-SM3>»Conien 

i 

PRINT  ••" 
PRINT  "- 
PRINT  "YOUR  ENTRIES  HAVE  PRODUCED  THE  FOLLOWING 


!  25 . 4mm  -  1  in 
!  4.448221615N  - 


lib 


(HP  PROGRAMING  BOOK) 


CONVERSION  DATA: 


;Slm(l);"NEWTON/VOLT   <";S1( 1 ) ;"Lbs/V)" 
<";SK2):"in/V)" 


PRINT  ••" 

PRINT  "FOR  LOAD 

PRINT  — 

PRINT  "FOR  COD:  ";Slm<2> :"mm/V 

PRINT  "" 

PRINT  "FOR  DEFLECTION:  ";Slm(3)  ;"mm/V  ("  ;SK3)  ;"in/V)" 

PRINT  "" 

PRINT  "WOULD  YOU  LIKE  A  HARD  COPY?" 

INPUT  HS 

IF  HS-"N"  OR  HS-"NO"  THEN  3570 
f 

PRINTER  IS  706  !  TOGGLES  PRINTER 


61 


3500  PRINT  -• 

3510  PRINT  -       CONVERSION  DATA" 

3520  PRINT  "FOR  L0AD:";Slm<1 ):"N/V   <";S1<1 > ;"Lbs/V>" 

3530  PRINT  "FOR  COD: " :Slm<2>  :"mm/V   (";SK2> :  "in/V>" 

3540  PRINT  "FOR  DEFLECTION: " ;Slm<3> :"mm/V   <"  :S1<3>  ;"  m/V>" 

3550  PRINTER  IS  1            !  TOGGLES  PRINTER 

3560  ! 

3570  PRINT  CHRSM2) 

3580  ! 

3590  RETURN 

3600 

3610 

3620 

3630  Slope:    !  LINEAR  LEAST  SQUARES  FIT 

3640  PRINT  CHRS(12)        !  CLEAR  CRT 

3650  !  INITIALIZE  VARIABLES 

3660  Sumx-0     !  SUM  OF  X  VALUES 

3670  Sumx2-Q    !  SUM  OF  X  SQUARED  VALUES 

3680  Swny-Q     !  SUM  OF  Y  VALUES 

3690  Sumy2"0    !  SUM  OF  Y  SQUARED  VALUES 

3700  Sujnxy-0    !  SUM  OF  THE  PRODUCT 

3710  FOR  I-K1  TO  K3    !  K1 -INDEX  OF  FIRST.  K3-INDEX  OF  LAST  DATA  POINT  USED, 

3720  ! 

3730  X-Codm<I> 

3740  Y-Lodm(I) 

3750  ! 

3760  Siunx-Sumx+X 

3770  Sumy-Smny+Y 

3780  Sumx2-Sujnx2+X»X 

3790  Suii»y2-Sujny2*-Y»Y 

3800  Smnxy-Sumxy+X»Y 

3810  NEXT   I 

3820  K3-K3-P2               !   USED   TO  CALCULATE  NUMBER  OF   DATA   POINTS  WITHIN   BOUNDS. 

3830  K4-K3-K1+1                    !    K4-T0TAL   #   DATA   POINTS  USED    IN   CALC 

3840  QaI-Sumx2-Suox«Sunx/K:4 

3850  Qa2-Sumxy-Sujnx*Sujny/K4 

3860  Qa3-Sumy2-Suj»y-Suiny/K4                  •<.- 

3870  ! 

3880  M-Qa2/Qa1            !  SLOPE 

3890  ! 

3900  IF  M<0  THEN  4190   !  IF  SLOPE  <  0,  BAD  DATA 

3910  !  FROM  DTNSRDC  PROGRAM 

3920  Q1«<K4»Sumx2-Su«x~2>/<K4-<K4-1>> 

3930  Q2-<K4»Sumy2-Suiny*2)/(K4»<K4-1  )> 

3940  Q3-<K4-1 )»<Q2-M»M-Q1 >/<K4-2> 

3950  Q4-SQR<Q3/(Q1»<K4-1))) 

3960  ! 

3970  B-<Suay-M-Sunx>/K4     !  INTERCEPT 

3980  ! 

3990  Correl-M»SQR<Qt/Q2) 

4000  ! 

4010  Qa4-SQR<ABS<<Sufliy2-B»Sui»y-M»Sufl»xy>/<K4-2))) 

4020  Signab-Qa4/SQR<ABS<Qa1>> 

4030  Sigmaffl-Signab»SQR(Sunx2/K4) 

4040  Stiffm-M                                                !    STIFF   -   SLOPE 

4050  Conplm-I/M                  !    COD  COMPLIANCE    -    1 /SLOPE 

4060  Con»pl-Complm«<4.448221615>/25.4      !    CONVERTS    TO      in/LB 

4070  ! 

4080  !    STUDENT    t   TEST 

4090  U4-Q4 
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4'00  T9-I.96 

41  10  IF  K4>120  THEN  4170 

4120  T9-2.1 

4130  IF  K4<20  THEN  4170 

4140  N9-<4-2 

4150  T9-N9/<-.617109+N9». 510239) 

4160  ! 

4170  RETURN 

4130  ! 

4190  PRINT  "POOR  DATA--  SLOPE  <  0" 

4200  M-333333.3 

4210  Correi«0 

4220  ! 

RETURN 


RETURNS 
RETURNS 


4230 

4240 

4250 

42S0 

4270  Params:     !  THIS  ROUTINE  UILL  DETERMINE  THE  CURRENT  MAX  AND  MIN  LOAD  FOR 

FATIGUE  PRECRACKING 


4280 

4290 

4300 

4310 

4320 

4330 

4340 

4350 

4360 

4370 

4380 

4390 

4400 

4410 

4420 

4430 

4440 

4450 

ACKING 

4460 

4470 

4430 

4490 

4500 

4510 

4520 

4530 

4540 

4550 

4560 

4570 

4580 

4590 

4600 

4610 

4620 

4630 

4640 


i 

!  NOTE  THE  USE  OF  ENGLISH  UNITS  DUE  NTS  READOUT  BEING  IN  POUNDS, 

!  INITIALIZE  VARIABLES 

Fau-0 

Pmax-0 

Pmin-0 

Deip-0 

Kmaxp»Q 

Knunp»Q 

Dellf-0 

Dellce-0 
; 

6pta-.6*U 

Dela-Opta-Al 

Dela50-Dela»< 
f 

Pmax-<4/3)»<Btnin-<W-Al)"2»Fs/<2»Hs>>    !  FROM  ASTM  E8T3 

IF  Flag3-1  THEN  Pmax-< .4)^Pmax        !  MAX  LOAD  DURING  LAST  50%  OF  PRECR 

(ASTM  E813) 

Pmin-Pmajt-.95»Pmaj»      !  FROM  ASTM  E813 

Delp-Ptnax-Pmin 

DETERMINE  Del*   BASED   ON  P  USING  ASTM  E399 

5«<1.99-<A1/W)*<1-AI/W)«<2.15-3.93»AI/W+2.7«<A1/U)'2>) 


!  OPTIMUM  CRACK  LENGTH 
!  EXTENSION  REQUIRED 
,5>+A   !  CRACK  LENGTH  AT  50%  OF  EXTENSION 


Faun-3-<Al/H) 
; 

Faud-2»<1+2»A1/U)»<!-Al/W)~1 .5 
; 

Fau"Faun/Faud 
f 

i 

Kmaxp-<Proax«2»Hs/(Bnun*U~T  .5))»Faw  ! 

Kmino-<Pmin«2»Hs/(Bmin»t>»"1  .5)  >»Fau  ! 

Delk-Kmaxp-Knunp 
f 

I  DETERMINE  Delke  BASED  ON  MODULUS  FROM  ASTM  E813 

IF  Flag3-0  THEN  Delke-40000         !  DURING  EARLY  EXTENSION 

IF  Flag3-1  THEN  Delke-< .001 )•£      ♦  DURING  LAST  507.  OF  PRECRACKING 

F.  ASTM  E399  4  E8131 

4650   PRINT  "" 


FROM  ASTM  E399 
FROM  ASTM  E399 


[RE 
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4660  PRINT  "BASED  ON  CURRENT  INFORMATION:" 

4670  PRINT  "" 

4680  PRINT  "CURRENT  CRACK  LENGTH  »";Al;"in" 

4690  PRINT  "OPTIMUM  aO-" ;Gpca:"in" 

4700  PRINT  "EXTENSION  REQUIRED  -":Deia 

4710  PRINT  "CRACK  LENGTH  AT  507.  OF  EXTENSION  -"  ;Dela50  ;"in" 

4720  PRINT  "Pmax  -":Pmax:"Lbs  Pmn  -"  ;Pnun:"Lbs" 

4730  PRINT  "KmaxD  -";Kmaxp:"   Kminp  -";Kminp 

4740  PRINT  "DELTA_K  -"  ;Delk  ;"m  '  .5" 

4750  PRINT  "DELTA_KE  -" ;0elke ;"in* .5" 

4760  PRINT  "FAW  -";Fau      !  SERVES  AS  A  CHECK 

4770  PRINT  MM 

4780  PRINT    "IF   Delk    >    DeJIce.    Pmax    MUST    BE    DECREASED" 

4790  IF   Deik<Delke   THEN   4900  !    SKIP    IF   OKAY 

4800  !    COMPUTE   NEW  Pmax 

4810  Deip1-<Delke/Fau/<2»Hs)«Bnun-W~1  .5) 

4820  Pmax1-Delpl/.95 

4830  Pmin1-Pmax1-.95*Pmax1 

4840  Kmax  !  -  <  Pmax  1  «2*Hs/  <  Bm  i  n»v4*  1.5))  »Fau* 

4850  PRINT    "THE   NEW   Pmax    BASED   ON    Delke    IS :" ;Pmax 1 ;"Lbs" 

4860  PRINT    "NEW   Pmn   -":Pnun !  :"Lbs" 

4870  PRINT   "NEW   DELTA_P   -" ;Delp1 ;"Lbs" 

4880  PRINT   "NEW  Kmax    «";Kraaxt 

4890  ! 

4900  PRINT  "   DO  YOU  WANT  A  HARD  COPY?" 

4910  INPUT  HS 

4920  IF  HS-"N"  OR  HS-"NO"  THEN  5080 

4930  PRINTER  IS  706       *  TOGGLES  PRINTER 

4940  PRINT  "BASED  ON  CURRENT  INFORMATION  FOR  ";DS;":" 

4950  PRINT  "" 

4960  PRINT  "CURRENT  CRACK  LENGTH  -":Al;"in" 

4970  PRINT  "OPTIMUN  CRACK  LENGTH  «";Opta:"in" 

4980  PRINT  "EXTENSION  REQUIRED  -";Dela:"in" 

4990  PRINT  "CRACK  LENGTH  AT  502  OF  EXTENSION  •" ;Dela50 ;"in" 

5000  PRINT  "Pmax  -";Pmax;"Lbs  Pmin  -" ;Pmin;"Lbs" 

5010  PRINT  "Delp  -";Delp;"Lbs  Kmax  -":Kmaxp 

5020  PRINT  "DELTA_K  -";Oelk :"in* .5  DELTA_KE  -";Delke;"in- .5" 

5030  IF  Delk<Deike  THEN  5090 

5040  PRINT  "" 

5050  PRINT  "NEW  LOADS  BASED  ON  DELTA_KE  ARE:" 

5060  PRINT  "Pmax  »";Pmax1 :"Lbs  Pmin  »" ;Pmin1 ;"Lbs" 

5070  PRINT  "DELTA_P  »";Delpl ;"Lbs   Kraax  -";Kmax1 

5080  PRINT  "" 

5090  PRINTER   IS   1  !    TOGGLES  PRINTER    TO  CRT 

5100  PRINT   "" 

5110  PRINT  "MAKE  THE  NECESSARY  ADJUSTMENTS  BEFORE  CONTINUING." 

5120  PRINT  "WHEN  COMPLETED  PRESS  'CONTINUE'." 

5130  INPUT  CS 

5140  ! 

5150  RETURN 

5160  ! 

5170  » •«••••••••••*•••••••••••••»«••••«•••,•••«• 

5180  ! 

5190  A_table:       !  PRODUCES  TABULATED  VALUES  FOR  LOAD  vs  COD  VALUES  AT  VARIOU 

S  CRACK  LENGTHS 

5200  ! 

5210  PRINT  "BEFORE  YOU  CAN  RUN  THIS  SUBROUTINE  YOU  MUST  ENTER:" 

5220  PRINT  "      1 .   SPECIMEN  DATA" 

5230  PRINT  "      2.   CONVERSION  DATA" 
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5240  PRINT  "" 

5250  PRINT  "ARE  YOU  READY?" 

5260  INPUT  HS 

5270  IF  HS-"N"  OR  H$-"NO"  THEN  5550    !  RETURNS  TO  If -KEY  SELCETION 

5280  ! 

5290  Beff-Bmax-<Bmax-Bmin)~2/Bmajt             !    IAW   ASTM   E313-31 

5300  PRINT   "BEFF   -";Beff 

5310  ! 

5320  !  USING  CRACK  LENGTH  EQUATION  FROM  DTNSRDC  PROGRAM  CALC  CRACK.  LENGTH  FOR 

VARIOUS  VALUES  OF  U 

5330  A00-. 998265 

5340  A1— 3.81662 

5350  A2--1. 80596 

5360  A3-32.31038 

5370  A4-- 44.15665 

5380  A5--52. 67876 

5390  FOR  U-.07  TO  .155  STEP  .005 

5400  Ar-A00*A1»U+A2»U~2+A3»U*3+A4«U~4+A5»U*5 

5410  CIrl'Ar»W 

5420  PRINTER  IS  706    !  TOGGLES  TO  THERMAL  PRINTER 

5430  PRINT  "U  -":U 

5440  PRINT  "a/U  -":Ar:"in  CRACK  LENGTH  -":Ckl  :"liT 

5450  PRINT  "LOAD  (lbs)      COD  <in>       COD  (Volts)" 

5460  FOR  Lodl-200  TO  1800  STEP  100 

5470  C-Beff-W»£*2/Hs 

5480  Cod1»Lod1-<(1-U)/U)~2/C 

5490  Codv-Cod1/SK2) 

5500  PRINT  Lodt.Codl.Codv 

5510  NEXT  Lodl 

5520  NEXT  U 

5530  PRINTER  IS  1     »  RETURNS  TO  USING  CRT  AS  PRINTER 

5540  ! 

5550  RETURN 

5560 

5570 

5580 

5590  End_precrac»- :        !  END  OF  PRECRACKING 

5600  ! 

5610  PRINT  "     PRECRACKING  IS  COMPLETED!" 

5620  PRINT  "" 

5630  PRINT  "FINAL  CRACK  LENGTH  -";A0<D2-1> 

5640  PRINT  "CORRELATION  OF  DATA  -";Correl 

5650  PRINT  — 

5660  PRINT  "    REMOVE  SPECIMEN  FROM  MTS." 

5670  ! 

5680  END 
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APPENDIX  E 


J  INTGL  LISTING 


to 

20 

30 

40 

50 

60 

70 

80 

90 

100 

110 

120 

130 

140 

150 

160 

170 

180 

130 

200 

2T0 

220 

230 

240 

250 

260 

270 

280 

290 

300 

310 

320 

330 

340 

350 

360 

370 

380 

390 

400 

410 

420 

430 

440 

450 

460 

470 

480 

490 

500 

510 

520 

530 

540 

550 

560 

570 

580 

590 


J-INTEGPAL  FRACTURE  TOUGHNESS  PROGRAM 
FOR  THREE  POINT  BEND  SPECIMEN.   THIS 
PROGRAM  INCORPORATES  THE  USE  OF  THE 
FOLLOWING  EQUIPTMENT: 
I •  HP-9826  COMPUTER 

2.  HP-3497A  DATA  ACQUISITION  SYSTEM 

3.  HP-3437A  SYSTEM  VOLTMETER 

4.  HP-2671G  THERMAL  PRINTER 

5.  HP-7225B  PLOTTER 

THIS  PROGRAM  WAS  WRITTEN  BY: 

LCDR  W.  K.  TRITCHLER.  USN 
VERSION  t.O     DATE  11/33 


COM  /J_mtgi/  F1 

DIM  Lod'8000) .Def 1 (8000) .CodOOOO) 

NOTE:  Lod.  Defl  AND  Cod  ARE  ACTUALLY  DIMENSIONED  IN  Initiai_a  ROUTINE 

Lod<I)  -  LOAD 

Def  I  (I)  -  DEFLECTION 

Coa<I>  -  CRACK  OPENING  DISPLACEMENT  (COD) 

DIM  A0C50) .Amin0(50> . Amax0(50> .Comol0<50> ,CorrelO(50> 
DIM  Area9( 100) .Aid  00) .Coef (6) .G<6> 
DIM  J9< 100 >.Delta_a(1 00 > .rtdef 1(100) 
DIM  DataSd  :4)C10] 

IS  706    !  THERMAL  PRINTER  ADDRESS 
!  DATA  ACQUISITION  ADDRESS 


DUMP  DEVICE 
Dvm-709 
CLEAR  Dvm 
PRINT  CHRSU2) 
PRINT  "•• 
PRINT  " 
PRINT  M 
PRINT  " 
PRINT  "" 
PRINT  " 


!  CLEARS  CRT 

-J- INTEGRAL  TEST  PROGRAM- 
VERSION  1  .0" 
WRITTEN  BY  LCDR  W.K.  TRITCHLER" 

PLEASE  SELECT  A  Ic-KEY" 


Begin:  ! 
ON  KEY 
ON  KEY 
ON  KEY 
ON  KEY 
ON  KEY 
ON  KEY 
ON  KEY 
!<»)  ON 
!  (SEE 
ON  KEY 


0  LABEL  "SPEC_DATA"  GOSUB  Spec_data 

1  LABEL  "CONV  DATA"  GOSUB  Conv  data 


2  LABEL  "DAT_FILE"  GOSUB  Data_file 

3  LABEL  "PLOTTER"  GOSUB  Setup_plotter 

4  LABEL  "GCLEAR"  GOSUB  Gclr 

5  LABEL  "aO"  GOSUB  Initial.a 
7  LABEL  "START"  GOSUB  Mam_prog 

KEY  8  LABEL  "R-CURVE"  GOSUB  R_curve 
NOTES  ON  <•)  IN  R_curve  ROUTINE 
9,4  CALL  Stop_test 


Spin:     !  WAITS  FOR  KEY  INPUT 
GOTO  Spin 
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600  Gel r:      !  GRAPHICS  SCREEN  CLEAR 

610  !  THIS  ROUTINE  ALLOWS  SINGLE  KEY  STROKE  TO  REMOVE  A  GRAPHICS  DISPLAY  FROM 

CRT 

620  GCLEAR 

630  RETURN   !  RETURNS  TO  <-KEY  SELECTION 

640  • 

650  * •••••••••»••*••••••*••*•«•**•••••«•»•• 

650  ! 

670  Spec.data:    !  INPUTS  REQUIRED  SPECIMEN  DATA 

660  PRINT  CHRS(12)       !  CLEARS  CRT 

630  PRINT  "INPUT  SPECIMEN  IDENTIFICATION  (ID):" 

700  INPUT  DS 

7!0  PRINT  "" 

720  PRINT  "J-INTEGRAL  TEST  for  ";DS 

730  PRINT  "" 

740  PRINT  "INPUT  SPECIMEN  DATA  USING  ENGLISH  UNITS  (LBS.      INCHES,  etc.)" 

750  PRINT  "" 

760  PRINT  "WHAT  IS  THE  SPECIMEN'S  WIDTH.  NOTCH  LENGTH,  Bmax .   Bmin  (THICKNESS) 

?" 

770  INPUT  W.A.Bmax.Bnun 

780  PRINT  "" 

790  PRINT  "WHAT  IS  THE  SPECIMEN'S  HALF  SPAN  AND  POISSON'S     RATIO?" 

800  INPUT  Hs.Nu. 

8t0  PRINT  "" 

820  PRINT  "WHAT  IS  THE  SPECIMEN'S  ELASTIC  MODULUS.  YIELD      STRENGTH  AND  ULT 

IMfiTE  STRENGTH?" 

820  INPUT  E.Ys.Uts 

840  PRINT  CHRS(12)  J  CLEARS  CRT 

850  Fs-(Ys+Uts)/2 

860  PRINT  "USING  ENGLISH  UNITS  (in. Ibs.psi)" 

870  PRINT  "THE  FOLLOWING  HAS  BEEN  ENTERED  for:  ";D£ 

880  PRINT  "WIDTH  «";W:"    NOTCH  LENGTH  »":A 

890  PRINT  "Bmax  »";Bmax:"        Bmin  -";Bmin 

900  PRINT  "HALF  SPAN  -":Hs:"   POISSON'S  RATIO  »";Nu 

910  PRINT  "YIELD  STRENGTH  -";Ys:"  ULTIMATE  TS  -":Uts 

920  PRINT  "ELASTIC  MODULUS  -";E;"   FLOW  STRESS  -";Fs 

930  PRINT  -' 

940  PRINT  "THIS  PROGRAM  ALLOWS  NO  EASY  RECOVERY,  IF  ANY  DATA  IS  INCORRECT  YOU 

MUST  RE-ENTER 

950  PRINT  "EVERYTHING." 

960  PRINT  "" 

970  PRINT  "IS  ALL  DATA  CORRECT?" 

980  INPUT  H$ 

990  IF  HS-"NO"  OR  HS-"N"  THEN  670   !  REPEATS  Spec_data 

1000  PRINT  "" 

1010  PRINT  "DO  YOU  WANT  A  HARD  COPY?" 

1020  INPUT  HS 

1030  IF  HS-"N"  OR  H$-"N"  THEN  1180 

1040  ! 

1050  PRINTER  IS  706         !  TOGGLES  PRINTER  TO  THERMAL  PRINTER 

1060  PRINT  "" 

1070  PRINT  "USING  ENGLISH  UNITS  <in.Ibs.psi>  THE  FOLLOWING  INFORMATION" 

1080  PRINT  "HAS  BEEN  ENTERED  FOR:  ";DS 

1090  PRINT  "" 

1100  PRINT  "WIDTH  -":W:"       NOTCH  LENGTH  -";A 

1110  PRINT  "Bmax  -":Bmax;"      Bmin  -";Bmin 

1120  PRINT  "HALF  SPAN  »":Hs;"  POISSON'S  RATIO  »";Na 

1130  PRINT  "YIELD  STRENGTH  -":Ys:"  ULTIMATE  TS  -";Uts 

1140  PRINT  "ELASTIC  MODULUS  -";£:"  FLOW  STRESS  -";Fs 

1150  PRINT  •"' 
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1 7  60 

1170 
1180 
1190 
1200 
1210 
1220 
1230 
1240 
1250 
12S0 
1270 
1280 
1290 
1300 
1310 
1320 
1330 
DS:" 
1340 
1350 
1360 
1370 
1380 
1390 
1400 
1410 
1420 
1430 
t440 
1450 
T460 
1470 
1480 
1490 
1500 
1510 
1520 
1530 
1540 
1550 
1560 
1570 
1580 
1590 
1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
HP- 
1690 
1700 
1710 
1720 
1730 


PRINTER  IS  I 
t 

PRINT  CHR3<12) 
i 

!  CONVERSION  TO  METRIC 

Conien-25.4 

W»W*Conlen 

A-A»Conlen 

Bmax * Bmax *Co n I en 

Bnu  n» Bmi n»Con len 

Hs*Hs»Conlen 

Conpress-6. 89475293 

E=*E*Conpress 

Ys"Ys»Conpress 

Uts»Uts*Conpres3 

Fs*Fs*Conpress 

» 

PRINT  "CONVERTING  THIS 


!  TOGGLES  PRINTER  BACK  TO  CRT 


•  CLEARS  CRT 


UNITS 
!  1  in 


25.40mm 


!  Jpsi  -  6.39475293kPa  (HP  PROGRAMING  BOOK) 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
INPUT 


"UIDTH  •";«:" 
"Bmax  "":Bmax:" 
"HALF  SPAN  •":Hs 
"YIELD  STRENGTH 
"ELASTIC  MODULUS 


DATA  TO  METRIC  UNITS  (mm.N.KPa/ 

NOTCH  LENGTH  «":A 
Bmin  *";Bmm 
;"    POISSON'S  RATIO  -";Nu 
■"lYs:"  UTILMATE  TS  »";Uts 

-";E;"  FLOW  STRESS  »";Fs 


for  ";DS:"  YIEL 


'DO  YOU  WISH  TO  HAVE  A  PRINT  OUT  OF  THIS  DATA?" 


IF 
» 


GS-"N"  OR  GS-"N0"  THEN  1540 


PRINTER  IS  706  !  TOGGLES  TO  THERMAL  PRINTER 

PRINT  "CONVERTING  THIS  DATA  TO  METRIC  UNITS  (mm.N.KPa)  FOR: 
'WIDTH  -";W:"       NOTCH  LENGTH  »";A 
'Bmax  »":Bmax:"     Bmin  "":Bmin 
•HALF  SPAN  -";Hs:"    POISSON'S  RATIO  -";Nu 


PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT  " 

PRINTER 

PRINT  CHRSU2) 

RETURN 
» 


";DS 


"YIELD  STRENGTH  »";Ys:" 
"ELASTIC  MODULUS  »";E;" 

IS  1 


ULTIMATE  TS 
FLOW  STRESS 


;Uts 
;Fs 


!  TOGGLES  PRINTER  TO  CRT 
!  CLEAR  CRT 


♦  RETURNS  TO  k-KEY  SELECTION 


Setup _plotter:  »  SET  UP  INTERACTIVE  PLOTTER 
!  YOU  CAN  CHOOSE  BETWEEN  THREE  DIFFERENT  PLOTS 
!  <P  vs  COD.  P  vs  DEFLECTION.  OR  J  vs  DELTA  A) 


PRINT  CHRSU2) 
DIM  TS[801,XSt80],YS[80] 
GCLEAR 
GINIT 

GRAPHICS  ON 

PLOTTER  IS  3," INTERNAL" 
7225B  PLOTTER 


!  CLEARS  CRT 


!  PLOTS  ON  CRT  FOR  VIEWING  BEFORE  SENDING  TO 


PRINT  "TYPE  OF  PLOT:" 

PRINT  "    1 .  LOAD  vs  COD" 

PRINT  "    2.  LOAD  vs  DEFLECTION" 

PRINT  "    3.  J  vs  CRACK  EXTENSION" 

PRINT  "ENTER  THE  *  GF  THE  PLOT  YOU  DESIRE" 
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1740  INPUT  C2 

1750  PRINT  "" 

1760  PRINT  "" 

1770  IF  C2-1  THEN  1810 

1780  IF  C2  =  2  THEN  1860 

1790  IF  C2=3  THEN  1910 

1800  PRINT  "" 

1310  TS-"PLOT  of  LOAD  vs  COD  for  "ADS 

1820  PRINT  TS 

1830  XS="COD  (mm)" 

1840  Y$="L0AD  <N>" 

1350  GOTO  1950 

I860  TS="PL0T  of  LOAD  vs  DEFLECTION  for  "ADS 

1870  PRINT  TS 

1880  XS="DEFLECTI0N  (mm)" 

1890  Y$="L0AD  (N)" 

1900  GOTO  1950 

1910  TS="PL0T  of  J  vs  CRACK  EXTENSION  for  "ADS 

1920  PRINT  TS 

1930  XS-"DELTA_a  (mm)" 

1940  YS="J      <lcJ/nT2>" 

1950  PRINT   ••" 

1960  PRINT  "USING  METRIC  UNITS  (mm.N.ic J/m '2.etc>  :" 

1970  PRINT  "ENTER  HERE  THE  Xnun.Xmax  .  Ymin  .  Ymax" 

1980   INPUT  XKC2)  .X2(C2>.Y1(C2).Y2(C2) 

1990  X3=X2(C2)-XKC2) 

2000  Y3»Y2(C2)-YT<C2) 

2010  VIEWPORT  13. 100»RATI0. 10.90 

2020  WINDOW  X1(C2).X2(C2).Y1(C2).Y2<C2> 

2030  PRINT  "" 

2040  PRINT  "INPUT  DELTA_X  AND  DELTA_Y  LABEL  INTERVAL" 

2050  PRINT  "FOR  THE  TWO  AXES:" 

2060    INPUT  X4.Y4 

2070  PRINT  "" 

2080  PRINT  "DO  YOU  WISH-  TO  HAVE  A  GRID  OVERLAY?" 

2090  INPUT  G3 

2100  IF  GS="N0"  OR  G$-"N"  THEN  2140 

2110  PRINT  "" 

2120  PRINT  "INPUT  GRID  LINES  PER  DELTA  LABEL:" 

2130  INPUT  B0 

2140  PRINT  "PRESS  'ENTER'  TO  PLOT" 

2150  INPUT  H$ 

2T60  PRINT  CHRS(12)    !  CLEARS  CRT  < 

2170  IF  GS-"N0"  OR  GS-"N"  THEN  2310 

2180  FOR  I-X4/B0+XKC2)  TO  X2<C2)  STEP  X4/B0 

2190     MOVE  I.YKC2) 

2200     DRAW  I,Y2(C2) 

2210  NEXT  I 

2220  FOR  I»Y4/B0+YKC2>  TO  Y2<C2)  STEP  Y4/B0 

2230     MOVE  X1(C2),I 

2240     DRAW  X2(C2).I 

2250  NEXT  I 

2260  GOTO  2320 

2270  IF  G$-"Y"  OR  G$-"YES"  THEN  2320 

2280  MOVE  X2<C2).Y1(C2) 

2290  DRAW  X2(C2) . Y2(C2) 

2300  DRAW  X1(C2).Y2<C2> 

2310  CLIP  ON 

2320  AXES  X4.Y4.X1(C2>.YKC2> 

2330  CLIP  OFF 
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2340 
2350 
2260 
2370 
2330 
2390 
2400 
2410 
2420 
2430 
2440 
2450 
2460 
2470 
2480 
24S0 
2500 
2510 
2520 
2530 
254G 
2550 
2560 
2570 
2580 
2590 
2600 
2610 
2620 
2630 
2640 
2650 
2660 
2670 
2630 
2690 
2700 
2710 
2720 
2730 
2740 
2750 
2760 
2770 
2780 
2790 
2800 
2810 
2820 
2830 
2840 
2850 
2860 
2870 
2880 
2390 
2900 
2910 
2920 
2930 


VIEWPORT  I3.!00*RATI0.5.90 

WINDOW  X1 <C2> .X2(C2),Y1 <C2>,Y2<C2> 

LQRG  4 

CSIZE  3. .6        !  LABELS  X  DIGITS 

FOR  X-XKC2)  TO  X2(C2)  STEP  X4 

CLIP  OFF 

MOVE  X.YKC2) 

CLIP  OFF 

LABEL  USING  "K";X 
NEXT  X 

VIEWPORT  4.100-RATIO,!0,30 
WINDOW  XKC2)  ,X2(C2).YKC2)  .Y2(C2> 
LORG  2  !  LABELS  Y  DIGITS 

FOR  Y-YKC2)  TO  Y2<C2)  STEP  Y4 

CLIP  OFF 

MOVE  XKC2),Y 

CLIP  OFF 

LABEL  USING  "K";Y 
NEXT  Y 

VIEWPORT  O.RATIO«IOO, 0.100 
WINDOW  XI <C2> ,X2(C2> ,Y1 (C2> , Y2(C2) 
CCIZE  5 
LORG  6  !  PLOTS  TITLE 

MOVE  X2<C2)/2.Y2(C2> 

LABEL  TS 
DEG 

LDIR  90 
CSIZE  4.. 7 

MOVE  XKC2>,Y2<C2>/2 

LABEL  YS        !  PLOTS  Y-AXIS  LABEL 
LORG  4 
LDIR  0 

MOVE  X2<C2)/2.YKC2> 

LABEL  XS        !  PLOTS  X-AXIS  LABEL 
» 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT  "" 

PRINT  " 

PRINT  " 

PRINT  " 

INPUT  FS 

IF  F$«"N"  OR  F$ 


•  ••• 


IS  THE  PLOT  (AXES  AND  LABELS)  OK?" 

NOTE:   IF  YOU  ANSWER  YES  OR  Y" 
PLOT  WILL  BE  PRODUCED" 
ON  EXTERNAL  PLOTTER." 


NO"  THEN  GOSUB  Setup_plotter 
GCLEAR  !  CLEARS  PLOT  FROM  CRT 

GOSUB  Ext_piotter    !  SENDS  PLOT  TO  PLOTTER 


!  REDO  PLOTTER  SETUP 


PENUP 

MOVE  0,0 
i 

RETURN 
» 


!  RETURNS  TO  If -KEY  SELECTION  OR  R  curve  ROUTINE 


Ext_plotter: 
i 


!  DUMPS  PLOTTER  SETUP  TO  EXTERNAL  HP-7225B  PLOTTER 
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29<iu  jLj:r£R  is  ?05."hpgl" 

2950  VIEWPORT  13. 1 00-RATIO . 1 0 .30 

2960  WINDOW  XUC2>  .X2<C2>.YMC2>  .  Y2<C2> 

2970  IF  GS-'W  OR  GS-"N"  THEN  3030 

2980  FOR  I-X4/B0+XKC2)  TO  X2<C2>  STEP  X4/B0 

2990     MOVE  I.YKC2) 

3000     DRAW  I.Y2<C2> 

3010  NEXT  I 

3020  FOR    I-Y4/B0  +  YKC2)    TO   Y2(C2)    STEP    Y4/B0 

3030     MOVE  X1(C2),I 

3040     DRAW  X2(C2),I 

3050  NEXT  I 

3060  GOTO  3120 

3070  IF  6t-"Y"  OR  GS=MYES"  THEN  3120 

3080  MOVE  X2(C2),YKC2) 

3090  DRAW  X2(C2),Y2<C2> 

3100  DRAW  X1(C2),Y2<C2) 

31 10  CLIP  ON 

3120  AXES   X4,Y4.X1<C2>,YKC2> 

3130  CLIP   OFF 

3U0  VIEWPORT    13. T00-RATIO.5.90 

3150  WINDOW  X1(C2).X2(C2),Y1(C2) ,Y2(C2) 

3160  LORG  4 

3170  CSIZE  3.. 6  !  SIZE  OF  DIGITS 

3130  FOR  X-XKC2)  TO  X2(C2)  STEP  X4 

3' 90     CLIP  OFF 

3200     MOVE  X.YKC2) 

3210     CLIP  OFF 

3220     LABEL  USING  "K";X     !  LABELS  X-AXIS 

3230  NEXT  X 

3240  VIEWPORT  4, I00-RATIO. 10 .90 

3250  WINDOW  X1 <C2) ,X2(C2> , Y1 <C2> , Y2(C2) 

3260  LORG  2 

3270  FOR  Y-YKC2)  TO  Y2<C2)  STEP  Y4 

3280     CLIP  OFF 

3290     MOVE  X1(C2),Y  " 

3300     CLIP  OFF 

3310     LABEL  USING  "K";Y     !  LABELS  V-AXIS 

3320  NEXT  Y 

3330  VIEWPORT  0. 100»RATIO,0 .  100 

3340  WINDOW  X1(C2),X2(C2),Y1(C2).Y2(C2> 

3350  CSIZE  5  !  TITLE  SIZE 

3360  LORG  6 

3370     MOVE  X2(C2)/2.Y2<C2) 

3380     LABEL  TS  !  PRINTS  TITLE 

3390  DEG 

3400  LDIR  90 

3410  CSIZE  4,. 7 

3420     MOVE  X1(C2),Y2(C2)/2 

3430     LABEL  YS  !  PRINTS  Y-AXIS  TITLE 

3440  LORG  4 

3450  LDIR  0 

3460     MOVE  X2(C2)/2.Y1(CC) 

3470     LABEL  XS  !  PRINTS  X-AXIS  TITLE 

3480  PENUP 

3490  MOVE  0,0 

3500  ! 

3510  RETURN'        !  RETURNS  TO  Setup_plotter 

3520  ! 

3530  !  #«MM«HMMM«HMHWMMMM«HIWHI««MWHi 
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3540  ! 

3550  Imtiai_ai    !  ESTIMATES  INITIAL  CRACK  LENGTH 

3560  8etf-8max-<8max-Bmin>  2/Bmax     !  USES  B( EFFECTIVE)  FROM  ASTM  ESI  3 -81 

3570  D2-1  !  INITIALIZES  COUNTER 

3580  PRINT  CHRS<12>        f  CLEARS  CRT 

3590  PRINT  "BEFORE  rOU  CAN  RUN  THIS  ROUTINE  YOU  MUST:" 

3600  PRINT  "     1.  ENTER  SPECIMEN  DATA" 

3610  PRINT  "     2.  ENTER  CONVERSION  DATA"   . 

3620  PRINT  "     3.  SET  UP  PLOTTER" 

3630  PRINT  "     4.  CREATE  DATA  FILE" 

3640  PRINT  "" 

3650  PRINT  "ARE  YOU  READY?" 

3660  INPUT  HS 

3670  IF  HS-"N"  OR  HS-"NO"  THEN  5070   !  RETURNS  TO  KEY  SELECTION  SCREEN 

3680  PRINT  CHRS(12>  !  CLEARS  CRT 

3690  PRINT  "THIS  ROUTINE  WILL  OBTAIN  THE  INITIAL  CRACK  LENGTH  AFTER  FATIGUE  PRE 

-CRACKING." 

3700  PRINT  — 

3710  PRINT  "YOU  MAY  LOAD  AND  UNLOAD  THE  SPECIMEN  AS  OFTEN  AS   YOU  LIKE  —  BUT  Y 

OU  MUST  REMAIN  "; 

3720  PRINT  "ON  THE  LINEAR      PORTION  OF  THE  CURVE!" 

3730  PRINT  "" 

3740  PRINT  "NOTE:   THE  MIS  SHOULD  BE  IN  'STROKE'  CONTROL" 

3750  PRINT  "       AT  A  VERY  SLOW  INVERTED  RAMP  FUNCTION" 

3760  PRINT  "       RATE  OR  YOU  MAY  DO  THIS  WITH  MANUAL  LOAD" 

3770  PRINT  "       CONTROL." 

3730  PRINT  "" 

3750  ! 

3800  PRINT  "WOULD  YOU  LIKE  AN  ESTIMATE  OF  Pma*  <Lbs>7" 

3810  INPUT  ES 

3320  IF  ES-"N"  OR  E3-"N0"  THEN  4060   !  SKIPS  Pmax  CALC 

3830  ! 

3840  PRINT  "YOU  MUST  PROVIDE  AN  ESTIMATE  OF  CURRENT  CRACK      LENGTH" 

3850  PRINT  "INPUT  YOUR  ESTIMATE  IN  INCHES." 

3360  INPUT  Aest 

3870  !  CONVERT  ALL  DIMENSIONS  BACK  TO  ENGLISH  UNITS. 

3880  !  aO  IS  SUFFIX  ON  ALL  VARIABLES  FOR  THIS  CALC. 

3f.30  BeffaO-BefF/Conlen 

3900  WaO-W/Conlen 

39t0  Blig-<WaO-Aest>         !  ESTIMATED  LENGTH  OF  UNBROKEN  LIGAMENT  IN  INCHES 

3920  FsaU-Fs/Conpress  9 

3930  HsaO-Hs/Conlen 

3940  PlaO-.533»BeffaO»Blig*2*FsaO/(2»HsaO) 

3950  PRINT  "AN  ESTIMATE  OF  Pmax  IS";PlaO 

3960  PRINT  "Pmax  ESTIMATES  THE  POINT  AT  WHICH  THE  CURVE" 

3970  PRINT  "BECOMES  NONLINEAR,  IAW  ASTM  E813-81" 

3980  ! 

3990  PRINTER  IS  706   !  TOGGLES  PRINTER 

4000  PRINT  "" 

4010  PRINT  "GIVEN  a»" ; Aest :" inches  THEN:" 

4020  PRINT  "   Pmax-"  ;Pla0;  "Lbs  FROM  ASTM  E813-81" 

4030  PRINT  '•" 

4040  PRINTER  IS  1 

4050  ! 

4060  DIM  Lod(8000),Def K8000) ,Cod(8000) 

4070  ! 

4080  !  ASSIGN  DEFL  CORRECTION  COEF  USED  IN  Data_acq  AND  Take_data  ROUTINES 

4090  !  THIS  USES  A  4th  ORDER  POLYNOMIAL  FIT  OF  RECORDED  DATA  ON  THE  MTS  SYSTEM' 

S  COMPLIANCE 

4100  ! 
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4  I  !  (j 

4120 

4130 

4140 

4150 

4160 

4T70 

4180 

4190 

4200 

4210 

4220 

4230 

4240 

4250 

4260 

4270 

4280 

4290 

4300 

4310 

4320 

4330 

ARK 

4340 

4350 

4360 

4370 

4380 

4390 

4400 

4410 

4420 

4430 

4440 

4450 

4460 

4470 

4480 

4490 

4500 

4510 

4520 

4530 

4540 

4550 

4560 

4570 

4580 

4590 

4600 

4610 

4620 


tor r* I > 

Corr<2> 

Corr(3) 

Corr<4) 

Corr(5> 
« 

PRINT 

PRINT 

INPUT 

PRINT 

K-0 

I-l 

Flag! 


■5.798E-5 

■4.7419E-6 
■-2.803E-9 
'1  .471E-12 
-2.7741E-I6 


"WHEN  YOU  ARE  READY  TO  BEGIN  PRESS  THE  ENTER'  KEY 

HS 

CHRSC12) 

!  USED 

!  USED 
0    !  NOTE; 


!  CLEARS  CRT 
AS  A  COUNTER 
AS  INDEX 
FLAG1  AND  FLAG2 


ARE  RESET 


ESTABLISHES  ONE  LOAD-UNLOAD  DATA  TAKING  LOOP 


!   RESETS  FLAG2 


Flag2-0    !  IN  Data_acq  ROUTINE 
WHILE  Flagl-0 

GOSUB  Data_acq 
END  WHILE 
FIag2-0 
WHILE  Flag2»0 

GOSUB  Data_acq 
END  WHILE 
PENUP 

!  THIS  SECTION  CAN  BE  USED  TO  VIEW  THE  COLLECTED  DATA  BY  REMOVING  THE  REM 
CHARACTER  <!>. 


!  FOR  I-l  TO  K 

!  PRINT  USING 

!  IMAGE 

!  NEXT  I 

!  WAIT  5 

KI  - 1 

K3-K 

GOSUB  Slope 

CorrelO(D2) 
i 


3780 
X.3D.4X.6D 


■Correl 


I.Lod(I)  .Def  KD.CodCD 
,D,4X.3D.7D.7X,3D.7D 


SLOPE 
SAVES 


•  J /COD 
CURRENT 


COMPLIANCE 
CORRELATION 


COEFFICIENT  IN  AN  ARRAY 


Compl0<D2)-Corapl 
i 

GOSUB  Crack_length 
» 

A0(D2)-Crckl 
\ 

!      FIND   957.  CONFIDENCE    INTERVAL  FOR   aO 
Stiff-Stiff+T9»U4 
GOSUB  Crack_length 
AmaxO<D2)-CrcH 
Stiff-Stiff-2*T9»U4 
GOSUB   Craclc_length 
AnnnO<D2)-Crckl 


!    SAVES  CURRENT   COMPLINACE  VALUE    IN   AN   ARRAY 


!    USES  COD   COMPLIANCE 

!    SAVE   THE  LAST  aO   MEASUREMENT   OBTAINED    IN  AN   ARRAY 


j 

PRINT 
PRINT 
PRINT 
PRINT 


CHRSU2) 

"THE  FOLLOWING 


!  CLEARS  CRT 
INFORMATION  HAS 


BEEN  OBTAINED: 


NT  #";D2 


'INITIAL  CRACK  LENGTH  <a0)  -"  :  A0<D2> : "mm 


FROM  MEASUREME 


4630 
4640 

mm" 
4650 
4660 
4670 


PRINT 
PRINT 

PRINT 
PRINT 
PRINT 


"WITH  95%  CONDIFENCE  INTERVAL  ON  a  FROM"; Ami nO<D2) ;"  TO" ;AmaxO<D2) 

Iff** 

"USING  COD  COMPLIANCE  OF  :  "  ;ComplO(D2)  :"n>m/N" 

"WITH  A  CORRELATION  OF :" ;CorrelO(D2) : "FROM" :K; "DATA  POINTS." 
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4630 

4690 

4700 

4710 

4720 

4730 

4740 

4750 

4760 

4770 

4780 

4730 

4800 

4810 

4320 

4830 

4840 

4850 

4860 

4870 

4830 

4890 

TORE. 

4900 

4910 

4920 

4930 

4940 

4950 

4960 

4970 

4980 

4950 

5000 

5010 

5020 

5030 

5040 

a  rs 

5050 
5060 
5070 
5080 
5090 
5100 
5110 
5120 
5130 
5140 
5150 
5160 
5170 
5180 
5190 
5200 
5210 
5220 
5230 
5240 


PRINT  ••" 

PRINT  "DO 
INPUT  HS 
IF  HS-"N" 


PRINTER 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT  " 

PRINTER 
i 

D2-D2+1 
PRINT  "" 
PRINT  "DO 
INPUT  GS 

IF  GS-"Y" 
t 

PRINT  "INPUT 


YOU  WANT  A  HARD  COPY?" 
OR  HS-"NQ"  THEN  4830 

TOGGLES  PRINTER 


IS  706 


"THE  FOLLOUING 
"INITIAL  CRACK 


» 


INFORMATION  HAS  BEEN 
LENGTH  -":A0<D2>  :"min 


OBTAINED:" 

FROM  MEASUREMENT  *";D2 


"WITH  957.  CONFIDENCE  INTERVAL  FROM" : Ami n0<D2) 
"USING  COD  COMPLIANCE  OF : " :Comp 10< D2> ; "mm/N" 
"WITH  A  CORRELATION  OF : " ;CorreiO(D2) : "FROM" ;K 


"  T0";AmaxG<D2> ;"mn' 
"DATA  POINTS." 


IS  1 


!  TOGGLES  PRINTER 

!  INCREMENT  COUNTER 


YOU  WANT  ANOTHER  aO  MEASUREMENT?" 
OR  G3«"YES"  THEN  4060 


THE  NUMBER  <#>  OF  THE  aO  MEASUREMENT  YOU 


UISH  TO  USE  AND 


'YOU 


PRINT  "    NOTE: 
PRINT  " 

INPUT  D2 
; 

A1-A0<D2> 
!  NOTE:  A1 
f 

PRINT  "" 

PRINT 

PRINT 

INPUT 

♦  NOTE 
t 

IF  PS-"N"  OR 

K2-0 

STORED 

GOfvJB  Stor^ 
» 

PRINT  CHRSC12) 

f 

j 

RETURN 


ONLY  THIS 
STORED  ON 


SET  OF 
DISK." 


DATA  ARE" 


PS 


!  ASSIGNS  A1  TO  DESIRED  A0(D2)  VALUE 
IS  USED  IN  Mam_prog  AS  INITIAL  CRACK  LENGTH  VALUE 


HAVE  SELECTED  MEASUREMENT  *";D2 
ARE  YOU  SURE?" 


ONLY  THIS  aO  MEASUREMENT  IS  STORED  ON  DISK 

PS-"NO"  THEN  4890    •  ♦  ALLOWS  CORRECTION 

!  K2  IS  UNLOADING  COUNTER  THIS  WILL  ENSURE  CORRECT  DAT 


data 


-*••  WRITES  ON  DISK 


!  RETURNS  TO  Ic-KEY  SELECTION 


♦  DATA  ACQUISITION  ROUTINE 
Dvm;"VR5  AF1  AL3" 


READS  LOAD 


Data_acq 
OUTPUT 
f 

OUTPUT  Dvm:"AIt  VT1" 
ENTER  Dvm:Lod<I) 
Lod<I)-Lod(I)»SK1) 
t 

OUTPUT  Dvm;"AI2  VT1" 
ENTER  Dvm:Cod<I) 
Cod(I)-<Cod<I))»SI<2) 
!  NOTE:  COD  IS  CALIBRATED 


READS  COD 


ONLY  FORM  0.475  TO  0.525  INCHES  (MAX  COD  IS  0 . 1 


5  INCHES) 
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5250  ! 

5260  OUTPUT  Dvm:"AI3  VT1"     !  READS  DEFLECTION 

5270  ENTER  Dvm:Defl< I> 

5230  Def  i<I)-Oef  KD-SK3) 

5290  ! 

5300  GOSUB  Correct_defl      \    DETERMINE  DEFL  CORRECTION  FACTOR 

5310  ! 

5320  DefKD-DefKD-Totcor 

5330 

5340 

5350  !     PLOT  THE  DATA  POINT 

5360  VIEWPORT  13 . 1 00»RATIO . 1 0 .90 

5370  UINDOW  X! <C2> .X2(C2) , Yl (C2) , Y2<C2> 

5380  IF  C2-2  THEN  5410 

5390  PLOT  Cod<I),Lod<I) 

5400  GOTO  5420 

5410  PLOT  Defl<I),Lod<I> 

5420  ! 

5430  !  TAKE  DATA  FOR  ONE  COMPLETE  LOADING  AND  UNLOADING 

5440  IF  1-1  THEN  5520         !  SKIPS  TEST  ON  THE  FIRST  DATA  POINT 

5450  ♦ 

5460  !  TEST  DATA 

5470  Lodiff-Lod<I)-Lod(I-1) 

5480  IF   Lodiff<-25.0    THEN  Flagl-1         !    MUST   BE          UNLOADING 

5490  IF   Locuff>25.0    THEN  Flag2'1           !    MUST   BE        LOADING 

5500  Lodiff-O. 

5510  ! 

5520  I-I+t      !  INCREMENTS  INDEX 

5530  ! 

5540  K-K  +  1      !  INCREMENT  COUNTER 

5550  ! 

5560  RETURN             !    RETURNS   TO   Imtial_a  ROUTINE 

5570  ! 

5580  !  .——.»—»••—•———•*•»—.•**.••• 

5590  Correct_defI:      !  THIS  ROUTINE  CALCULATES  THE  DEFLECTION  CORRECTION  FACT 

OR 

5600  ! 

5610  !  CORRECTION  FACTOR  IS  REQUIRED  TO  ACCOUNT  FOR  SYSTEM  COMPLIANCE 

5620  !  A  4TH  ORDER  POLYNOMIAL  FIT  HAS  PREVIOUSLY  DETERMINED. 

5630  !                                                   f 

5640  !     NOTE:   •••••••••••••••••••••• 

5650  !  THE  5  CORRECTION  COEFFICIENTS  CORR(I)  HERE  ASSIGNED  IN  Irutial_a  ROUTIN 

E  TO  SPEED  UP  PROG. 

5G60  !       ••••••••••••—••••...•..••• 

5670  ! 

5680  Lodp-Lod<I)/Conlb    !  CONVERTS  CURRENT  LOAD  TO  Lbs  FOR  CALC 

5690  Defied  )-Lodp 

5700  FOR  J-2  TO  4 

5710  Def lc(J)-Def lc(J-1)»Def lc<1> 

5720  NEXT  J 

5730  Totcor-Corr<1)      ♦  PREVIOUSLY  ASSIGNED  AS  NOTED. 

5740  FOR  J*1  TO  4 

5750  Totcor-Totcor+Corr<J+1)»Def Lc(J> 

5760  NEXT  J 

5770  ! 

5780  Totcor-Totcoi-Conien    !  CONVERTS  TOTAL  CORRECTION  FACTOR  FROM  (in)  TO  < 

mm) 

5790  ! 

5800  ! 

5810  RETURN             !    RETURNS    TO   Data  acq   OR    Take_data  ROUTINE 
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5820  f 

5830    !••«••••• .•—.. •*•*••••««• 

5840  ! 

5850  Craclc_length:  !  CALCULATES  CRACK  LENGTH  USING  HUDAK  (DTNSRDC  PROG)  [PROF.  J 

.  JOYCE'S  EQUATION] 

5860  !   NOTE:  THIS  ROUTINE  USES  COD  COMPLIANCE  METHOD  VICE  DEFLECTION  COMPLIA 

NCE. 

5370  FOR  1-1  TO  5 

5880  READ  Coef(I) 

5890  NEXT  I 

5900  Comol-1/Stiff 

5910  C»2«U*8eff»E»Compl/Hs/1000 

5920  C-1/(SQR(C)+1 > 

5930  G<1)«C 

5940  FOR  1-2  TO  5 

5950  G(I)-G(I-1)»G(T) 

5960  NEXT  I 

5970  F-Coef(l) 

5980  FOR  1-1  TO  5 

5990  F-F+Coef(IH)-G(I) 

6000  NEXT  I 

6010  Crcki-F-U 

6020  RESTORE 

6030  DATA  0 . 998265,-3.81662 , -  1 .30596 ,32 . 31 033 , -44 . 1 5665 , -52 .67876 

6040  ! 

6050  RETURN                  !    RETURNS   TO    Imtial_a  OR   Main_prog 

6060 

6070 

6080 

6090  Conv_data:    !  INPUTS  CONVERSION  DATA  FOR  MEASURING  DEVICES 

6100  PRINT  CHRS(12)          !  CLEARS  CRT 

6110  PRINT  "INPUT:  LOAD  RANGE.  STRAIN  RANGE.  STROKE  RANGE" 

6120  PRINT  "    (ENTER  IN  PERCENT  ie.  10,20,50,100)" 

6130  INPUT  Lr.Cr.Sr 

6140  PRINT  "" 

6150  PRINT  "LOAD  RANGE  -  ";Lr 

6160  PRINT  "STRAIN  (COD)  RANGE  -  ":Cr 

6170  PRINT  "STROKE  (DEFLECTION)  RANGE  -  ";Sr 

6180  Lmax-20000     !  MAX  LOAD  RANGE  IN  (lbs). 

6190  Craax-.200      !  MAX  COD  RANGE  IN  (in). 

6200  Smax-5.0       !  MAX  STROKE  RANGE  IN  (in). 

6210  V-10          !  MAX  VOLTAGE  IN  (volts). 

6220  ! 

6230  !   MAXIMUM  OPERATING  POINTS  ARE: 

6240  Olmajt-Lmajc-Lr/100 

6250  Ocmajc-CmajfCr/100 

6260  Osnax-SmajfSr/lOO 

6270  ! 

6280  !  CONVERSION  CONSTANTS  ARE:  (ENGLISH  UNITS) 

6290  Sle(1)=Olmax/V 

6300  Sle(2)-0craajt/V 

6310  Sle(3)-0smax/V 

6320  !  CONVERTING  TO  METRIC: 

6330  Conlen-25.4            !  25.4mm  -  1  in 

6340  Conib-4. 448221615      !  4.44822T615N  -  1  lb      (HP  PROGRAMING  BOOK) 

6350  ! 

6360  !  Sled, 2. 3)  IS  USED  TO  DENOTE  CONVERSION  OF  VOLTAGE  TO  ENGLISH  UNITS. 

6370  Sl(1)-Sle(1)»Conlb«(-1 .0)    !  COMPRESSIVE  LOAD  IS  NEGATIVE  VOTLAGE 

6380  Sl(2)-Sle(2)»Conien 

6390  SK3)-Sle(3)»Conlen«(-1  .0)    !  COMPRESSIVE  STROKE  SUPPLIES  A  NEGATIVE  VOL 

TAGE 
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5400  ! 

6«4:u  PRINT  "" 

6^20  PRINT  •"' 

6430  PRINT  "YOUR  ENTRIES  HAVE  PRODUCED  THE  FOLLOWING         CONVERSION  DATA: 

6440  PRINT  "" 

6450  PRINT  "FOR  LOAD:  ":SK  I)  : "NEWTON/VOLT   ("  ;Sle<  I )  :"Lbs/V>" 

6460  PRINT  '- 

6470  PRINT  "FOR  COD:  " ;Si <2> ;"mm/V  <M;Sle<2> :"in/V>" 

6480  PRINT  "M 

6490  PRINT  "FOR  DEFLECTION:  ";SM3>  ;"m/V   <"  ;Sle(3) :  "in/V)" 

6500  PRINT  "M 

6510  PRINT  "WOULD  YOU  LIKE  A  HARD  COPY?" 

6520  INPUT  HS 

6530  IF  HS-"N"  OR  HS-"NO"  THEN  6630 

6540  ! 

6550  PRINTER  IS  706  !  TOGGLES  PRINTER 

6560  PRINT  "" 

6570  PRINT  "       CONVERSION  DATA" 


"Lbs/V)" 

"in/V)" 

Sle<3);"in/V) 


6580  PRINT  "FOR  LOAD :" ;S1< 1 ) ;"N/V   <";Sle<1> 

6590  PRINT  "FOR  COD: " :S1 (2) :"mm/V    <":SIe<2> 

6600  PRINT  "FOR  DEFLECTION  :  "  ;SL  (3)  :"n.m/V   <" 

6610  PRINTER  IS  i            !  TOGGLES  PRINTER 

6620  ! 

6630  PRINT  CHRSU2) 

6640  ! 

6650  RETURN         !  RETURNS  TO  k-KEY  SELECTION 

6660 

6670 

6630 

6690  Slope:    !  LINEAR  LEAST  SQUARES  FIT 

6700  PRINT  CHRS< 12)        !  CLEAR  CRT 

6710  !  INITIALIZE  VARIABLES 

6720  P2-0       !  USED  AS  A  COUNTER 

6730  Sumx-0     !  SUM  OF  X  VALUES 

6740  Smnx2-0    !  SUM  OF  X  SQUARED  VALUES 

6750  Sumv-0     !  SUM  OF  Y  VALUES 

6760  Sumy2-0    !  SUM  OF  Y  SQUARED  VALUES 

6770  Sumxy-0    !  SUM  OF  THE  PRODUCT 

6780  FOR  I-K1  TO  K3    !  KI-INDEX  OF  FIRST,  K3-INDEX  OF  LAST  DATA  POINT  USED 

6790  !  TYPE  OF  PLOT 

6800  IF  C2-1  THEN  6840         !  LOAD  vs  COD 

6810  IF  C2-2  THEN  6880         !  LOAD  vs  DEFL 

6320  .     IF  C2-3  THEN  6920         !  J  vs  DELTA_a 

6830  ! 

6840  X»Cod<I) 

6850  Y-Lod(I) 

6860  GOTO  6990 

6870  ! 

6880  X-DefKI) 

6890  Y-Lod<I) 

6900  GOTO  6990 

6910  » 

6920  X-Delta_a<I) 

6930  Y-J9<D 

6940  » 

6950  !  CHECK  IF  DATA  POINT  IS  WITHIN  BOUNDS 

6960  IF  Delta_a<IXJ9<I>/(2»Fs>  +  .15  THEN  7050     !  CHECKS  LOWER  LIMIT 

6970  IF  Deita_a<I)>J9<I)/<2»Fs>+1.5  THEN  7050     !  CHECKS  UPPER  LIMIT 

6980  ! 

6990  Swnx-Sujnx+X 
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7H0U 
7010 
7020 
7030 
7040 
7050 
7060 
7070 
7080 
7090 
7100 
7110 
7120 
7130 
7140 
7!50 
7160 
7170 
7180 
7130 
7200 
7210 
7220 
7230 
7240 
7250 
7260 
7270 
7280 
7290 
7300 
7310 
7320 
7330 
7340 
7350 
7360 
7370 
7380 
7390 
7400 
7410 
7420 
7430 
7440 
7450 
7460 
7470 
VER 


GET  NEXT  DATA  POINT 
INCREMENTS  COUNTER  EACH  TIME 


Sumy  •Sumy  +  T' 

Sumx2-Sumx2+X»X 

S  am  y  2  *  S  urn  y  2  ♦  Y  *  Y 

SuLmxy"Sumxy*X»Y 

GOTO  7060      ! 

P2-P2+1        ! 
NEXT  I 

K3-K3-P2      !  USED  TO  CALCULATE  NUMBER  OF  OATA 
K4-K3-K1+1        !  K4-T0TAL  *  DATA  POINTS  USED 
Qal  -Sumx2-Sumx»Sumx/K4 
Qa2*Sumxy-Sumx»Sumy/K4 
Qa3-Sumy2-Sumy*Sumy/K4 
» 

M-Qa2/Qal  !  SLOPE 

r 

IF  M<0  THEN  7430   !  IF  SLOPE  <  0 ,  BAD  DATA 

!   FROM  DTNSRDC  PROGRAM 

G1»<K4«Sumx2-Sumx"2)/(K4«<K4-1)) 

Q2-(K4»Sumy2-Sumy',2>/(K4»<K4-1  )> 

G3-<K4-t )»<Q2-M-M»01 )/(K4-2) 

Q4-SQR<Q3/<Q1*<K4-1)>) 

» 

B-<Sumy-M-Sumx)/K4     !  INTERCEPT 
i 

Correl; 


A  POINT  IS  ELIMINATED 

POINTS  UITHIN  BOUNDS. 
IN  CALC 


M*SQR<Ql/Q2) 


Qa4«SGR<ABS<  <Sumy2-B»Sumy-M»Sumxy>/(K4-2) ) ) 

Sigmab-Qa4/SQR<ABS<Qa1  >  > 

S i gmam-S igmao*SOR < Sumx2/K4 ) 

Stiff-M  !  STIFF  -  SLOPE 

Compi-1/M        !  COD  COMPLIANCE  -  1 /SLOPE 

; 

!  STUDENT  t  TEST 

U4-Q4 

T9-1.96 

IF  K4>120  THEN  7410  !  RETURNS 

T9-2.1 

IF  K4<20  THEN  7410  !  RETURNS 

N9"K4-2 

T9-N9/<-.617109+N9-. 510239)  <, 

t 

RETURN 
♦ 

PRINT  "POOR  DATA—  SLOPE  <  0" 
M-333333.3 


Correl-I 
» 

RETURN 
CALLED  IT) 


!  RETURNS  TO  ImUal.a.  Main_prog  or  R_curve  ROUTINE  (WHICHE 


file:     !  ESTABLISHES  DATA  FILE 
"THE  FILE  YOU  CREATE  CAN  BE  READ  BY 


THE 


•J  DATA"  PROGRAM 


7480 

7490 

7500 

7510  Data_ 

7520 

7530 

7540   !  NOTE:  ANY  CHANGES  MADE  TO  THIS  ROUTINE  MAY  REQUIRE  SIMILAR  CHANGES  TO 

J_DATA"  PROGRAM  DUE 

7550    !  TO  DISK  FORMATTING  CONSULT  HP  MANUAL  FOR  ASSISTANCE. 

7560    ! 

7570        PRINT  CHRSU2)  !    CLAERS   CRT 
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7580   PRINT  "THIS  ROUTINE  HILL  CREATE  A  DATA  FILE  GN  DISK  FOR  THE  TEST." 

7590   PRINT  ••" 

7S00   PRINT  "INSERT  AN  INITIALIZED  DATA  DISK  AMD  ENTER  FILE      NAME  (SAMPLE  I 

D  IS  SUGGESTED >." 

7GI0    INPUT  NameS 


7520 

PRINT 

'••• 

7B30 

PRINT 

"IS  THIS  A  NEW  FILE?" 

7640 

INPUT 

OS 

7650 

IF  QS- 

"N"  OR  Q$-"NO"  THEN  7630 

7660 

I 

7670 

CREATE 

BDAT  NameS. 6, 40000 

7630 

ASSIGN 

«File  TO  NajneS 

7690 

» 

7700 

PRINT 

PtM 

7710 

PRINT 

"FOR  RECORD  PURPOSES  THE  FOLLOWING  DATA  IS  ALSO      RECORDED:" 

"720 

PRINT 

••  •• 

7730 

PRINT 

"ENTER  THE  DATE  (MM-DAY-YEAR) :" 

.'740 

INPUT 

DataSd) 

7750 

PRINT 

MM 

7760 

PRINT 

"ENTER  THE  TIME  (HH:MM):" 

7770 

INPUT 

DataS<2) 

7780 

PRINT 

"" 

7730 

PRINT 

"ENTER  THE  OPERATOR'S  NAME:" 

7800 

PRINT 

(UP  TO  10  CHARACTERS)" 

7810 

INPUT 

DataS<3) 

7820 

PRINT 

••" 

7830 

PRINT 

"ENTER  SPECIMEN  ID:" 

7840 

PRINT 

(UP  TO  10  CHARACTERS)" 

7850 

INPUT 

DataS(4) 

7860 

!  URITE  ON  DISK 

7870 

FOR  I- 

1  TO  4 

7880 

OUTPUT  *File;DataS(I)[1 ;10] 

7890 

NEXT  I 

7900 

• 

7910 

BEEP 

7920 

PRINT 

MM 

7930 

PRINT,- 

"DATA  FILE  HAS  BEEN  CREATED." 

7940 

; 

"950 

; 

"960 

RETURN 

!  RETURNS  TO  Ic-KEY  SELECTION 

7970 
7980 
7990 

9000 

Stor  data:    !  WRITES  DATA  ON  DISK 

8010 

» 

■ 

8020 

t  THE 

ITEMS  STORED  BY  THIS  ROUTINE  CAN  BE 

8030 

!  RECALLED  USING  THE  "J  DATA"  PROGRAM 

8040 

!  CARE 

.  MUST  BE  TAKEN  IF  ANY  CHANGES  ARE  MADE  TO  THIS  ROUTINE  SUCH  THAT  SI 

MILAR 

CHAGES 

ARE 

3050 

!  MADE 

.  TO  THE  "J  DATA"  PROGRAM. 

8060 

f 

3070 

IF  K2>0  THEN  8160        !  K2  IS  UNLOADING  COUNTER 

8080 

OUT 

PUT  (5Fiie;A0(D2),Amin0(D2),Amaj(0(D2).Compl0(D2)  ,CorrelO(D2)  .K 

8090 

t 

8100 

PRINT 

MM 

8110 

PRINT 

MM 

8120 

D2 

8130 

PRINT 

"INITIAL  CRACK  LENGTH  (aO)  INFORMATION  FROM        MEASUREMENT  *"; 

PRINT 

"IS  NOW  IN  DISK  FILE  ";NameS 

8140 

» 
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3150        IF   K2-0    THEN  8220  !    SKIP    IF    THIS   UAS      aJ   MEASUREMENT 

3160  OUTPUT   <?File:K2,AMK2>.Delta_a<K2>  .K.L.Correl  .Compl  .Hmm.Amax  .Area9<K2 

)    J9( K'*) 

3170  OUTPUT    3File:MdefKK2),Codiff  ,K0,KI 

8180 

3190 

3200 

8210 

8220 

3230 

8240 

3250 

8260 

8270 

3280 

3290 

3300  Take  data: 


PRINT  "" 

PRINT  "UNLOADING  INFORMATION  IS  NOW  IN  DISK  FILE  ";NameS 

» 

BEEP 

WAIT  3     !  ALLOWS  A  PAUSE  SO  MESSAGE  IS  PRESENTED  ON  THE  CRT 


RETURN 


!  RETURNS  TO  Imtial_a  OR  Main_prog 


!  THIS  ROUTINE  IS  USED  TO  TAKE  DATA  DURING  THE  ACTUAL  TEST, 


3310 

3220 

3330 

8340 

3350 

3360 

8370 

8380 

8390 

8400 

8410 

8420 

8430 

8440 

3450 

NCE 

8460 

3470 

8480 

3490 

8500 

3510 

3520 

8530 

3540 

3550 

8560 

8570 

8580 

8590 

8600 

8610 


r-K 

OUTPUT  Dvm:"AI1  VTV 
ENTER  Dvm:Lod(I) 
Lod(I>-Lod(I>»Si(1> 
{ 

OUTPUT  Dvm:"AI2  VT!" 
ENTER  Dvm;Cod(I) 
Cod(I>«Cod(I)«SK2> 
t 

OUTPUT   Dvm;"AI3  VTT" 
ENTER   Dvm:DefKI> 
Defl(I>-Defl(I)«Sl(3> 
» 

GOSUB  Correct   defl 


!    INDEX 

!    READS   LOAD 


!    READS  COD 


!    READS  DEFLECTION 


!    DETERMINE   CORRECTION  FACTOR   FOR  SYSTEM  COMPLIA 


Def  KD-OefKD-Totcor 
♦ 

!        PLOT   THE  DATA  POINT 
VIEWPORT    13. 100-RATIQ, 10,90 
WINDOW   XI  <C2)  ,X2(C2>  ,  YKC2)  ,Y2(C2> 
IF  C2-2   THEN  8550 

PLOT  Cod<I),Lod<I> 
GOTO   8560 

PLOT   DefKI).Lod<I> 
; 
j 

RETURN     !  RETURNS  TO  Nain_prog 
j 

f  MMMWMHNMMWWMMWHHHHHMH 


8620  J_caic: 

TNSRDC) 

8630    ! 

8640 

8650 

8660 

8670 

8680 

3690 

3700 

8710 


!  CALCULATES  J  WITH  CRACK  GROWTH  CORRECTION  ERNST  MODEL  (USED  BY  D 


E5-2 

♦   CALCULATE  J  FOR  FIRST  (1st)  UNLOADING 

IF  K2>1  THEN  8690 

J6-J7+E5/(W-Al(K2))-Area9(K2> /Brain 

GOTO  8730 

JS-J7*E5/(W-Al(K2)WArea9(K2)-Area9(K2-1>>/Bmin 

!  IF  DELTA_a  <  0.254mm  DO  NOT  USE  CRACK  GROWTH  CORRECTION 

IF  Delta  a<K2X.254  THEN  8730 
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8720   J6-J6*<!-1/<W-AKK2)>»<A1<K2)-A1<K; 
STM  E813-81) 


-I)))    !  CORRECTS  FOR  CRACK  GROWTH  <A 


8730 

8740 

8750 

8760 

8770 

8780 

3790 

8800 

8810 

8820 

8830 

8840 

8850 

EDING 

8S60 

8870 

8380 

8890 

8900 

8910 

8920 

8930 

8940 

8950 

89G0 

8970 

8980 

8990 

9000 

OKE* 

9010 

9020 

9030 

9040 


J7=J6 
i 

J9(K2>= 

RETURN 
i 


J6        !  STORES  CURRENT  J  VALUE  IN  J9  ARRAY. 
!  RETURNS  TO  Main_prog 


Main_prog:   !  THIS  IS  THE  MAIN  PROGRAM  FOR   THE  J-INTEGRAL  TEST, 


PRINT  CHRSU2)      !  CLEARS  CRT 

PRINT  "" 

PRINT  "ENSURE  YOU  HAVE  COMPLETED  THE  FOLLOWING  ITEMS 


BEFORE  PROCE 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
INPUT 


1.  ENTERED  SPECIMEN  DATA." 

2.  ENTERED  CONVERSION  DATA." 

3.  SET  UP  INTEACTIVE  PLOTTER." 

4.  OBTAINED  INITIAL  CRACK  LENGTH." 

p 

ARE  YOU  READY?" 


HS 

IF  HS="NO"  OR 
PRINT  "" 


HS="N"  THEN  10530   !  RETURNS  TO  WAIT  FOR  k-KEY  SELECTION 


PRINT 
PRINT 
INPUT 

Lod5  = 
t 

PRINT 
CONTROL 
PRINT 
PRINT 
PRINT 
PRINT 

TS  MACHINE." 

9050 

9060 

9070 

9080 

9090 

9100 

9110 

9120 

9130 

9140 

9150 

9160 

9170 

9180 


"WHAT  IS  YOUR  EXPECTED  MAXIMUM  LOAD  IN  NEWTONS?" 

Pmax 

l»Pmax     !  USED  TO  SIGNAL  TURN  AROUND  FROM  UNLOADING 

;;the  test  should  START  WITH  NEAR  zero  <0>  load  WITH  THE  MTS  IN  'STR 

"ensure  this  is  done." 

"when  you  are  ready  press  'enter'  or  'continue'  key  and  start  the  m 


INPUT  HS 
PRINT  CHR$(12> 


!  CLEARS  CRT 


INITIALIZE  VARIABLES 
Area-0    !  AREA  UNDER  LOAD  vs  DEFLECTION  CURVE   DTNSRDC  (18) 
Area9(1)-0  *  AREA  BENEATH  SPECIFIC  LOD  vs  DEFL  (DTNSRDC  I9(K2)> 


K=K 

I-K 

K0=K 

K1=0 

K2»0 

K3  =  0 

L-1 


CURRENT  INDEX  (AFTER  OBTAINING  INITIAL  CRACK  LENGTH) 

INDEX  OF  FIRST  TEST  DATA  POINT 

WHERE  TEST  DATA  STARTS 

LOWER  INDEX  FOR  SLOPE  CALC 

A  COUNTER      "    (K2)   NUMBER  OF  UNLOADINGS 

UPPER  INDEX  FOR  SLOPE  CALC 


A  COUNTER   DTNSRDC  (L) 


UNLOADING 


9190 
9200 
9210 
9220 
9230 
9240 
9250 


Lod1  = 
Lod2= 
Lod5= 
Lod9= 
Cod1  = 
Cod2= 


■0    !  TEMP  LOAD 
■0    »  TEMP  LOAD 
■ . 1 »Pmax   ! 
■0    !  MAX  LOAD 
■0    !  TEMP  COD 
■0    !  TEMP  COD 


Cod4=0 


!  TEMP  COD 


(PI) 
(P2) 
(P5) 
(P9) 
(D1) 
(D2) 
(C4> 


DATA  POINTS  TAKEN  DURING 


USED  IN  AREA  CALC 

USED  FOR  TURN  AROUND  SIGNAL 
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9260 

9270 

9280 

9290 

9300 

9310 

9320 

9330 

9340 

9350 

93G0 

9370 

9380 

9390 

9*00 

9410 

9420 

9430 

9440 

9450 

9460 

9470 

9480 

9490 

9500 

ADING 

9510 

9520 

9530 

9540 

9550 

9560 

9570 

9580 

9590 

9600 

9610 

9620 

9630 

9640 

9650 

9660 

9670 

9680 

9690 

9700 

9710 

9720 

9730 

9740 

9750 

9760 

9770 

9780 

9790 

9300 

9810 

9820 

9830 


MAX  COD  "  <D9> 
COD  DIFFERENCE  "  (D3) 
TEMP  DEFL      "    <D4)   USED  IN  AREA  CALC 

TEMP  DEFL      "    <D5>    

JI  FROM  CURRENT  UNLOAD 
JI  FROM  PREVIOUS  UNLOAD  DTNSDRC  (J7> 
FLAG  <0=LOADING.  1 'UNLOADING) 
FLAG  (F1=999  ENDS  TEST) 


Cod9-0 
Codiff»0 

Def 11=0 

Defl2=0 

J6=0 

J7=0 

11=0 

F1  =0 

; 

f 

PRINT  "" 
i 

IF  Ft =999  THEN  GOSUB  End_test 
; 

GOSUB  Take_data 

Defl6=Defl(I) 

Lod7=Lod<I) 

Lod2=Lod(I) 

Defl2=DefKI) 

Cod2=Cod<I) 

I  TEST  CONDITIONS 
IF  (K-K0K20  THEN 
Codiff=Cod2-Cod1 
IF  Codiff<-3.0E-4 


OF  DATA  POINT  JUST  TAKEN 

10430   !  CALC  AREA  4  TAKE  MORE  DATA 

THEN  9540   !  IF  CODIFF  IS  <-)  CHECK  IF  LOADING  OR  UNLO 


t0430   ! 

THEN  9600 

; 

10340   ! 

t 

!  MUST 
!  TEST 
0290 


IF  11=0  THEN 

IF  Cod2>Cod9 

GOTO    10340 

IF   I1O0    THEN 

Lod9=Lod2 

Cod9=Cod2 

Defl9=Defl2 

11=1 

GOTO    10340 

IF  L<5   THEN    1 

Kl-K-L 

K3=K1+L 

PENUP 

; 

K2=K2+1 
f 

GOSUB   Slope 

Stiff2=Stiff 

GOSUB   Crack_length        ! 

AKK2)=Crcicl  ! 

Delta_a<K2)=Al(K2)-A1 

} 

I  NOTE:  A1  -  INITIAL  CRACK  LENGTH;  ASSIGNED  IN  Initial_a  ROUTINE, 


IF  LOADING  CALC  AREA  &   TAKE  MORE  DATA 
!  IF  TRUE  CHECK  IF  GOOD  UNLOADING 

DISREGARD  UPPER  27.  OF  UNLOADING  CURVE 

IF  UNLOADING  CHECK  IF  ON  UPPER  27.  OF  CURVE 

IF  LOADING  MAKE  CURRENT  LOAD  =  MAXIMUM  LOAD 

MAKE  CURRENT  COD  =  MAX  COD 

MAKE  CURRENT  DEFL  =  MAX  DEFL 

BE  UNLOADING  . 

IF  UPPER  27.  OF  UNLOADING 
!  BAD  UNLOADING,  SET  11=1 
!  SETS  LOWER  INDEX  FOR  SLOPE  CALCULATION. 
!  SETS  UPPER  INDEX  FOR  SLOPE  CALCULATION. 


!  INCREMENT  UNLOADING  COUNTER 
!  IF  GOOD  UNLOADING 


CALC  CRACK  LENGTH 
CURRENT  CRACK  LENGTH 

!  DETERMINE  CRACK  GROWTH 


Crckl4=Crclcl 

PRINT   " 

!    PRINT    "" 

PRINT  "CRACK  LENGTH  -";A1<K2>  ;'W 

PRINT  "DELTA_a  =" ;Delta_a<K2) ; "mm" 
PRINT  "TOTAL  *   DATA  POINTS  -  ":K 
PRINT  "#  DATA  POINTS  THIS  UNLOADING 
PRINT  "CORRELATION  =";Correl 
PRINT  "COD  COMPLIANCE  =" ;Compl :"mm/N" 


UNLOADING  *";K2;"  — " 


'";L 
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9840 

3850 

9860 

9870 

9880 

9890 

9900 

9910 

9920 

9930 

9940 

9950 

9960 

9970 

9980 

9990 

10000 

10010 

10020 

10030 

10040 

10050 

10060 

10070 

10080 

10090 

10100 

10110 

10120 


Area3<K2>=Area 
• 

!  CALCULATE  957.  CONFIDENCE  INTERVAL 

Stiff=Stiff+T9*U4 

GOSUB   Craclc_length 

Amax*=Crck 1 

Stiff=Stiff-2»T9»U4 

GOSUB  Crack_length 

Stxff=Stiff2 

Amin=Crcl:  1 

PRINT   "" 

PRINT  "357.   CONFIDENCE  INTERVAL 
i 

GOSUB  J_calc   !  CALC  J  (USES  J-ERNST 
t 

lidefl<K2>-Defl<K>     !  MAX  DEFLECTION 
» 

!  PRINT  "" 

PRINT  "AREA  -"  ;Area9(K2)  ;  "N»mro" 

"J  =  ";J9<K2):"KJ/nr2" 

"LOAD  LINE  DEFLECTION  -"  ;Mdef  KK2) 


ON  CRACK  LENGTH  FROM  " ;Amin ;"TO";Amax 
MODEL) 


PRINT 
PRINT 
PRINT 

PRINT 

? 

t 

GOSUB 
j 

I  DATA 


mm 


"CHANGE  IN 
"K0-";K0;" 


COD 


=";Codiff 
";K1 


mm 


Prt_results  ! 

POINTS- CAN  BE  PRESENTED 


PRINT 
» 


PRINT  "  PT  #    LOAD(N) 

FOR  I=<K1-1>  TO  K 
!  PRINT  USING  91 70 : I.Lod< I) ,Def 1< I) ,Cod< I) 
!  IMAGE  4D,4X,5D.D,6X,DD.6D.7X,DD.6D 

NEXT  I 

WAIT  5        !  PAUSE  TO  ALLOW  VIEWING 


STORES  INFORMATION  ON  DISK 


HARACTER  (!) 
10130 
10140 
10150 
10160 
10170 
10180 
10190 
10200 
10210 
10220 
10230 
10240 
10250 
10260 
10270 
10280 
10290 
10300 
10310 
K 

10320 
10330 
10340 
RE  DATA 
10350  L=L+1 
10360   IF  Lod2>Lod9 
10370  BEEP        ! 
10380  BEEP        ! 
10390  BEEP        ! 
10400  GOTO  10440   ! 
10410   ! 
10420   !   CALCULATE 


PRINTS  RESULTS  ON  PRINTER 

FOR  REVIEW  BY  SIMPLY  REMOVING  THE  REMARK  C 

DEFLECTION(mm)    COD(mm)" 


GOSUB  Stor_data 
I 

; 

IF  Crclcl4-A1>0   THEN    10290 

PRINT   "" 

PRINT  "    NO  CRACK  EXTENSION  THIS  POINT." 

PRINT  "" 

11=0 

L=1 

!  IF  L-1  THEN  K1=K   !  SETS  LOWER  INDEX  FOR  SLOPE  CALCULATION  AT  CURRENT 


GOTO  10360 

IF  Cod2<Cod1  THEN  10440   ?   GET  MORE  DATA 

IF  Lod2>.98»Lod9  THEN  10440    !  DISREGARD 


UPPER  27.  OF  UNLOADING  &  GET  MO 


!  INCREMENT  DATA  POINT  COUNTER 

Lod5  THEN  10440      !  IF  TRUE  GET  MORE  DATA 

SENDS  OUT  BEEP 

AS  TURN  AROUND 

SIGNAL  TO  BEGIN  RELOADING 

RESUME  TAKING  AND  TEST  DATA 

AREA  UNDER  LOAD  vs  DEFLECTION  CURVE 
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10430 
10440 
10450 
10460 
10470 
10480 
10490 
10500 
10510 
10520 
10530 
10540 
10550 
10560 
10570 
10580 
10590 
10600 
10610 
10620 
10630 
10640 
10650 
10660 
10670 
10680 
10690 
10700 
10710 
10720 
10730 
10740 
10750 
10760 
10770 
10780 
10790 
10800 
10810 
10820 
10830 
10840 
10850 
10860 
10870 
10880 
10890 
10900 
10910 
10920 
10930 
10940 
10950 
10960 
10970 
10980 
10990 
11000 
11010 
11020 


Area=Area+(Def 12-Def 11 >»<Lod2+Lod1 >/2 

Lod1=Lod2 

K-K+1       !  INCREMENT  TOTAL  •   DATA  POINTS  COUNTER 

Cod1=Cod2 

Defl1=DefI2 

t 

GOTO  9380    !  LOOP  BACK  CHECK  IF  TEST  IS 

OVER,  IF  NOT  TAKE  MORE  DATA. 


RETURN     !  NEVER  REACH  THIS  BUT  ROUTINE  WOULD  NOT  BE  VALID  WITHOUT  IT 

jwiMHW WW *tHHW»M»«W««« 

Prt_resnlts:      !  PRINTS  TEST  INFORMATION  WHILE  PROGRAM  IS  RUNNING 
PRINTER  IS  706   !  TOGGLES  PRINTER 
PRINT  M" 

PRINT  "      UNLOADING  *";K2;" " 

PRINT  "•' 

PRINT  "CRACK  LENGTH  =":AL<K2>  :"mm    DELTA  _a  =  "  ;Delta_a<K2>  ;"mm" 

PRINT  "95%  CONFIDENCE  INTERVAL  ON  CRACK  LENGTH  IS  FROM" ; Am  in;"  TO";Amax 

PRINT  "M 

PRINT  "TOTAL  *  DATA  POINTS  =";K;"   *   DATA  POINTS  THIS  UNLOADING  =";L 

PRINT  "CORRELATION  =";Correl 

PRINT  "COD  COMPLIANCE  =" ;Compl ;"mm/N" 

PRINT  "" 

PRINT  "AREA  -" : Area9(K2) ;"N«mra" 

PRINT  "J  =  ";J9<K2>;"KJ/nT2" 

PRINT  "MAX  LOAD  LINE  DEFLECTION  =";Mdef 1 <K2) ;"mm" 

PRINT  "CHANGE  IN  COD  -" ; Cod if f ;"mm" 

f 

PRINT  "" 

PRINTER  IS  1    !  TOGGLES  PRINTER  TO  CRT 
f 

RETURN         !  RETURNS  TO  Main_prog 

**•*«*••  WHHHHHHM4HHHWMMWMHHNHHM  »•»■»■»»•*•»» 

Stor_arrays:   !  STORES  ARRAYS  (Lod.Def I ,Cod)  ON  DISK 

THE  ARRRAYS  CAN  BE  RECALLED  LATER  USING 
THE  "J_DATA"  PROGRAM. 

OUTPUT  9File;Lod<»),END 

OUTPUT  ^File;Defl(*),END 

OUTPUT  3File;Cod<»),END 
» 

BEEP 

f 

♦ 

RETURN 


!  RETURNS  TO  End  test 


End_test : 
t 

PRINT  CHRS<12) 

\ 

\ 

GOSUB   Stor_arrays 

f 


!    ENDS   TEST 


84 


11030 

11040 

11050 

11060 

11070 

11080 

11090 

11100 

11110 

11120 

11130 

11140 

11150 

11160 

11170 

11180 

11190 

11200 

11210 

11220 

11230 

11240 

11250 

IT 

11260 

11270 

11280 

11290 

11300 

11310 

11320 

11330 

11340 

11350 

11360 

11370 

11380 

11390 

11400 

11410 

11420 

11430 

11440 

11450 

11460 

11470 

11480 

11490 

11500 

11510 

11520 

11530 

11540 

11550 

11560 

11570 

11580 

11590 

11600 

11610 


MM 
I«M 


PRINT  "DATA  HAS  BEEN  STORED  ON  DISK  IN  FILE  ";Name5 
I 

GOSUB  R_curve 
i 

PRINT 
PRINT 
PRINT 
PRINT 
PRINT  '•" 
PRINT  "" 
PRINT  "" 
PRINT  " 
PRINT  "" 
PRINT  — 
PRINT  "M 
PRINT  "" 
PRINT  "" 
PRINT  " 
I 

t 

GOTO  12590 
; 

RETURN   !  NEVER  REACH  THIS  POINT  BUT  ROUTINE  WOULD  NOT  BE  VALID  WITHOUT 


TEST  IS  OVER!!' 


TEST  IS  OVER!' 


TEST  IS  OVER!" 


!  GO  TO  END  STATEMENT 


R_curve:   !  J  vs  DELTA_a  DATA  IS  PRESENTED  FOR  PLOTTING 
i 

PRINT  "" 
t 

!  ALL  (»)  STATEMENTS  WILL  ALLOW  DIRECT  INPUT  OF  Delta_a<I)  vs  J9<I>  DATA 

!  YOU  MUST  REMOVE  [!<»>]  FROM  BEGINNING  OF  EACH  LINE  TO  ACTIVATE 

!(-)  PRINT  "INPUT  THE  *  OF  UNLOADINGS" 

!<»>  INPUT  K2 

!(»)  PRINT  "ENTER  DELTA_A< I) , J9< I)" 

PRINT  "UNLOADING  #   DELTA_a(mm)     J(lcJ/m"2)" 

FOR  1=1  TO  K2 

!(»)  INPUT  Delta_a(I>,J9<I> 

PRINT  USING  1 1420:1, Delta_a< I) , J9< I) 

IMAGE  6X,3D,4X,DD.6D,5X,6D.D 
NEXT  I 
PRINT  "" 

PRINT  "DO  YOU  WANT  A  HARD  COPY?" 
INPUT  HS 

IF  H$="N"  OR  HS="NO"  THEN  11680 
; 

PRINTER  IS  706        !  TOGGLES  TO  THERMAL  PRINTER 
PRINT  "UNLOADING  *  DELTA_a<mm)    J<kJ/nT2)" 

FOR  1-1  TO  K2 

PRINT  USING  11 420; I.Del ta_a< I), J9(I> 
NEXT  I 
PRINTER  IS  1    !  RETURNS  TO  CRT  AS  PRINTER 

*)  PRINT  "ARE  ALL  OKAY?" 

»)  INPUT  GS 

•)  IF  G$-"Y"  THEN  11582 

*)  PRINT  "INPUT  I  (INDEX)  OF  BAD  #" 

»>  INPUT  LI 

*)  PRINT  "INPUT  VALUES  OF  DELTA  A<I>  AND  J9(D" 
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11620 

11630 

11640 

11650 

11660 

A 

11670 

11630 

11690 

11700 

11710 

11720 

11730 

'  OR  ' 

11740 

11750 

11760 

11770 

11780 

11790 

11800 

11810 

11820 

11830 

11840 

11850 

11860 

11870 

11880 

11890 

11900 

11910 

11920 

11930 

11940 

11950 

11960 

11970 

11980 

11990 

12000 

12010 

12020 

12030 

12040 

12050 

12060 

12070 

12080 

12090 

12100 

12110 

12120 

12130 

12140 

12150 

12160 

12170 

12180 

12190 


!(*>  INPUT  Delta_a(LD,J9<Ll> 

!(-)  PRINT  "ARE  YOU  OK" 

!<»)  INPUT  EeeS 

!(»)  IF  Eee$="N"  THEN  GOTO  11574 

!  END  OF  <*)  CHANGES  REQUIRED  FOR  DIRECT  INPUT  OF  DELTA_A(I>  vs  J9<I)  DAT 

t 

PRINT  "" 

PRINT  "DO  YOU  WANT  A  PLOT  OF  J  vs  DELTA_a?" 
INPUT  RS 

IF  R$="N"  OR  R$="NO"  THEN  12550 
PRINT  "" 

PRINT  "WHEN  THE  PLOTTER  IS  READY  (NEW  GRAPH  PAPER  INSTALLED)  PRESS  'ENTER 
CONTINUE'  KEY" 
INPUT  CS 
» 

GOSUB  Setup_plotter 

f  PLOT  THE  DATA  THE  DATA  POINTS 
VIEWPORT  13,100*RATIO,10,90 
WINDOW  XI (C2) .X2CC2) , Y1 <C2) , Y2(C2) 
FOR  1  =  1  TO  K2 

X=Delta_a<I> 

Y=J9(I) 

MOVE  X-X3».005,Y+Y3«.005 

DRAW  X+X3*.005 

MOVE  X+X3-.005.Y+Y3 

DRAW  X-X3».005,Y-Y3 
NEXT  I 

!  PLOT  THE  BLUNTING  LINE 
VIEWPORT  13,100»RATI0,10,90 
WINDOW  X1 <C2> ,X2<C2) , Y1 <C2) , Y2(C2) 
Bly=Y2(C2)-Y3/8 
Blx=Bly/(2»Fs)»1000 
MOVE  0,0 
DRAW  Blx.Bly 
t 

Blminx=Blx+.15 
MOVE  .15,0 
LINE  TYPE  4 
DRAW  Blminx.Bly 
t 

Blmaxx=Blx+! .5 
MOVE  1  .5,0 
LINE  TYPE  4 
DRAW  Blmaxx.Bly 
PENUP 

K1  =  1 


■Y3».005 
,005 
005 


DRAWS  A 
SMALL  X 
FOR  DATA 
POINT 


!  TAKES  LINE  TO  7/8  OF  Ymax 
!  DELTA  a  -  J/2»Fs 


!  DRAWS  LOWER  BOUNDS 


!  DRAWS  UPPER  BOUNDS 


K3=K2 
GOSUB  Slope 


!  LOWER  INDEX  FOR  LEAST  SQUARES  FIT 
!  UPPER  INDEX 


!  PLOT  THE  LINE  OBTAINED 
VIEWPORT  13.100»RATI0,10,90 
WINDOW  XKC2)  ,X2(C2)  ,YKC2>  ,Y2<C2) 

Xplt=X2<C2)-X3/8    !  TAKES  LINE  TO  7/8  OF  Xmax 

Yplt=Xplt*M+B 

MOVE  O.B 

LINE  TYPE  6 
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DRAW  Xplt.Yplt 

PENUP 
PRINT  "" 

PRINT  "ONLY  DATA  UITHIN  BOUNDS  IS  USED" 
PRINT  "FOR  LEAST  SQUARES  FIT." 
PRINT  "" 

PRINT  "SLOPE-" ;M;"  INTERCEPT-" ;B 
PRINT  "CORREL-":Correl 
PRINT  "NUMBER  OF  DATA  POINTS  USED:";K4 
♦ 

PRINTER  IS  706        !  TOGGLES  PRINTER 

PRINT  "" 

PRINT  "" 

PRINT  "DATA  FROM  J  vs  DELTA_a  PLOT" 

PRINT  "" 

PRINT  "SLOPE  =";M;"   Y- INTERCEPT  =";B 

PRINT  "CORRELATION  =";Correl;"  FR0M";K4: 

PRINTER  IS  1  !  TOGGLES  PRINTER 


12200 
12210 
12220 
12230 
12240 
12250 
12260 
12270 
12280 
12290 
12300 
12310 
12320 
12330 
12340 
12350 
12360 
12370 
12380 
12390 
12400 
12410 
12420 
12430 
12440 
12450 
12460 
12470 
12480 
12490 
12500 
12510 
12520 
12530 
12540 
12550 
12560 
12570 
12580 
12590 
12600 
12610 
12620 
12630  SUB  Stop_test 


DATA  POINTS. 


FIND  THE  INTERSECTION,  HENCE  J<Tc). 
SOLVE  BOTH  EQUATIONS  SIMULTANEOUSLY 
J  -  M*(DELTA_a)+B  =  <DELTA_a)*(2»Fs>  HENCE: 

J1c-<B/M>/(1/M-<1000/(2»Fs))> 

PRINT  "" 

PRINT  "     J1c  =";J1c 

PRINTER  IS  706         !  TOGGLES  PRINTER 

PRINT  "" 

PRINT  "      Jlc  «";J1c 

PRINT  "" 

PRINT  "" 

PRINT  "  TEST  IS  OVER!" 

PRINTER  IS  1  !  TOGGLES  PRINTER 

GOTO  12100      !  GO  TO  END  STATEMENT 
RETURN       !  RETURNS  TO  End  test  ROUTINE 


END 


!  END  OF  TEST  AND  PROGRAM 


MHHHHHHHHHHHNHIMHM  •»•»  ■»»■»■»  ■»■»»■»■»  ■*-»•»■»»•»» -IHHHH* 


!  STOPS  TEST 


12640 

12650 

ROGRAM 

12660 

12670 

12680 

12690 


k9-KEY  ACTS  AS  PROGRAM  INTERRUPT 
!  Ic9-KEY  IS  NOT  LABELED,  BUT  IS  ASSIGNED  A  HIGHER  PRIORITY  TI  INTERRUPT  P 


t 

COM  /J_intgl/  F1 

F1-999 

i 


!  SETS  STOP  FLAG 


12700  SUBEND 
12710   !«h»« 


►  -»♦•»-*»-»■»-»-»-»-•-»■»•»-»-«■■»■»■» 
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APPENDIX   F 


J   DATA   LISTING 


J.DATA   PROGRAM 
THIS   PROGRAM   WILL   ALLOW    YOU    TO   READ   A    DATA 
FILE  CREATED   BY   "J   INTGL".      YOU   MUST  KNOU 
HOW   MANY   UNLOAD INGS   WERE   RECORDED   DURING 
THE   TEST   TO  USE   THIS  PROGRAM. 

PROGRAM  WAS  WRITTEN   BY: 
LCDR   U.K.    TRITCHLER 
VERSION    t.O  DATE    lt/83 


10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

HO 

120 

130  DIM  Data$<1 :4)[t01 

140  DIM  Lod<8000).Def 1(8000), Cod<3000) 

150  DIM  AK100).Area9<?00>,Delta_a<100> 

160  DIM  Mdefl<100),J9(100) 

170  ! 

180  PRINT  CHRS(!2>     f  CLEARS  CRT 

190  ! 

200  DUMP  DEVICE  IS  706 

210  PRINT  "WHAT  DATA  FILE  DO  YOU  WISH  TO  READ?" 

220  INPUT  NameS 

230  ASSIGN  jFiIe  TO  NameS 

240  FOR  1-1  TO  a 

250  ENTER  SFile:DataS(I)C 1:101 

260  NEXT  I 

270  PRINT  "*• 

280  PRINT  "J-INTEGRAL  TEST  OF  ";Data5<4> ;"BY :  ";DataS<3> 

290  PRINT  "" 

300  PRINT  "DATE:  ";DataS(1>:"     TIME:  ";DataS<2> 

310  ENTER  6File;Crdcl,A«un.Amax.Co(npl,Correl,K 

320  PRINT 

330  PRINT 

340  PRINT  ••" 

350  PRINT  "aO«";Crclcl,"WITH  95Z  CONDIFENCE  INTERVAL  FROM"  ;Arain;"    TO":Ar»ax 

360  PRINT  "" 

370  PRINT  "FROM  COMPLIANCE  -  ";Compl                      ., 

380  PRINT  "" 

390  "PRINT  "CORRELATION-";Correl;"FROM";K;"DATA  POINTS" 

400  PRINT   "" 

410  PRINT  "DO  YOU  WANT  A  HARD  COPY  OF  THIS  INFORMATION?" 

420  INPUT  G$ 

430  IF  G$-"N"  OR  G$-"NO"  THEN  450 

440  DUMP  ALPHA 

450  PRINT  "" 

460  PRINT  "INPUT  THE  NUMBER  OF  UNLOADINGS  THAT  WERE  RECORDED  DURING  THE  TEST." 

470  INPUT  CI           !  WILL  BE  USED  AS  COUNTER 

480  FOR  1-1  TO  CT 

490  ENTER   $File;K2,AKK2>  ,Del  ta_a<K2)  .K.L.Correl  .Compl  ,Anun,Amax.Area9<K2>  . 

J9<K2> 

500  ENTER  9fx le:Mdef 1(K2> .Codlf f.KO .Kl 

510  GOSUB   Prt_resuits             !    DISPLAYS    INFO 

520  PRINT  "" 

530  PRINT  "DO  YOU  WANT  A  HARD  COPY?" 

540  INPUT  Q$ 

550  IF  QS-"N"  OR  Q$-"NO"  THEN  580 

560  PRINTER  IS  706 

570  GOSUB  Prt  results    !  WILL  NOW  PRINT  INFO 

580  NEXT  I 

590  PRINT  "DO  YOU  WANT  A  LISTING  OF  THE  DATA  POINTS?" 
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•f  M 


600 

610 

620 

630 

640 

650 

660 

670 

680 

690 

700 

710 

720 

730 

740 

750 

760 

770 

730 

790 

300 

810 

820 

830 

840 

850 

860 

DISKE 

870 

880 

890 

900 

910 

920 

930 

940 

950 

960 

970 

980 

990 

1000 

1010 

1020 

1030 

1040 

1050 

1060 

1070 

1030 

1090 

1100 

1110 

1120 

1130 

1140 

1150 

1160 

1170 

1180 


COD  (mm)" 


INPUT  ES 

IF  ES-"N"  OR  ES-"NO"  THEN  830 

PRINT  "" 

PRINT  "PT   *    LOAD  (N>     DEFLECTION  (mm) 

ENTER  *F 1 1 e : Lod( » > , Defl < • > , Cod ( » > 

FOR  1-1  TO  K 

PRINT  USING  670:1. Lod( I) .Def 1 ( I> ,Cod( I) 

IMAGE  X , 40 . 3X , 50 . D , 8X , 0 . 6D . 3X , D . 6D 
NEXT  I 
PRINT  "" 

PRINT  "DO  YOU  WANT  A  HARD  COPY?" 
INPUT  HS 

IF  HS-"N"  OR  HS-"NO"  THEN  830 
PRINTER  IS  706 

PRINT  "DATA  FROM  aO  FOR: ";DataS(4) ;"BY:  ";Data5<3> 
PRINT  "" 

PRINT  "TESTED  ON:  ":DataS( 1 ) ;"TIME:  ";Data$(2) 
PRINT  "" 

PRINT  "PT   #    LOAD  (N)     DEFLECTION  (mm)    COD  (mm)" 
FOR  1-1  TO  < 

PRINT  USING  670:1. Lod(I), Def  HI). Cod(I) 
NEXT  I 
PRINTER  IS  1 

PRINT  CHRSU2)  !  CLEARS  CRT 

PRINT  "DO  YOU  WISH  TO  PLOT  ANY  DATA?" 
INPUT  PS 

IF  ES-"Y"  OR  ES-"YES"  THEN  880   !  OTHERWISE  MUST  NOW  READ  DATA  POINTS  OFF 
TTE 
ENTER  *F 1 1 e : Lod ( • > . Def 1 ( • > , Cod < » ) 
IF  PS-"Y"  OR  PS-"YES"  THEN  GOSUB  Setap_plotter 
PRINT  "DO  YOU  WANT  ANOTHER  PL0T7" 
INPUT  AS 

IF  AS-"Y"  OR  AS-"YES"  THEN  GOSUB  Setup_plotter 
PRINT  "DO  YOU  WANT  ANOTHER  PLOT?" 
INPUT  FS 

IF  FS-"Y"  OR  FS-"YES"  THEN  GOSUB  Setup_plotter 
PRINT  "    PROGRAM  IS  OVER!!!" 
PRINT  "" 
PRINT  "" 
PRINT  "" 
PRINT 
PRINT 
PRINT  " 
PRINT  "" 
PRINT  "" 
PRINT  "" 
PRINT  **" 
PRINT  "" 

PRINT  " 
» 

GOTO  1290 
j 

» 

Prt_resalts: 
PRINT  "" 

PRINT  "FROM  UNLOADING  *" :K2: "RECORDED" ;L: "DATA  POINTS" 
PRINT  "COD  COMPLIANCE  -" :Compl :"mm/N   CORRELATION  -";Correl 
PRINT  "CRACK  LENGTH  -" ;AI (K2) :"mra   DELTA_A  -" :Del ta_a(K2) :"mm" 
PRINT  "95Z  CONFIDENCE  INTERVAL  FROM" ;Amin: "mm  TO" ; Amax ; "mm" 


»•*• 

!••* 


PROGRAM  IS  OVER! 


PROGRAM  IS  OVER!" 


!  DISPLAYS  OR  PRINTS  UNLOADING  INFORMATION 
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;!90  PRINT  "AREA  -"  ;Area:3<;<2>  :"N-nm   J  ■"  :  J9<<2>  ;  "KPa*m" 

1200  PRINT  "DATA  POINT  INDEXES  FROM" :K1 ;"TO" ;K :"  WERE  USED 

1210  PRINT  ,,M 

1220  PRINTER  IS  1 

1230  J 

1240  RETURN 

1250  J 

1 260  !  ••••••••••••••••— ••••«••••*••«  •••••«•••— 

1270  ! 

1280  Setup _plotter:      !  SET  UP  PLOTTER 

1290  !  YOU  CAN  CHOOSE  BETWEEN  THREE  DIFFERENT  PLOTS 

1300  PRINT  CHRS(12)      »  CLEARS  CRT 

1310  DIM  TSC80J.XSC801.YSC80] 

1320  GCLEAR 

1330  GINIT 

1340  GRAPHICS  ON 

1350  PLOTTER  IS  3. "INTERNAL"    !  PLOTS  ON  CRT 

1360  PRINT  "TYPE  OF  PLOT:" 

1370  PRINT  "    1.   LOAD  vs  COD" 

1380  PRINT  "    2.   LOAD  vs  DEFLECTION" 

1390  PRINT  "    3.   J  vs  CRACK  EXTENSION" 

1400  PRINT  "ENTER  THE  e   OF  THE  PLOT  YOU  DESIRE" 

1410  INPUT  C2 

1420  PRINT  "" 

1430  PRINT  "" 

1440  IF  C2-1  THEN  1470 

1450  IF  C2-2  THEN  1520 

1460  IF  C2-3  THEN  1570 

1A70  TS«"PLOT  of  LOAD  vs  COD  for  "4DataS<4> 

1480  PRINT  TS 

1490  XS-"COD  (mm)" 

1500  Y$-"LOAD  <N)" 

1510  GOTO  1610 

1520  TS-"PLOT  of  LOAD  vs  DEFLECTION  for  "4DataS(4> 

1530  PRINT  TS 

1540  XS-"DEFLECTION  (mm)" 

1550  YS-"L0AD  <N)"  ■„. 

1560  GOTO  1610 

1570  TS-"PLOT  OF  J  vs  CRACK  EXTENSION  for  ,,4DataS<4> 

1580  PRINT  TS 

1590  X$-"DELTA_a  <mm>" 

1600  YS-"J   <kJ/m-2>" 

1610  PRINT  "" 

1620  PRINT  "USING  METRIC  UNITS  (mm.N.k J/m*2.etc.  > :" 

1630  PRINT  "ENTER  THE   Xmin.Xmax , Ymin, Ymax" 

1640  INPUT  X1(C2),X2<C2),YKC2),Y2(C2) 

1650  X3-X2<C2)-X1(C2) 

1660  Y3-Y2<C2)-YKC2) 

1670  VIEWPORT  13. 100»RATI0, 10 ,90 

1680  WINDOW  XKC2>,X2(C2),Y1(C2),Y2<C2) 

1690  PRINT  "" 

1700  PRINT  "INPUT  DELTA  X  AND  DELTA_Y  LABEL  INTERVAL" 

1710  PRINT  "FOR  THE  TWO~AXES:" 

1720  INPUT  X4.Y4 

1730  PRINT  "M 

1740  PRINT  "DO  YOU  WISH  TO  HAVE  A  GRID  OVERLAY?" 

1750  INPUT  GS 

1760  IF  G$-"N"  OR  G$«"NO"  THEN  1800 

1770  PRINT  '"• 

1780  PRINT  "INPUT  GRID  LINES  PER  DELTA  LABEL:" 
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t 790  INPUT  BO 

1800  PRINT  "PRESS  'ENTER'  OR  'CONTINUE'  KEY  TO  PLOT" 

1810  INPUT  HS  !  DUMMY  VARIABLE 

1820  PRINT  CHRSU2)    !  CLEARS  CRT 

1830  IF  GS-"N"  OR  G3-"NO"  THEN  1970 

1340  FOR  I-X4/B0+XKC2)  TO  X2<C2)  STEP  X4/B0 

1850     MOVE  I.YKC2) 

1860     DRAW  I.Y2<C2> 

1870  NEXT  I 

1880  FOR    I-Y4/B0  +  YKC2)    TO   Y2<C2)    STEP   Y4/B0 

1890  MOVE   XKC2KI 

1900     DRAW  X2(C2),I 

1910  NEXT  I 

1920  GOTO  1970 

1930  IF  GS-"YM  OR  GS-"YES"  THEN  1970 

1940  MOVE  X2(C2).Y1(C2) 

1950  DRAW  X2(C2),Y2(C2) 

1960  DRAW  XKC2),Y2(C2> 

1970  CLIP  ON 

1980  AXES  X4,Y4,XKC2>,YKC2> 

1990  CLIP   OFF 

2000  VIEWPORT  13. 100-RATI0.5.90 

2010  WINDOW  X1(C2),X2(C2),Y1(C2),Y2(C2) 

2020  LORG  4 

2030  CSIZE  3,. 6 

2040  FOR  X-XKC2)  TO  X2<C2)  STEP  X4 

2050     CLIP  OFF 

2060     MOVE  X.YKC2) 

2070     LABEL  USING  "K";X       !  PLOTS  X-DIGITS 

2080  NEXT  X 

2090  VIEWPORT  4. 100-RATIO. 10 ,90 

2100  WINDOW  X1(C2),X2<C2)fYKC2),Y2<C2) 

2T10  LORG  2 

2120  FOR  Y-YKC2)  TO  Y2CC2)  STEP  Y4 

2130     CLIP  OFF 

2T40     MOVE  X1(C2),Y 

2150     LABEL  USING  "K";Y 

2160  NEXT  Y 

2170  VIEWPORT  0, 100-RATIO.O . 100 

2180  WINDOW  X1^2),X2<C2),YKC2),Y2<C2) 

2190  CSIZE  5 

2200  LORG  6         !  PLOTS  TITLE 

2210      MOVE  X2(C2)/2,Y2<C2) 

2220      LABEL  T$ 

2230  DEG 

2240  LDIR  90 

2250  CSIZE  4f.7 

2260     MOVE  XKC2),Y2(C2>/2 

2270     LABEL  Y$      !  PLOTS  Y-AXIS  LABEL 

2280  LORG  4 

2290  LDIR  0 

2300     MOVE  X2<C2)/2,YKC2) 

2310     LABEL  X$      !  PLOTS  X-AXIS  LABEL 

2320  ! 

2330  PRINT  "" 

2340  PRINT  "- 

2350  PRINT  "" 

2360  PRINT  •"• 

2370  PRINT  "       IS  THE  PLOT  (AXES  AND  LABELS)  OK?' 

2380  PRINT  "" 
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2390  PRINT  "  NOTE:   IF  YOU  ANSWER  YES  OR  Y" 

2400  PRINT  "  PLOT  WILL  BE  SENT  TO" 

2410  PRINT  "  EXTERNAL  PLOTTER" 

2420  INPUT  FS 

2^30  IF  FS-"NM  OR  FS-"NO"   THEN   GOSUB   Setup_plotter      f    REDO   PLOTTER   SETUP 

2440  ! 

2450  GCLEAR      !  CLEARS  PLOT  FROM  CRT 

2460  ! 

2470  G0SU8  Ext_plotter 

2480  ! 

2490  PENUP 

2500  MOVE  0,0 

25  tO  ! 

2520  GOSUB  Plot.data 

2530  ! 

2540  RETURN 

2550  ! 

2560  !  »—•••«.•*••*•—*•••——•••••••—•••••— 

2570  ! 

2580  Ext_plotter:    !  DUMPS  PLOTTER  SETUP  ON  EXTERNAL  HP-7225B  PLOTTER 

2590  PLOTTER  IS  705,"HPGL" 

2600  VIEWPORT  !3. 100-RATIO. 10 .90 

2610  WINDOW  X1(C2).X2(C2).Y1<C2),Y2<C2) 

2520  IF  GS-"N"  OR  G3-"NO"  THEN  2730 

2630  FOR  I-X4/B0+XKC2)  TO  X2<C2)  STEP  X4/B0 

2640     MOVE  I.YKC2) 

2650     DRAW  I,Y2<C2> 

2660  NEXT  I 

2670  FOR  I-Y4/B0  +  YKC2)  TO  Y2(C2>  STEP  Y4/B0 

2680     MOVE  XKC2),I 

2690     DRAW  X2<C2),I 

2700  NEXT  I 

2710  GOTO  2760 

2720  IF  GS-"Y"  OR  GS-^YES"  THEN  2760 

2730  MOVE  X2<C2),YKC2) 

2740  DRAW  X2<C2> , Y2(C2) 

2750  DRAW  X1 <C2) , Y2<C2) 

2760  CLIP  ON 

2770  AXES   X4.Y4,XKC2>,YKC2> 

2780  CLIP  OFF 

2790  VIEWPORT  13, 100«RATT0,5,90 

2800  WINDOW  X1<C2),X2(C2),Y1(C2),Y2(C2) 

2810  LORG  4 

2820  CSIZE  3,. 6 

2830  FOR  X-XKC2)  TO  X2(C2)  STEP  X4 

2840     CLIP  OFF 

2850     MOVE  X.YKC2) 

2860     LABEL  USING  "K";X 

2870  NEXT  X 

2880  VIEWPORT  4, 100-RATIO, 10 .90 

2890  WINDOW  X1  (C2>  ,X2<C2)  ,  YHC2) ,  Y2CC2) 

2900  LORG  2 

2910  FOR  Y-YKC2)  TO  Y2(C2>  STEP  Y4 

2920     CLIP  OFF 

2930     MOVE  XKC2),Y 

2940     LABEL  USING  "K";Y 

2950  NEXT  Y 

2960  VIEWPORT  0  .  I0O-RATIO.0 . 100 

2970  WINDOW  X1 <C2> ,X2(C2) , Y1 <C2> , Y2<C2) 

2980  LORG  2 
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2990 

CSIZE  5 

3000 

LORG  6 

3010 

MOVE  X2(C2)/2.Y2<C2) 

3020 

LABEL  TS 

3030 

DEG 

3040 

LDIR  90 

3050 

CSIZE  4,. 7 

3060 

MOVE  XHC2).Y2(C2)/2 

3070 

LABEL  YS 

3080 

LORG  4 

3090 

LDIR  0 

3100 

MOVE  X2(C2>/2,Y1(C2) 

3110 

LABEL  XS 

3120 

» 

3130 

PENUP 

3  MO 

MOVE  0.0 

3150 

f 

3160 

RETURN 

3170 
3180 
3190 

3200 

Piot  data:       !  PLOTS  TO  DATA  SELECTED 

3210 

t 

3220 

VIEWPORT  13. 100»RATIO. 10.90 

3230 

WINDOW  XI (C2> ,X2(C2) ,Y1 (C2> , Y2(C2) 

3240 

t 

3250 

IF  C2-I  THEM  3290 

3260 

IF  C2-2  THEN  3340 

3270 

IF  C2-3  THEN  3390 

3280 

i 

3290 

FOR  1-1  TO  K 

3300 

PLOT  Cod(I),Lod<I) 

3310 

NEXT  I 

3320 

GOTO  3430 

3330 

t 

3340 

FOR  1-1  TO  K 

3350 

PLOT  DefKI).Lod(I) 

3360 

NEXT  I 

3370 

GOTO  3430 

3380 

i 

3390 

FOR  1-1  TO  CI 

3400 

PLOT  Delta  a<K2).J(K2) 

3410 

NEXT  I 

3420 

i 

3430 

PENUP 

3440 

BEEP 

3450 

RETURN       !  RETURNS  TO  SETUP_PLOTTER 

3460 

3470 

3480 

3490 

END 
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APPENDIX  G 


PRECRACK  OUTPUT 


USING  ENGLISH  UNITS  <in.lbs,psi)  THE  FOLLOWING  INFORMATION 
HAS  BEEN  ENTERED  FOR:  HY80-5B 

UIDTH  -  .9469        NOTCH  LENGTH  =  .444 
Bmax  •  .496       Bmin  -  .394 
HALF  SPAN  '    1.894   PQISSON'S  RATIO  »  .3  ■ 
YIELD  STRENGTH  -  86600  ULTIMATE  TS  -  100000 
ELASTIC  MODULUS  -  2.9E+7  FLOW  STRESS  -  93300 

CONVERTING  THIS  DATA  TO  METRIC  UNITS  (mra.N.KPa)  FOR:  HY80-5B 

WIDTH  -  24.05126       NOTCH  LENGTH  -  11.2776 

Bmax  -  12.5984     Bnun  -  10.0076 

YIELD  STRENGTH  »  597085.603738   ULTIMATE  TS  -  689475.293 

ELASTIC  MODULUS  -  I .9994783497E+8   FLOW  STRESS  -  643280.448369 

BASED  ON  CURRENT  INFORMATION  FOR  HY80-5B: 

CURRENT  CRACK  LENGTH  -  .444  in 

OPTIMUN  CRACK  LENGTH  -  .56314  in 

EXTENSION  REQUIRED  -  .12414  in 

CRACK  LENGTH  AT  507.  OF  EXTENSION  -  .50607  in 

Pmax  -  3272.42651647  Lbs  Pmm  -  163.621325823  Lbs 

Delp  -  3108.80519064  Lbs  Kmax  -  32563.7466931 

DELTA_K  -  78435.5593585  in*. 5  DELTA_KE  -  40000  in~ .5 

NEW  LOADS  BASED  ON  DELTA_KE  ARE: 

Pmax  -  1668.84843723  Lbs  Pnun  -  83.4424218615  Lbs 

DELTA_P  -  1585.40601537  Lbs   Kmax  -  42105.2631579 


CONVERSION  DATA 
FOR  LOAD:-1 779.288646  N/V   (-400  Lbs/V> 
FOR  COD:  .0508  mm/V   (  .002  in/V) 
FOR  DEFLECTION:-12.7  mm/V   (-.5  in/V> 
U  -  .07 
a/W  -  .732186117982  in  CRACK  LENGTH  -  .693307035117  in 


LOAD  (lbs; 

*  COD  (in) 

COD  (Volts) 

200 

.00256290169302 

1.28145084651 

300 

.00384435253953 

1.92217626976 

400 

.00512580338604 

2. 56290 1693ri2 

500 

.00640725423255 

3.20362711627 

600 

.00768870507906 

3.84435253953 

700 

.00897015592557 

4.48507796278 

800 

.0102516067721 

5.12580338604 

900 

.0115330576186 

5.76652880929 

1000 

.0128145084651 

6.40725423255 

U00 

.0140959593116 

7.0479796558 

1200 

.0153774101581 

7.68870507906 

1300 

.0166588610046 

8.32943050231 

1400 

.0179403118511 

8.97015592557 

1500 

.0192217626976 

9.61088134832 

1600 

.0205032135441 

10.2516067721 

1700 

.0217846643907 

10.8923321953 

1800 

.0230661152372 

tl.  5330576186 

U  -  .075 

a/W  -  .713968763392  in  CRACK 

LENGTH  -  .676057022056 

LOAD  (lbs) 

COD  (in) 

COD  (Volts) 

200 

.00220863052208 

1.10431526104 

300 

.00331294578312 

1  .65647289156 

400 

.00*41726104416 

2.20863052208 

500 

.00552157630521 

2.7607881526 

in 


f  5  J 
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USING   ENGLISH   UNITS    (in.ibs.psi)    THE   FOLLOWING    INFORMATION 
HAS   3EEN   ENTERED   FOR:    HY80-5B 

UIDTH    -    .9469  NOTCH  LENGTH   -    .444 

Bmax    ■    .496  Bmin    -    .394 

HALF   SPAN    =■    1.394     POISSON'S  RATIO    ■    .3 
YIELD   STRENGTH    -   86600      ULTIMATE    TS    -    100000 
ELASTIC  MODULUS   -   2.3E+7      FLOW  STRESS   -   93300 

CONVERTING    THIS   DATA    TO   METRIC  UNITS    <mm.N.KPa>    FOR:    HY80-5B 

UIDTH   -   24.05126  NOTCH  LENGTH    -    11.2776 

Bmax    ■    12.5984  Bmin    -    10.0076 

YIELD   STRENGTH   -   597085.603733        ULTIMATE    TS    -   689475.293 

ELASTIC  MODULUS   -    1 .9994783497E+8        FLOW   STRESS    »   643280.448369 

BASED  ON   CURRENT    INFORMATION   FOR   HY80-5B: 

CURRENT   CRACK   LENGTH   -    .444    in 

OPTIMUM  CRACK   LENGTH   =■    .56314    in 

EXTENSION   REQUIRED   -    .'12414    in 

CRACK   LENGTH   AT   507.   OF  EXTENSION   -    .50607    in 

Pmax    -   3272.42651647   Lbs     Pmin   -    163.621325823   Lbs 

Delp    -   3108.80519064  Lbs     Kmax    •   82563.7466931 

DELTA_K    »    78435.5593585    iiT.5      DELTA_KE   -    40000    in'. 5 

NEW  LOADS      BASED  ON   DELTA_KE   ARE: 

Pmax    -    1668.84843723  Lbs     Pmin   -   83.4424218615   Lbs 

DELTA_P   »    1585. 4C60 1537  Lbs        Kmax   -   42105.2631579 

CONVERSION  DATA 
FOR   LDAD:-1 779.288646  N/V        <-400   Lbs/V) 
FOR  COD:    .0508  mm/V        <    .002    m/V> 
FOR   DEFLECTION: -12. 7   mm/V        (-.5    in/V> 

THE  FOLLOWING    INFORMATION  HAS   BEEN  OBTAINED: 

CRACK   LENGTH   -    .477215195165    in   FROM   MEASUREMENT   *    1 

WITH   957.C0NFIDENCE    INTERVAL  FROM    .47792545433      TO    .476507410981    in 

USING   COD   COMPLIANCE   OF:    2.841 0431 1 054E-6    in/Lb   -    1 .629051 73472E-5   mm/N 

WITH  A  CORRELATION  OF:    .999771386472  FROM    106      DATA  POINTS. 

THE  FOLLOWING    INFORMATION  HAS  BEEN  OBTAINED: 

CRACK   LENGTH   -    .477198279382    m   FROM   MEASUREMENT   #   2 

WITH  95ZC0NFIDENCE   INTERVAL  FROM    .478064054454      TO    .476336178981    in 

USING   COD   COMPLIANCE  OF:    2.84076126094E-6    in/Lb   -    1 .63037750213E-5  mn/N 

WITH  A  CORRELATION  OF:    .99959257929  FROM    124      DATA  POINTS. 

BASED  ON   CURRENT   INFORMATION  FOR  HY80-5B: 

CURRENT   CRACK   LENGTH   •    .477198279382    in 

OPTIMUN   CRACK   LENGTH   -    .56314    in 

EXTENSION   REQUIRED   »    .0909417206182    in 

CRACK  LENGTH   AT  507.  OF  EXTENSION   »    .489470860309   in 

Pmax    -   2854.63728597   Lbs      Pmin    -    142.731864298   Lbs 

Deip    -   2711.30542167   Lbs      Kmax    •   80310.2753096 

DELTA_K    -    76294. 7G15441    in". 5      DELTA_KE    -   40000    m    .5 

NEW  LOADS      BASED   ON   DELTA_KE   ARE: 

Pmax    -    1496.63606161    Lbs      Pmin    -    74.8313030804    Lbs 

DELTA  P   -    1421.30425353   Lbs        Kmax    =   ^2105 .2631579 
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".-£  r.jL_3wulG  INFORMATION  -HS  SEEN  JB'^I.'itD: 

CRACK  LENGTH  -  .  4393151  ill  398  in  FROM  MEASUREMENT  *  1 

Wi'TH  357.C0NFIDENCE  INTERVAL  FROM  .500434044632   TO  .498202259284  t-n 

USING  COD  COMPLIANCE  OF:  3.23884557474E-6  in/Lb  -  1 .361 83031 428E-5  mm/N 

WITH  A  CORRELATION  OF:  .399477933275  FROM  34  DATA  POINTS. 

THE  FOLLOWING  INFORMATION  HAS  BEEN  OBTAINED: 

CRACK  LENGTH  -  .499773214641  in  FROM  MEASUREMENT  *   2 

WITH  357.CONFIDENCE  INTERVAL  FROM  .50177796077   TO  .497787999797  in 

USING  COD  COMPLIANCE  OF:  3 .24775972936E-6  in/Lb  -  1  .8769892409E-5  mn/N 

WITH  A  CORRELATION  QF:  .999019979378  FROM  57  OATA  POINTS. 

BASED  ON  CURRENT  INFORMATION  FOR  HY80-5B: 

CURRENT  CRACK  LENGTH  -  .499773214641  m 

OPTIMUN  CRACK  LENGTH  -  .56314  in 

EXTENSION  REQUIRED  -  .0633667853593  in 

CRACK  LENGTH  AT  507.  OF  EXTENSION  -  .47818339268  in 

Pmax  -  1034.73227996  Lbs  Pmin  -  51.736612998  Lbs 

Delp  -  982.995665363  Lbs  Km*x  -  31462.7520142 

DELTA.K  -  29889.5146985  in'. 5  DELTA_KE  -  29000  in". 5 

NEW  LOADS  BASED  ON  DELTA_KE  ARE: 

Pmaa  -  1003.9351936  Lbs  Pmin  -  50.1967596798  Lbs 

DELTA_P  -  953.738433917  Lbs   Kraax  -  30526.3157895 

THE  FOLLOWING  INFORMATION  HAS  BEEN  OBTAINED: 

CRACK  LENGTH  -  .523238325069  in  FROM  MEASUREMENT  *  1 

WITH  95ZC0NFIDENCE  INTERVAL  FROM  .524952826657   TO  .521538109085  in 

USING  COD  COMPLIANCE  OF:  3. 74604099481 E-6  in/Lb  -  2. 16181567648E-5  mm/N 

WITH  A  CORRELATION  OF:  .999t83456903  FROM  61   DATA  POINTS. 

THE  FOLLOWING  INFORMATION  HAS  BEEN  OBTAINED: 

CRACK  LENGTH  -  .519260644793  in  FROM  MEASUREMENT  #  2 
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WITH  95%C0NFIDENCE  INTERVAL  FROM  .521257634776  TO  .51728301025  in 
USING  COD  COMPLIANCE  OF:  3 .65545467572E-6  in/Lb  ■  2 . 11 309694446E-5  mm/N 

WITH  A  CORRELATION  OF:  .998905179767  FROM  61   DATA  POINTS. 

i 

THE  FOLLOWING  INFORMATION  HAS  BEEN  OBTAINED: 

CRACK  LENGTH  =  .520164834896  in  FROM  MEASUREMENT  #  3 

WITH  957.C0NFIDENCE  INTERVAL  FROM  .521732160051   TO  .51860945694  in 

USING  COD  COMPLIANCE  OF:  3.67581279234E-6  in/Lb  -  2 . 1 1927866299E-5  nWN 

WITH  A  CORRELATION  OF:  .999107091656  FROM  79  DATA  POINTS. 

3ASED  ON  CURRENT  INFORMATION  FOR  HY80-5B: 

CURRENT  CRACK  LENGTH  -  .520164834396  in 

OPTIMUN  CRACK  LENGTH  =  .56314  in 

EXTENSION  REQUIRED  -  .0479751651041  in 

CRACK  LENGTH  AT  507.  OF  EXTENSION  -  .467987582552  xn 

Pmax  -  942.504620395  Lbs  Pnun  -  47.1252310197  Lbs 

Delp  -  895.379389375  Lbs  Kmax  -  30831.6643108 

DELTA_K  -  29290.0816653  in". 5  DELTA_KE  -  29000  in". 5 

NEW  LOADS  BASED  ON  DELTA_KE  ARE: 

Pmax  -  933.170289649  Lbs  Pmin  -  46.6585144825  Lbs 

DELTA_P  -  886.511775167  Lbs   Kmax  -  30526.3157895 

THE  FOLLOWING  INFORMATION  HAS  BEEN  OBTAINED: 

CRACK  LENGTH  -  .557542180253  in  FROM  MEASUREMENT  #  1 

WITH  95/.C0NFIDENCE  INTERVAL  FROM  .559342565208  TO  .555757596637  in 

USING  COD  COMPLIANCE  OF:  4 .6522934609E-6  in/Lb  -  2 .6377228564E-5  mm/N 

WITH  A  CORRELATION  OF:  .999026349545  FROM  60  DATA  POINTS. 

THE  FOLLOWING  INFORMATION  HAS  BEEN  OBTAINED: 

CRACK  LENGTH  -  .556985314799  in  FROM  MEASUREMENT  #  2 

WITH  957.C0NFIDENCE  INTERVAL  FROM  .559603140017   TO  .554400760935  in 

USING  COD  COMPLIANCE  OF:  4.63555555373E-6  in/Lb  -  2. 6922753751 8E-5  mm/N 

WITH  A  CORRELATION  OF:  .998286279816  FROM  51   DATA  POINTS. 

THE  FOLLOWING  INFORMATION  HAS  BEEN  OBTAINED: 

CRACK  LENGTH  -  .55746232304  in  FROM  MEASUREMENT  *   3 

WITH  95XCONFIDENCE  INTERVAL  FROM  .558944733365  TO  .555990641885  in 

USING  COD  COMPLIANCE  OF:  4.64988856512E-6  in/Lb  -  2.68079074153E-5  mm/N 

WITH  A  CORRELATION  OF:  .999219985095  FROM  70  DATA  POINTS. 

BASED  ON  CURRENT  INFORMATION  FOR  HY80-5B: 

CURRENT  CRACK  LENGTH  -  .55746232304  in 

OPTIMUN  CRACK  LENGTH  -  .56814  in 

EXTENSION  REQUIRED  -  .0106776769602  in 

CRACK  LENGTH  AT  50%  OF  EXTENSION  -  .44933883348  in 

Pmax  -  784.95101444  Lbs  Pmin  -  39.247550722  Lbs 

Delp  -  745.703463718  Lbs  Kmax  -  29594.1425953 

DELTA_K  -  28114.4354655  in*. 5  DELTA_KE  -  29000  in". 5 

THE  FOLLOWING  INFORMATION  HAS  BEEN  OBTAINED: 

CRACK  LENGTH  -  .556873246958  in  FROM  MEASUREMENT  #  1 

WITH  35 /.CONFIDENCE  INTERVAL  FROM  .558630143876  TO  .555131398279  in 

USING  COD  COMPLIANCE  OF:  4.6321961 2139E-6  in/Lb  =■  2 .57532481 877E-5  mm/N 
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WfTH  A  CORRELATION  OF:  .399!  74 1  71419  FROM  54   DATA  POINTS. 

BASED  ON  CURRENT  INFORMATION  FOR  HY30-58: 

CURRENT  CRACK  LENGTH  -  .556873246958  in 

CPTIflUN  CRACK  LENGTH  -  .56314  in 

EXTENSION  REQUIRED  -  .0112667530417  in 

CRACK  LENGTH  AT  507.  OF  EXTENSION  -  .449633376521  in 

Pmax  -  787.327495262  Lbs  Pmin  -  39.3663747631  Lbs 

Delp  -  747.961120499  Lbs  Kmax  ■  29614.5202835 

DELTA_K  -  28133.3037694  in- .5  DELTA.KE  -  23000  in"  .5 

THE  FOLLOWING  INFORMATION  HAS  BEEN  OBTAINED: 

CRACK  LENGTH  »  .560000734225  in  FROM  MEASUREMENT  «  1 

UITH  95XC0NFIDENCE  INTERVAL  FROM  .56106384425   TO  .55893721925  in 

USING  COD  COMPLIANCE  OF:  4 . 7270940221 8E-6  in/Lb  -  2 . 71 807983847E-5  mm/N 

WITH  A  CORRELATION  OF:  .99933646638  FROM  101   DATA  POINTS. 

THE  FOLLOWING  INFORMATION  HAS  BEEN  OBTAINED: 

CRACK  LENGTH  -  .562208107534  m  FROM  MEASUREMENT  *  2 

WITH  957.CONFIDENCE  INTERVAL  FROM  .56243540602   TO  .560988130658  in 

USING  COD  COMPLIANCE  OF:  4 . 79552996551E-6  in/Lb  -  2 . 76035025025E-5  mm/N 

WITH  A  CORRELATION  OF:  .99948124671  FROM  67  DATA  POINTS. 

BASED  ON  CURRENT  INFORMATION  FOR  HY30-5B: 

CURRENT  CRACK  LENGTH  *  .562208107534  in 

OPTIMUN  CRACK  LENGTH  -  .56314  in 

EXTENSION  REQUIRED  -  .00593139246594  in 

CRACK  LENGTH  AT  507.  OF  EXTENSION  -  .446965946233  in 

Pmax  -  765.936366326  Lbs  Pmin  -  38.2963183163  Lbs 

Delp  -  727.63954301  Lbs  Kmax  -  29428.399829 

DELTA_K  -  27957.4548376  in  .5  DELTA_KE  ■  29000  in*  .5 

THE  FOLLOWING  INFORMATION  HAS  BEEN  OBTAINED: 

CRACK  LENGTH  =  .564151471472  in  FROM  MEASUREMENT  *  1 

WITH  35/.CONFIDENCE  INTERVAL  FROM  .565674835605  TO  .562639453926  in 

USING  COD  COMPLIANCE  OF:  4 .85680462792E-6  in/Lb  •  2.801 124373Z9E-5  mm/M 

WITH  A  CORRELATION  OF:  .999133273453  FROM  72  DATA  POINTS. 

THE  FOLLOWING  INFORMATION  HAS  BEEN  OBTAINED: 

CRACK  LENGTH  -  .568361917225  in  FROM  MEASUREMENT  #  2 

WITH  957.CONFIDENCE  INTERVAL  FROM  .569305393318   TO  .567422811142  in 

USING  COD  COMPLIANCE  OF:  4 .99294827909E-6  m/Lb  -  2 .86883916896E-5  mm/N 

WITH  A  CORRELATION  OF:  .999601)697254  FROM  84  DATA  POINTS. 

PRECRACKING  IS  COMPLElED! 

FINAL  CRACK  LENGTH  -  .563361917225 
CORRELATION  OF  DATA  -  .999601697254 

REMOVE  SPECIMEN  FROM  MTS. 
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APPENDIX  H 


J  INTGL  OUTPUT 


THREE    POINT   BEND   SPECIMEN 
TEST   REPORT 


SPECIMEN    I. P.    H  Y SO- SB     MATERIAL     /*//  SO 
Data     I Z  J 13 1  9,3 Operator    "7*  >  -rc^f  <-  g>? 


ORIENTATION   T~  S> 
Taat   Taap     2.)  C 


T 

Bmax 


Bnin 

i 


W    -       0.  9+(g  9~  A    -       Q.4-44-  b    «        Q.  303.9  "  L    ,      «4,  5*  " 


f.g<?4' 


Bmax    »     0'49G,0' 


Bain    -        g.^P4 
Ya    «       8b(e>CO   ■/*+*■  uta    •     /COOOC^QU^i  Fa    *        933C)<Z  jau^ 


Poiaaon'a    ratio    ■        Q*  *3 


Elaatic   Modulua    •      X-9  </Q    fL+Z 


Consents: 

NY  BP ft**r.     /Vfr.TXL         F«a-      P*r.u*.*-rvA h/fcfcfl . 

■P«ccc«^g*       q4        ta.  hits;         3*-  Tcs-r    .«,*j      ^//,i/fl3  ^  T,.'/^^ 
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UCLMG  ENGLISH  UNITS  Cin.Ibs.psi>  THE  FOLLOWING  INFORMATION 
HAS  BEEN  ENTERED  FOR:  HY30-5B 

WIDTH  »  .9469        NOTCH  LENGTH  -  .444 
Bmax  -  .496       Bmin  -  .394 
HALF  SPAM  -  1.394  POISSON'S  RATIO  -  .3 
YIELD  STRENGTH  -  86600  ULTIMATE  TS  -  100000 
ELASTIC  MODULUS  -  2.9E+7  FLOW  STRESS  -  93300 

CONVERTING  THIS  DATA  TO  METRIC  UNITS  (mm.N.KPa)  FOR:  HY80-5B 

WIDTH  -  24.05126       NOTCH  LENGTH  =  11.2776 

Bmax  ■  12.5984     Bnun  -  10.0076 

HALF  SPAN  -  48.1076     POISSON'S  RATIO  -  .3 

YIELD  STRENGTH  -  5970S5 .603738   ULTIMATE  TS  -  639475.293 

ELASTIC  MODULUS  -  1  .9994783497E+8   FLOW  STRESS  =  643280.448369 

CONVERSION  DATA 

FOR  LOAD: -1779. 288646  N/V  (  400  Lbs/V) 

FOR  COD:  .254  mm/V   (  .01  in/V) 

FOR  DEFLECTION :-1  .27  mm/V  (  .05  in/V> 

GIVEN  a-  .568  inches  THEN: 

Pmax  =  895.290182804  Lbs  FROM  ASTM  E813-81 

THE  FOLLOWING  INFORMATION  HAS  BEEN  OBTAINED: 

INITIAL  CRACK  LENGTH  -  14.3875332557  mm  FROM  MEASUREMENT  *  1 

WITH  95/C  CONFIDENCE  INTERVAL  FROM  14.4079726404  TO  14.3671845862  n 

USING  COD  COMPLIANCE  OF:  2 .81 520064535E-5  mm/N 

WITH  A  CORRELATION  OF:  .999323843963  FROM  189  DATA  POINTS. 

THE  FOLLOWING  INFORMATION  HAS  BEEN  OBTAINED: 

INITIAL  CRACK  LENGTH  -  14.4326755548  mm  FROM  MEASUREMENT  *   2 

WITH  957.  CONFIDENCE  INTERVAL  FROM  14.4442813846  TO  14.42109582  mm 

USING  COD  COMPLIANCE  OF:  2.84829954475E-5  mm/N 

WITH  A  CORRELATION  OF:  .99984474937  FROM  134  DATA  POINTS. 
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INITIAL  CRACK  LENGTH  =  14.4455712237  mm  FROM  MEASUREMENT  *  3 

WITH  95'/.  CONFIDENCE  INTERVAL  FROM  14.4565506791   TO  14.434615129  mm 

USING  COD  COMPLIANCE  OF:  2 .8578451 4764E-5  mm/N 

WITH  A  CORRELATION  OF:  .999816522357  FROM  176  DATA  POINTS. 

UNLOADING  *  1   

CRACK  LENGTH  -  14.4487507815  mm    DELTA  _a  -  .016075226716  mm 

957.  CONFIDENCE  INTERVAL  ON  CRACK  LENGTH  IS  FROM  14.4693514291   TO  14.4282322251 

TOTAL  #  DATA  POINTS  -  632    #  DATA  POINTS  THIS  UNLOADING  -  128 

CORRELATION  =»  .999529434569 

COD  COMPLIANCE  -  2.8755529381 7E-5  mm/N 

AREA  -  711.365042717  N»mm 

J  -  14.8049806501  KJ/m"2 

MAX  LOAD  LINE  DEFLECTION  -  .24034834255  mm 

CHANGE  IN  COD  -  .00109982  mm 

UNLOADING  #  2   

CRACK  LENGTH  -  14.4711513219  mm    DELTA  _a  ■  .0384757670308  mm 

952  CONFIDENCE  INTERVAL  ON  CRACK  LENGTH  IS  FROM  14.4927078135  TO  14.4496847244 

TOTAL  0   DATA  POINTS  -  947    #  DATA  POINTS  THIS  UNLOADING  ■  132 

CORRELATION  -  .999467111958 

COD  COMPLIANCE  -  2. 8930791 3363E-5  mm/N 

AREA  -  1693.69098454  N*mm 

J  ■  35.2970225265  KJ/m'2 

MAX  LOAD  LINE  DEFLECTION  -  .382347556072  mm 

CHANGE  IN  COD  -  .0008382  mm 

UNLOADING  *  3  

CRACK  LENGTH  -  14.4875164202  mm    DELTA  _a  -  .0548408654168  mm 

95%  CONFIDENCE  INTERVAL  ON  CRACK  LENGTH  IS  FROM  14.5143346851   TO  14.460837178 

TOTAL  0   DATA  POINTS  -  1300    *   DATA  POINTS  THIS  UNLOADING  -  157 

CORRELATION  -  .999017363069 

COD  COMPLIANCE  -  2.90942709989E-5  mm/N 

AREA  =  3033.0003923  N«mm 

J  -  63.28381022  KJ/m'2 

MAX  LOAD  LINE  DEFLECTION  -  .543433314464  mm 

CHANGE  IN  COD  -  .0009398  mm 

UNLOADING  *   4   

CRACK  LENGTH  -  14.5414951598  mm    DELTA  _a  -  .108819604992  mm 

957.  CONFIDENCE  INTERVAL  ON  CRACK  LENGTH  IS  FROM  14.5601316376  TO  14.5229259715 

TOTAL  *  DATA  POINTS  -  1670    #  DATA  POINTS  THIS  UNLOADING  -  109 

CORRELATION  -  .999675611283 

COD  COMPLIANCE  -  2.944429881 1 8E-5  mm/N 

AREA  -  5208.32927275  N*mm 

J  -  108.998439815  KJ/m'2 

MAX  LOAD  lINE  DEFLECTION  -  .782918890293  mm 
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CHANGE  IN  COD  =■  .0007874  mm 

UNLOADING  *   5  

CRACK  LENGTH  -  14.7407494376  mm    DELTA  _a  -  .308073382773  mm 

957.  CONFIDENCE  INTERVAL  ON  CRACK  LENGTH  IS  FROM  14.7587196422  TO  14.7228419911 

TOTAL  0   DATA  POINTS  -  2060    *   DATA  POINTS  THIS  UNLOADING  -  95 

CORRELATION  -  .999731105207 

COD  COMPLIANCE  -  3. 10248861 726E-5  mm/N 

AREA  -  8216.7400787  N»mra 

J  -  169.858689852  KJ/nT2 

MAX  LOAD  LINE  DEFLECTION  »  1.10455056629  mm 

CHANGE  IN  COD  -  .0008636  mm 

UNLOADING  *   6  

CRACK  LENGTH  -  14.8152586049  mm    DELTA  _a  »  .382583050059  mm 

957.  CONFIDENCE  INTERVAL  ON  CRACK  LENGTH  IS  FROM  14.8481840789  TO. ,1 4. 782543466 

TOTAL  #  DATA  POINTS  -  2230    #  DATA  POINTS  THIS  UNLOADING  =  138 

CORRELATION  -  .998642018089 

COD  COMPLIANCE  -  3  . 1 771921946 1E-5  mm/N 

AREA  -  9001.97373669  N-mm 

J  -  135.342172083  KJ/m'2 

MAX  LOAD  LINE  DEFLECTION  -  1.18989063161  mm 

CHANGE  IN  COD  -  .0009906  mm 

UNLOADING  *  7  

CRACK  LENGTH  -  14.8695206824  mm    DELTA  _a  -  .436845127588  mm 

957.  CONFIDENCE  INTERVAL  ON  CRACK  LENGTH  IS  FROM  14.3849566722  TO  14.8541311238 

TOTAL  *   DATA  POINTS  -  2531    #  DATA  POINTS  THIS  UNLOADING  =•  147 

CORRELATION  -  .999677856393 

COD  COMPLIANCE  -  3.20855793284E-5  mm/N 

AREA  -  9678.11102926  N»mm 

J  -  198.876552682  KJ/m~2 

MAX  LOAD  LINE  DEFLECTION  -  1.26297067275  mm 

CHANGE  IN  COD  -  .0004572  mm 

UNLOADING  #  8  

CRACK  LENGTH  -  15.0148445285  mm    DELTA  _a  =  .582163973634  mm 

957.  CONFIDENCE  INTERVAL  ON  CRACK  LENGTH  IS  FROM  15.0416778385  TO  14.9881516089 

TOTAL  *  DATA  POINTS  -  2761    #  DATA  POINTS  THIS  UNLOADING  -  97 

CORRELATION  -  .999364204197 

COD  COMPLIANCE  -  3.34673392487E-5  mm/N 

AREA  -  10586.6765231  N»mm 

J  =  215.448776655  KJ/m'2 

MAX  L.CAD  LIME  DEFLECTION  =  1.36396114491  mm 
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CHANGE  IN  COD  -  .0006604  mm 

UNLOADING  *   9  

CRACK  LENGTH  -  15.1216765395  mm    DELTA  _a  ■  .689000984704  mm 

95/.  CONFIDENCE  INTERVAL  ON  CRACK  LENGTH  IS  FROM  15.1358606171   TO  15.1075318754 

TOTAL  *   DATA  POINTS  -  3042    #  DATA  POINTS  THIS  UNLOADING  -  119 

CORRELATION  -  .99977636544 

COD  COMPLIANCE  -  3.43343078736E-5  mm/N 

AREA  -  11709.4655627  N»mm 

J  -  237.699077286  KJ/m~2 

MAX  LOAD  LINE  DEFLECTION  -  1.49020183752  mm 

CHANGE  IN  COD  -  .0005842  mm   " 

UNLOADING  *  10   

CRACK  LENGTH  -  15.2624877891  mm    DELTA  _a  -  .829812234314  mm 

957.  CONFIDENCE  INTERVAL  ON  CRACK  LENGTH  IS  FROM  15.2760292626  TO  15.2489823604 

TOTAL  *   DATA  POINTS  -  3394    *   DATA  POINTS  THIS  UNLOADING  -  167 

CORRELATION  -  .999700915055 

COD  COMPLIANCE  -  3.56785340653E-5  mm/N 

AREA  -  12983.4083914  N»mm 

J  -  262.394826944  KJ/m"2 

MAX  LOAD  LINE  DEFLECTION  -  1.63760345064  mm 

CHANGE  IN  COD  -  .0007366  mm 

UNLOADING  *  1 1   

CRACK  LENGTH  -  15.3993548961  mm    DELTA  _a  -  .966679341294  mm 

957.  CONFIDENCE  INTERVAL  ON  CRACK  LENGTH  IS  FROM  15.4107739959  TO  15.3879615278 

TOTAL  #  DATA  POINTS  -  3867    #  DATA  POINTS  THIS  UNLOADING  -  204 

CORRELATION  -  .999734089719 

COD  COMPLIANCE  -  3. 70351235759E-5  mm/N 

AREA  -  14557.9277896  N»mm 

J  -  294.038000741  KJ/nT2 

MAX  LOAD  LINE  DEFLECTION  -  1.82262379631  mm 

CHANGE  IN  COD  -  .0004318  mm 

UNLOADING  *    12  

CRACK  LENGTH  -  15.669999544  mm    DELTA  _a  -  1.2373239892  mm 

957.  CONFIDENCE  INTERVAL  ON  CRACK  LENGTH  IS  FROM  15 .5847004127   TO  15.6553416287 

TOTAL  *  DATA  POINTS  -  4267    *   DATA  POINTS  THIS  UNLOADING  -  154 

CORRELATION  -  .999652892868 

COD  COMPLIANCE  =  A.00050024785E-5  mm/M 

AREA  =  16330.4160222  N»mm 

J  =  325.442585526  KJ/m  2 

MAX  l-OAD  LINE  DEFLECTION  =  2.03825076401  mm 
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CHANGE  IN  COD  =■  .0013462  mm 

UNLOADING  #  13  

CRACK  LENGTH  -  15.8781418346  mm    DELTA  _a  -  1.44546627973  mm 

95'/.  CONFIDENCE  INTERVAL  ON  CRACK  LENGTH  IS  FROM  15.3901071812  TO  15.8662051558 

TOTAL  *   DATA  POINTS  -  4710    *  DATA  POINTS  THIS  UNLOADING  •  208 

CORRELATION  -  .999677232181 

COD  COMPLIANCE  -  4.24377123766E-5  mm/N 

AREA  -  17933.7183158  N»mm 

J  -  355.360006494  KJ/m~2 

MAX  LOAD  LINE  DEFLECTION  =•  2.24213294335  mm 

CHANGE  IN  COD  =■  .000508  mm 

UNLOADING  *  14  

CRACK  LENGTH  -  16.1357477216  mm    DELTA  _a  -  1.7030721668  mm 

957.  CONFIDENCE  INTERVAL  ON  CRACK  LENGTH  IS  FROM  16.1463233402  TO  16.1251947229 

TOTAL  *   DATA  POINTS  •  5206    *   DATA  POINTS  THIS  UNLOADING  ■  239 

CORRELATION  -  .999696024845 

COD  COMPLIANCE  -  4.57508725867E-5  mm/N 

AREA  »  19646.0705407  N»mm 

J  ■  385.620902249  KJ/m*2 

MAX  LOAD  LINE  DEFLECTION  =  2.47449867471  mm 

CHANGE  IN  COD  -  .0004826  mm 

UNLOADING  *  15  

CRACK  LENGTH  -  16.4193583405  mm    DELTA  _a  -  1.98668278568  mm 

952  CONFIDENCE  INTERVAL  ON  CRACK  LENGTH  IS  FROM  16.4364102869  TO  16.4022658462 

TOTAL  »   DATA  POINTS  -  5648    #  DATA  POINTS  THIS  UNLOADING  -  194 

CORRELATION  -  .999324290353 

COD  COMPLIANCE  -  4.99261 698493E-5  mm/N 

AREA  -  21264.009056  N»mm 

J  -  412.083508069  KJ/nT2 

MAX  LOAD  LINE  DEFLECTION  -  2.70365472703  mm 

CHANGE  IN  COD  -  .0005334  mm 

UNLOADING  *    16   

CRACK  LENGTH  -  16.9028626935  mm    DELTA  _a  -  2.47018713865  mm 

957.  CONFIDENCE  INTERVAL  ON  CRACK  LENGTH  IS  FROM  16.91 71 1 T3874  TO  16.388656585 

TOTAL  *   DATA  POINTS  -  6102    *  DATA  POINTS  THIS  UNLOADING  -  162 

CORRELATION  -  .999564028818 

COD  COMPLIANCE  -  5.80350094828E-5  mm/N 

AREA  ■  22950.1834392  N»mm 

J  =  428.162914109  KJ/m~2 

MAX  LOAD  LIME  DEFECTION  =  2.96684800605  mm 
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DATA  FROM  J  vs  DELTA.a  PLOT 

SLOPE  -  154.148013927    Y-INTERCEPT  -  123.933075018 
CORRELATION  -  .99667195294   FROM  9  DATA  POINTS. 


J1c  ■ 

142.05774949 

TEST  IS  OVER! 

NLOADING  0 

DELTA  a<mm> 

J<KJ/nT2> 

1 

.016075 

14.8 

2 

.038476 

35.3 

3 

.054841 

63.3 

4 

.108820 

109.0 

5 

.308074 

169.9 

6 

•  .382583 

185.3 

7 

.436845 

198:9 

8 

.532169 

215.4 

9 

.689001 

237.7 

10 

.829812 

262.4 

11 

.966679 

294.0 

12 

1.237324 

325.4 

13 

1 .445466 

355.4 

14 

1.703072 

385.6 

15 

1.986683 

412.1 

16 

2.470187 

428.2 

DATA  FROM  J  vs  DELTA_a  PLOT 

SLOPE  -  164.148013927    Y-INTERCEPT  -  123.923075018 
CORRELATION  -  .99667195294   FROM  9  DATA  POINTS. 

Jlc  -  142.05774949 


TEST  IS  OVER! 
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APPENDIX  I 
J_DATA  OUTPUT 

WHAT  DATA  FILE  DO  YOU  WISH  TO  READ? 
J-INTEGRAL  TEST  OF  HY80-5B   BY:  KURT 
DATE:  12-13-83       TIME:  10:35 

a0  =  14.4326755548   WITH  95%  CONDIFENCE  INTERVAL 
FROM  14.4442813846    TO  14.42109582 

FROM  COMPLIANCE  =  2.34829954475E-5 

CORRELATION-  .99984474937  FROM  176  DATA  POINTS 

DO  YOU  WANT  A  HARD  COPY  OF  THIS  INFORMATION? 


FROM  UNLOADING  #  1  RECORDED  128  DATA  POINTS 

COD  COMPLIANCE  =  2.8755529331 7E-5  mm/N   CORRELATION  =  .999529434569 

CRACK  LENGTH  -  14.4487507815  mm   DELTA_A  -  .016075226716  mm 

957.  CONFIDENCE  INTERVAL  FROM  14.4693514291  mm  TO  14.4282322251  mm 

AREA  =  711.365042717  N»mm   J  -  14.8049806501  KPa*m 

DATA  POINT  INDEXES  FROM  504  TO  632  WERE  USED. 

FROM  UNLOADING  *  2  RECORDED  132  DATA  POINTS 

COD  COMPLIANCE  -  2 .89307913363E-5  mm/N   CORRELATION  -  .999467111958 

CRACK  LENGTH  =  14.4711513219  mm   DELTA_A  =  .0384757670308  mm 

957.  CONFIDENCE  INTERVAL  FROM  14.4927078135  mm  TO  14.4496847244  mm 

AREA  -  1693.69098454  N-mm   J  -  35.2970225265  KPa»m 

DATA  POINT  INDEXES  FROM  815  TO  947  WERE  USED. 

FROM  UNLOADING  *   3  RECORDED  157  DATA  POINTS 

COD  COMPLIANCE  -  2.90942709989E-5  mm/N   CORRELATION  =  .999017363069 

CRACK  LENGTH  =  14.4875164202  mm   DELTA_A  =  .0548408654168  mm 

957.  CONFIDENCE  INTERVAL  FROM  14.5143346851  mm  TO  .  14 .4608371 78  mm 

AREA  =  3033.0003923  N»mm   J  -  63.28381022  KPa*m 

DATA  POINT  INDEXES  FROM  1143  TO  1300  WERE  USED. 

FROM  UNLOADING  *  4  RECORDED  109  DATA  POINTS 

COD  COMPLIANCE  =  2 .944429881 18E-5  mm/N   CORRELATION  -  .999675611283 

CRACK  LENGTH  -  14.5414951598  mm   DELTA_A  =  .108819604992  mm 

957.  CONFIDENCE  INTERVAL  FROM  14.5601316376  mm  TO  14.5229259715  mm 

AREA  =  5208.32927275  N»mm   J  =  108.998439815  KPa»m 

DATA  POINT  INDEXES  FROM  1561  TO  1670  WERE  USED. 

FROM  UNLOADING  *   5  RECORDED  95  DATA  POINTS 

COD  COMPLIANCE  =  3 . 1 0248861 726E-5  mm/N   CORRELATION  =  .999731105207 
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CRACK  LENGTH  =  14.7407494376  mm   DELTA_A  =  .308073882773  mm 
95*/.  CONFIDENCE  INTERVAL  FROM  14.7587196422  mm  TO  14.7228413911  mm 
AREA  =  32! 6. 7400787  N-mm   J  =  169.858689852  KPa*m 
DATA  POINT  INDEXES  FROM  1965  TO  2060  HERE  USED. 

FROM  UNLOADING  *   6  RECORDED  138  DATA  POINTS 

COD  COMPLIANCE  -  3. 1 771921 9461 E-5  mm/N   CORRELATION  =  .998642018089 

CRACK  LENGTH  =  14.8152586049  mm   DELTA_A  «  .382583050059  mm 

95%  CONFIDENCE  INTERVAL  FROM  14.8481840789  mm  TO  14.782543466  mm 

AREA  =  9001.97373669  N»mm   J  -  185.342172083  KPa*m 

DATA  POINT  INDEXES  FROM  2142  TO  2280  HERE  USED. 

FROM  UNLOADING  *  7  RECORDED  147  DATA  POINTS 

COD  COMPLIANCE  =  3.20855798284E-5  mm/N   CORRELATION  =  .999677856393 

CRACK  LENGTH  =  14.3695206824  mm   DELTA_A  =  .436845127588  mm 

957.  CONFIDENCE  INTERVAL  FROM  14.8849566722  mm  TO  14.8541311238  mm 

AREA  =  9678. lit 02926  N*mm   J  =  198.876552682  KPa*m 

DATA  POINT  INDEXES  FROM  2384  TO  2531   HERE  USED. 

FROM  UNLOADING  *  8  RECORDED  97  DATA  POINTS 

COD  COMPLIANCE  -  3.34673392487E-5  mm/N   CORRELATION  =  .999364204197 

CRACK  LENGTH  -  15.0148445285  mm   DELTA_A  =  .582168973634  mm 

95%  CONFIDENCE  INTERVAL  FROM  15.0416778385  mm  TO  14.9881516089  mm 

AREA  =  10586.6765231  N*mm   J  =  215.448776655  KPa»m 

DATA  POINT  INDEXES  FROM  2664  TO  2761   HERE  USED. 

FROM  UNLOADING  *  9  RECORDED  119  DATA  POINTS 

COD  COMPLIANCE  -  3 .43343078736E-5  mm/N   CORRELATION  =  .99977636544 

CRACK  LENGTH  -  15.1216765395  mm   DELTA_A  -  .689000984704  mm 

95%  CONFIDENCE  INTERVAL  FROM  15.1358606171  mm  TO  15.1075318754  mm 

AREA  -  11709.4655627  N»mm   J  -  237.699077286  KPa-m 

DATA  POINT  INDEXES  FROM  2923  TO  3042  HERE  USED. 

FROM  UNLOADING  *  10  RECORDED  167  DATA  POINTS 

COD  COMPLIANCE  =  3.56785340653E-5  mm/N   CORRELATION  -  .999700915055 

CRACK  LENGTH  -  15.2624877891  mm   DELTA_A  =  .829812234314  mm 

957.  CONFIDENCE  INTERVAL  FROM  15.2760292626  mm  TO  15.2489823604  ram 

AREA  =  12983.4083914  N-mm   J  =  262.394826944  KPa*m 

DATA  POINT  INDEXES  FROM  3227  TO  3394  HERE  USED. 

FROM  UNLOADING  *   11  RECORDED  204  DATA  POINTS 

COD  COMPLIANCE  =  3.70351235759E-5  mm/N   CORRELATION  =  .999734089719 

CRACK  LENGTH  =  15.3993548961  mm   DELTA_A  =  .966679341294  mm 

957.  CONFIDENCE  INTERVAL  FROM  15.4107739959  mm  TO  15.3879615278  mm 

AREA  =  14557.9277896  N*mm   J  =  294.038000741  KPa*m 

DATA  POINT  INDEXES  FROM  3663  TO  3867  HERE  USED. 

FROM  UNLOADING  *    12  RECORDED  154  DATA  POINTS 

COD  COMPLIANCE  -  4.00050024785E-5  mm/N   CORRELATION  =  .999652892868 

CRACK  LENGTH  =  15.669999544  mm   DELTA_A  »  1.2373239892  mm 

957.  CONFIDENCE  INTERVAL  FROM  15.6847004127  mm  TO  15.6553416287  mm 

AREA  =  16330.4160222  M»mm   J  =  325.442585526  KPa*m 

DATA  POINT  INDEXES  FROM  4113  TO  4267   HERE  USED. 
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FROM  UNLOADING  *  13  RECORDED  208  DATA  POINTS 

COD  COMPLIANCE  =  4 .243771 23766E-5  mm/N   CORRELATION  =  .399677232131 

CRACK  LENGTH  =  15.8781418346  mm   DELTA_A  =  1.44546627973  mm 

952  CONFIDENCE  INTERVAL  FROM  15.8901071812  mm  TO  15.8662051558  mm 

AREA  =  17933.7183158  N«mm   J  =  355.360006494  KPa*m 

DATA  POINT  INDEXES  FROM  4502  TO  4710  WERE  USED. 

FROM  UNLOADING  #  14  RECORDED  239  DATA  POINTS 

COD  COMPLIANCE  =  4 .57508725867E-5  mm/N   CORRELATION  =  .999696024845 

CRACK  LENGTH  =  16.1357477216  mm   DELTA_A  =  1.7030721668  mm 

957.  CONFIDENCE  INTERVAL  FROM  16.1463233402  mm  TO  16.1251947229  mm 

AREA  =  19646.0705407  N»mm   J  =  385.620902249  KPa*m 

DATA  POINT  INDEXES  FROM  4967  TO  5206  WERE  USED. 

FROM  UNLOADING  #  15  RECORDED  194  DATA  POINTS 

COD  COMPLIANCE  =  4 .99261 698493E-5  mm/N   CORRELATION  =  .999324290353 

CRACK  LENGTH  =  16.4193583405  mm   DELTA_A  =  1.98668278568  mm 

95'/.  CONFIDENCE  INTERVAL  FROM  16.4364102869  mm  TO  16.4023658462  mm 

AREA  =  21264.009056  N»mm   J  =  412.083508069  KPa*m 

DATA  POINT  INDEXES  FROM  5454  TO  5648  WERE  USED. 

FROM  UNLOADING  *  16  RECORDED  162  DATA  POINTS 

COD  COMPLIANCE  -  5.80350094828E-5  mm/N   CORRELATION  =  .999564028818 

CRACK  LENGTH  =  16.9028626935  mm   DELTA_A  =  2.47018713865  mm 

957.  CONFIDENCE  INTERVAL  FROM  16.9171113874  mm  TO  16.888656585  mm 

AREA  =  22950.1834392  N»mm   J  =  428.162914109  KPa*m 

DATA  POINT  INDEXES  FROM  5940  TO  6102  WERE  USED. 
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DATA   FROM   all   r-"0R:HY80-5B        BY:    KURT 
TESTED   ON:    12-13-83      TIME:    10:35 


# 

LOAD  <N> 

DEFLECTION  <mm> 

COD  (mm) 

1 

2262.5 

.079424 

.060576 

-> 

2261.7 

.074706 

.060576 

■J 

2265.2 

.076740 

.060637 

4 

2277.3 

.079003 

.06 1003 

5 

2350.8 

.093469 

.063533 

6 

2409.0 

.094978 

.065014 

7 

2442.3 

.095827 

.065395 

8 

2463.4 

.096491 

.066429 

9 

2473.9 

.096963 

.066954 

JO 

2496.2 

.097378 

.067230 

11 

2505.2 

.097666 

.067627 

12 

2515.4 

.098165 

.067772 

13 

2526.1 

.093554 

.063207 

14 

2539.  S 

.099114 

.063511 

15 

2554.0 

.099864 

.068969 

16 

2570.4 

.100177 

.06341 1 

17 

2595.3 

.101411 

.070231 

18 

2624.1 

.102889 

.071107 

13 

2653.4 

.103775 

.072045 

20 

2683.5 

.104678 

.072903 

21 

2722.7 

.106054 

.073863 

22 

2746.9 

.106721 

.074519 

23 

2764.1 

.107052 

.075024 

24 

2779.2 

.108077 

.075443 

25 

2799.9 

.108725 

.075956 

26 

2820.2 

.109596 

.076655 

27 

2837.3 

.110482 

.077142 

29 

2854.0 

.110667 

.077617 

29 

2873.4 

.112035 

.073313 

30 

2939.6 

.114805 

.030284 

31 

2984.3 

.115696 

.081471 

J  J 


1! 


J 


6085 

5686.4 

2.954239 

2.460422 

6086 

5692.7 

2.954786 

2.460701 

6087 

5699.1 

2.955204 

2.461203 

6088 

5704.9 

2.955629 

2.461616 

6089 

5716.5 

2.956736 

2.462454 

6090 

5731  .8 

2.957790 

2.463394 

6091 

5747.5 

2.959093 

2.464359 

6092 

5761.2 

2.959789 

2.465172 

6093 

5768.6 

2.960445 

2.465730 

6094 

5777.0 

2.961089 

2.466442 

6095 

5785.5 

2.961730 

2.467051 

6096 

5792.8 

2.962262 

2.467534 

6097 

5803.0 

2.963643 

2.468321 

6098 

5812.6 

2.964142 

2.469007 

6099 

5818.6 

2.364692 

2.469540 

6100 

582^.1 

2.965503 

2.470175 

6101 

5830.0 

2.365801 

2.470582 

6102 

5836.3 

2.366848 

2.471293 
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